Evaluation of a bounded high order upwind scheme for 3D incompressible free surface flow computations

fit the context of normalized variable formulation (NVF) of Leonard and total variation diminishing (TVD) constraints of Harten. this paper presents an extension of it previous work by the authors for solving unsteady incompressible flow problems. The main contributions of the paper are threefold. First, it presents the results of the development and implementation of a bounded high order upwind adaptative QUICKEST scheme in the 3D robust code (Freeflow), for the numerical solution of the full incompressible Navier-Stokes equations. Second, it reports numerical simulation results for 1D hock tube problem, 2D impinging jet and 2D/3D broken clam flows. Furthermore, these results are compared with existing analytical and experimental data. And third, it presents the application of the numerical method for solving 3D free surface flow problems. (C) 2007 IMACS. Published by Elsevier B.V. All rights reserved,

Numerical study of the atmospheric flow over a coastal cliff

This work presents a numerical method suitable for the study of the development of internal boundary layers (IBL) and their characteristics for flows over various types of coastal cliffs. The IBL is an important meteorological occurrence for flows with surface roughness and topographical step changes. A two-dimensional flow program was used for this study. The governing equations were written using the vorticity-velocity formulation. The spatial derivatives were discretized by high-order compact finite differences schemes. The time integration was performed with a low storage fourth-order Runge-Kutta scheme. The coastal cliff (step) was specified through an immersed boundary method. The validation of the code was done by comparison of the results with experimental and observational data. The numerical simulations were carried out for different coastal cliff heights and inclinations. The results show that the predominant factors for the height of the IBL and its characteristics are the upstream velocity, and the height and form (inclination) of the coastal cliff. Copyright (C) 2010 John Wiley & Sons, Ltd.

Twofold adaptive partition of unity implicits

Partition of Unity Implicits (PUI) has been recently introduced for surface reconstruction from point clouds. In this work, we propose a PUI method that employs a set of well-observed solutions in order to produce geometrically pleasant results without requiring time consuming or mathematically overloaded computations. One feature of our technique is the use of multivariate orthogonal polynomials in the least-squares approximation, which allows the recursive refinement of the local fittings in terms of the degree of the polynomial. However, since the use of high-order approximations based only on the number of available points is not reliable, we introduce the concept of coverage domain. In addition, the method relies on the use of an algebraically defined triangulation to handle two important tasks in PUI: the spatial decomposition and an adaptive polygonization. As the spatial subdivision is based on tetrahedra, the generated mesh may present poorly-shaped triangles that are improved in this work by means a specific vertex displacement technique. Furthermore, we also address sharp features and raw data treatment. A further contribution is based on the PUI locality property that leads to an intuitive scheme for improving or repairing the surface by means of editing local functions.

Numerical solution of the PTT constitutive equation for unsteady three-dimensional free surface flows

This work deals with the development of a numerical technique for simulating three-dimensional viscoelastic free surface flows using the PTT (Phan-Thien-Tanner) nonlinear constitutive equation. In particular, we are interested in flows possessing moving free surfaces. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. The fluid is modelled by a Marker-and-Cell type method and an accurate representation of the fluid surface is employed. The full free surface stress conditions are considered. The PTT equation is solved by a high order method, which requires the calculation of the extra-stress tensor on the mesh contours. To validate the numerical technique developed in this work flow predictions for fully developed pipe flow are compared with an analytic solution from the literature. Then, results of complex free surface flows using the FIT equation such as the transient extrudate swell problem and a jet flowing onto a rigid plate are presented. An investigation of the effects of the parameters epsilon and xi on the extrudate swell and jet buckling problems is reported. (C) 2010 Elsevier B.V. All rights reserved.

Investigation by Combined Solid-State NMR and SAXS Methods of the Morphology and Domain Size in Polystyrene-b-Polyethylene Oxide-b-Polystyrene Triblock Copolymers

The microphase structure of a series of polystyrene-b-polyethylene oxide-b-polystyrene (SEOS) triblock copolymers with different compositions and molecular weights has been studied by solid-state NMR, DSC, wide and small angle X-ray scattering (WAXS and SAXS). WAXS and DSC measurements were used to detect the presence of crystalline domains of polyethyleneoxide (PEO) blocks at room temperature as a function of the copolymer chemical composition. Furthermore, DSC experiments allowed the determination of the melting temperatures of the crystalline part of the PEO blocks. SAXS measurements, performed above and below the melting temperature of the PEO blocks, revealed the formation of periodic structures, but the absence or the weakness of high order reflections peaks did not allow a clear assessment of the morphological structure of the copolymers. This information was inferred by combining the results obtained by SAXS and (1)H NMR spin diffusion experiments, which also provided an estimation of the size of the dispersed phases of the nanostructured copolymers. (C) 2009 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 48:55-64,2010

Three-dimensional, fully adaptive simulations of phase-field fluid models

We present an efficient numerical methodology for the 31) computation of incompressible multi-phase flows described by conservative phase-field models We focus here on the case of density matched fluids with different viscosity (Model H) The numerical method employs adaptive mesh refinements (AMR) in concert with an efficient semi-implicit time discretization strategy and a linear, multi-level multigrid to relax high order stability constraints and to capture the flow`s disparate scales at optimal cost. Only five linear solvers are needed per time-step. Moreover, all the adaptive methodology is constructed from scratch to allow a systematic investigation of the key aspects of AMR in a conservative, phase-field setting. We validate the method and demonstrate its capabilities and efficacy with important examples of drop deformation, Kelvin-Helmholtz instability, and flow-induced drop coalescence (C) 2010 Elsevier Inc. All rights reserved

Theoretical Characterization of Hydrogen Polyoxides: HOOH, HOOOH, HOOOOH, and HOOO

We have investigated the polyoxides HOOH, HOOOH, HOOOOH, and HOOO employing the CCSD(T) methodology, and the correlation consistent basis sets. For all molecules, we have computed fundamental vibrational frequencies, structural parameters, rotational constants, and rotation-vibration corrections. For HOOOH, we have obtained a good agreement between our results and microwave and infrared spectra measurements, although for the symmetric OO stretch some important differences were found. Heats of formation were computed using atomization energies, and our recommendation is as follows: Delta H degrees(f,298)(HOOOH) = -21.50 kcal/mol and Delta H degrees(f,298)(HOOOOH) = -10.61 kcal/mol. In the case of HOOO, to estimate the heat of formation, we have constructed three isodesmic reactions to cancel high order correlation effects. The results obtained confirmed that the latter effects are very important for HOOO. The new Delta H degrees(f,298)(HOOO) obtained is 5.5 kcal/mol. We have also calculated the zero-point energies of DO and DOOO to correct the experimental lower limit determined for the Delta H degrees(f,298)(HOOO). The Delta(Delta ZPE) decreases the binding energy of HOOO by 0.56 kcal/mol. Employing the latter value, the new experimental lower limit for Delta H degrees(f,298)(HOOO) is 3.07 kcal/mol, just 2.4 kcal/mol lower than our determination. We expect that the fundamental vibrational frequencies and rotational constants determined for HOOOOH and DOOOOD contribute to its identification in the gas phase. The vibrational spectrum of HOOOOH shows some overlapping with that of HOOOH thus indicating that one may encounter some difficulties in its characterization. We discuss the consequences of the thermochemical properties determined in this work, and suggest that the amount of HOOO present in the atmosphere is smaller than that proposed recently in this journal (J. Phys. Chem A 2007, 111, 4727).

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

High homologous gene conservation despite extreme autopolyploid redundancy in sugarcane

P>Modern sugarcane (Saccharum spp.) is the leading sugar crop and a primary energy crop. It has the highest level of `vertical` redundancy (2n = 12x = 120) of all polyploid plants studied to date. It was produced about a century ago through hybridization between two autopolyploid species, namely S. officinarum and S. spontaneum. In order to investigate the genome dynamics in this highly polyploid context, we sequenced and compared seven hom(oe)ologous haplotypes (bacterial artificial chromosome clones). Our analysis revealed a high level of gene retention and colinearity, as well as high gene structure and sequence conservation, with an average sequence divergence of 4% for exons. Remarkably, all of the hom(oe)ologous genes were predicted as being functional (except for one gene fragment) and showed signs of evolving under purifying selection, with the exception of genes within segmental duplications. By contrast, transposable elements displayed a general absence of colinearity among hom(oe)ologous haplotypes and appeared to have undergone dynamic expansion in Saccharum, compared with sorghum, its close relative in the Andropogonea tribe. These results reinforce the general trend emerging from recent studies indicating the diverse and nuanced effect of polyploidy on genome dynamics.

The effects of UV radiation on photosynthesis estimated as chlorophyll fluorescence in Zygnemopsis decussata (Chlorophyta) growing in a high mountain lake (Sierra Nevada, Southern Spain)

The effect of increased UV radiation on photosynthesis estimated as in vivo chlorophyll fluorescence i.e. optimal quantum yield (F(v)/F(m)) and electron transport rate (ETR) in the green filamentous alga Zygnemopsis decussata (Streptophyta, Zygnematales) growing in the high mountain lake ""La Caldera"" (Sierra Nevada, Spain) at 3050 m altitude was evaluated. Two sets of in situ experiments were conducted: (1) On July 2006, F(v)/F(m) was measured throughout the day at different depths (0.1, 0.25, 0.5 and 1 m) and in the afternoon. ETR and phenolic compounds were determined. In addition, in order to analyze the effect of UV radiation, F(v)/F(m) was determined in algae incubated for 3 days at 0.5m under three different light treatments: PAR+UVA+UVB (PAB). PAR+UVA (PA) and PAR (P). (2) On August 2007, F(v)/F(m) was determined under PAB, PA and P treatments and desiccation/rehydration conditions. F(v)/F(m) decreased in algae growing in surface waters (0.1 m) but also at 1 m depth compared to that at 0.5 in depth. The decrease of F(v)/F(m) at noon due to photoinhibition was small (less than 10%) except in algae growing at 1 m depth (44%). The maximal electron transport rate was 3.5-5 times higher in algae growing at 0.25-0.5 m respectively than that at 0.1 and 1 m depth. These results are related to the accumulation of phenolic compounds: i.e. the algae at 0.25-0.5 in presentedrespectively about a 3-5 times higher concentration of phenolic compounds than that of algae at 0.1-1 m depth. The protection mechanisms seem to be stimulated by UVB radiation, since F(v)/F(m) was higher in the presence of UVB (PAB treatment) compared to PA or P treatments. UVA exerts the main photoinhibitory effect, not Only at midday, but also in the afternoon. UVB radiation also had a protective effect in algae grown under desiccation conditions for three days. During re-hydration, the rapid increase of F(v)/F(m) (after 1 h) was higher in the UVB-grown algae than in algae grown under UVA radiation. After 5 h. F(v)/F(m) values were similar in algae submitted to desiccation/rehydration under PAB and P treatments as they were in the control (submerged algae). The combined effect of desiccation and UVA produced the greatest decrease of photosynthesis in Z. decussata. Thifs UVB, in contrast to other species, may support the recovery process. Z. decussata can acclimate to severe stress, conditions in this high mountain lake by the photoprotection mechanism induced by UVB radiation through dynamic photoinhibition and the accumulation of phenolic compounds (UV screen and antioxidant substances).

Structural and Biochemical Characterization of Peroxiredoxin Q beta from Xylella fastidiosa CATALYTIC MECHANISM AND HIGH REACTIVITY

The phytopathogenic bacterium Xylella fastidiosa is the etiological agent of various plant diseases. To survive under oxidative stress imposed by the host, microorganisms express antioxidant proteins, including cysteine-based peroxidases named peroxiredoxins. This work is a comprehensive analysis of the catalysis performed by PrxQ from X. fastidiosa (XfPrxQ) that belongs to a peroxiredoxin class still poorly characterized and previously considered as moderately reactive toward hydroperoxides. Contrary to these assumptions, our competitive kinetics studies have shown that the second-order rate constants of the peroxidase reactions of XfPrxQ with hydrogen peroxide and peroxynitrite are in the order of 107 and 106 M(-1) s(-1), respectively, which are as fast as the most efficient peroxidases. The XfPrxQ disulfides were only slightly reducible by dithiothreitol; therefore, the identification of a thioredoxin system as the probable biological reductant of XfPrxQ was a relevant finding. We also showed by site-specific mutagenesis and mass spectrometry that an intramolecular disulfide bond between Cys-47 and Cys-83 is generated during the catalytic cycle. Furthermore, we elucidated the crystal structure of XfPrxQ C47S in which Ser-47 and Cys-83 lie similar to 12.3 angstrom apart. Therefore, significant conformational changes are required for disulfide bond formation. In fact, circular dichroism data indicated that there was a significant redox-dependent unfolding of alpha-helices, which is probably triggered by the peroxidatic cysteine oxidation. Finally, we proposed a model that takes data from this work as well data as from the literature into account.

Continuous and High-Level In Vivo Delivery of Endostatin From Recombinant Cells Encapsulated in TheraCyte (R) Immunoisolation Devices

Endostatin (ES) is a potent inhibitor of angiogenesis and tumor growth. Continuous ES delivery of ES improves the efficacy and potency of the antitumoral therapy. The TheraCyte (R) system is a polytetrafluoroethylene (PTFE) semipermeable membrane macroencapsulation system for implantation of genetically engineered cells specially designed for the in vivo delivery of therapeutic proteins, such as ES, which circumvents the problem of limited half-life and variation in circulating levels. In order to enable neovascularization at the tissues adjacent to the devices prior to ES secretion by the cells inside them, we designed a scheme in which empty TheraCyte (R) devices were preimplanted SC into immunodeficient mice. Only after healing (17 days later) were Chinese hamster ovary cells expressing ES injected into the preimplanted devices. In another model for device implantation, the cells expressing ES where loaded into the immunoisolation devices prior to implantation into the animals, and the TheraCyte (R) were then immediately implanted SC into the mice. Throughout the 2-month study, constant high ES levels of up to 3.7 mu g/ml were detected in the plasma of the mice preimplanted with the devices, while lower but also constant levels of ES (up to 2.1 mu g/ml plasma) were detected in the mice that had received devices preloaded with the ES-expressing cells. Immunohistochemistry using anti-ES antibody showed reaction within the device and outside it, demonstrating that ES, secreted by the confined recombinant cells, permeated through the membrane and reached the surrounding tissues.

Least square projection: A fast high-precision multidimensional projection technique and its application to document mapping

The problem of projecting multidimensional data into lower dimensions has been pursued by many researchers due to its potential application to data analyses of various kinds. This paper presents a novel multidimensional projection technique based on least square approximations. The approximations compute the coordinates of a set of projected points based on the coordinates of a reduced number of control points with defined geometry. We name the technique Least Square Projections ( LSP). From an initial projection of the control points, LSP defines the positioning of their neighboring points through a numerical solution that aims at preserving a similarity relationship between the points given by a metric in mD. In order to perform the projection, a small number of distance calculations are necessary, and no repositioning of the points is required to obtain a final solution with satisfactory precision. The results show the capability of the technique to form groups of points by degree of similarity in 2D. We illustrate that capability through its application to mapping collections of textual documents from varied sources, a strategic yet difficult application. LSP is faster and more accurate than other existing high-quality methods, particularly where it was mostly tested, that is, for mapping text sets.

MODIS aerosol optical depth Retrievals with high spatial resolution over an urban area using the critical reflectance

[1] The retrieval of aerosol optical depth (Ta) over land by satellite remote sensing is still a challenge when a high spatial resolution is required. This study presents a tool that uses satellite measurements to dynamically identify the aerosol optical model that best represents the optical properties of the aerosol present in the atmosphere. We use aerosol critical reflectance to identify the single scattering albedo of the aerosol layer. Two case studies show that the Sao Paulo region can have different aerosol properties and demonstrates how the dynamic methodology works to identify those differences to obtain a better T a retrieval. The methodology assigned the high single scattering albedo aerosol model (pi o( lambda = 0.55) = 0.90) to the case where the aerosol source was dominated by biomass burning and the lower pi(o) model (pi(o) (lambda = 0.55) = 0.85) to the case where the local urban aerosol had the dominant influence on the region, as expected. The dynamic methodology was applied using cloud-free data from 2002 to 2005 in order to retrieve Ta with Moderate Resolution Imaging Spectroradiometer ( MODIS). These results were compared with collocated data measured by AERONET in Sao Paulo. The comparison shows better results when the dynamic methodology using two aerosol optical models is applied (slope 1.06 +/- 0.08 offset 0.01 +/- 0.02 r(2) 0.6) than when a single and fixed aerosol model is used (slope 1.48 +/- 0.11 and offset - 0.03 +/- 0.03 r(2) 0.6). In conclusion the dynamical methodology is shown to work well with two aerosol models. Further studies are necessary to evaluate the methodology in other regions and under different conditions.

Fragment mass distributions in the fission of heavy nuclei by intermediate- and high-energy probes

Recent experiments have shown that the multimode approach for describing the fission process is compatible with the observed results. Asystematic analysis of the parameters obtained by fitting the fission-fragment mass distribution to the spontaneous and low-energy data has shown that the values for those parameters present a smooth dependence upon the nuclear mass number. In this work, a new methodology is introduced for studying fragment mass distributions through the multimode approach. It is shown that for fission induced by energetic probes (E > 30 MeV) the mass distribution of the fissioning nuclei produced during the intranuclear cascade and evaporation processes must be considered in order to have a realistic description of the fission process. The method is applied to study (208)Pb, (238)U, (239)Np and (241)Am fission induced by protons or photons.