A continuously differentiable upwinding scheme for the simulation of fluid flow problems

This paper deals with the numerical solution of complex fluid dynamics problems using a new bounded high resolution upwind scheme (called SDPUS-C1 henceforth), for convection term discretization. The scheme is based on TVD and CBC stability criteria and is implemented in the context of the finite volume/difference methodologies, either into the CLAWPACK software package for compressible flows or in the Freeflow simulation system for incompressible viscous flows. The performance of the proposed upwind non-oscillatory scheme is demonstrated by solving two-dimensional compressible flow problems, such as shock wave propagation and two-dimensional/axisymmetric incompressible moving free surface flows. The numerical results demonstrate that this new cell-interface reconstruction technique works very well in several practical applications. (C) 2012 Elsevier Inc. All rights reserved.

Simulation results and applications of an advection bounded scheme to practical flows

This paper reports experiments on the use of a recently introduced advection bounded upwinding scheme, namely TOPUS (Computers & Fluids 57 (2012) 208-224), for flows of practical interest. The numerical results are compared against analytical, numerical and experimental data and show good agreement with them. It is concluded that the TOPUS scheme is a competent, powerful and generic scheme for complex flow phenomena.

Semilinear functional difference equations with infinite delay

We obtain boundedness and asymptotic behavior of solutions for semilinear functional difference equations with infinite delay. Applications to Volterra difference equations with infinite delay are shown. (C) 2011 Elsevier Ltd. All rights reserved.

Flow Visualization with Quantified Spatial and Temporal Errors Using Edge Maps

Robust analysis of vector fields has been established as an important tool for deriving insights from the complex systems these fields model. Traditional analysis and visualization techniques rely primarily on computing streamlines through numerical integration. The inherent numerical errors of such approaches are usually ignored, leading to inconsistencies that cause unreliable visualizations and can ultimately prevent in-depth analysis. We propose a new representation for vector fields on surfaces that replaces numerical integration through triangles with maps from the triangle boundaries to themselves. This representation, called edge maps, permits a concise description of flow behaviors and is equivalent to computing all possible streamlines at a user defined error threshold. Independent of this error streamlines computed using edge maps are guaranteed to be consistent up to floating point precision, enabling the stable extraction of features such as the topological skeleton. Furthermore, our representation explicitly stores spatial and temporal errors which we use to produce more informative visualizations. This work describes the construction of edge maps, the error quantification, and a refinement procedure to adhere to a user defined error bound. Finally, we introduce new visualizations using the additional information provided by edge maps to indicate the uncertainty involved in computing streamlines and topological structures.

Influence of the Substrate and Precursor on the Magnetic and Magneto-transport Properties in Magnetite Films

We have investigated the magnetic and transport properties of nanoscaled Fe3O4 films obtained from Chemical Vapor Deposition (CVD) technique using [(FeFe2III)-Fe-II(OBut)(8)] and [Fe-2(III)(OBut)(6)] precursors. Samples were deposited on different substrates (i.e., MgO (001), MgAl2O4 (001) and Al2O3 (0001)) with thicknesses varying from 50 to 350 nm. Atomic Force Microscopy analysis indicated a granular nature of the samples, irrespective of the synthesis conditions (precursor and deposition temperature, T-pre) and substrate. Despite the similar morphology of the films, magnetic and transport properties were found to depend on the precursor used for deposition. Using [(FeFe2III)-Fe-II(OBut)(8)] as precursor resulted in lower resistivity, higher M-S and a sharper magnetization decrease at the Verwey transition (T-V). The temperature dependence of resistivity was found to depend on the precursor and T-pre. We found that the transport is dominated by the density of antiferromagnetic antiphase boundaries (AF-APB's) when [(FeFe2III)-Fe-II(OBut)(8)] precursor and T-pre = 363 K are used. On the other hand, grain boundary-scattering seems to be the main mechanism when [Fe-2(III)(OBut)(6)] is used. The Magnetoresistance (MR(H)) displayed an approximate linear behavior in the high field regime (H > 796 kA/m), with a maximum value at room-temperature of similar to 2-3 % for H = 1592 kA/m, irrespective from the transport mechanism.

Ionic Transport and Water Vapor Uptake of Ammonium Exchanged Perfluorosulfonic Acid Membranes

Nafion membranes series N117 doped with ammonium, at different cation fractions (H+/NH4+), were investigated for ionic transport and water vapor uptake, for several water activities and temperatures. Ammonium cations change both properties of the polymer in a similar manner. Membrane ionic conductivity and water vapor uptake (lambda) decrease as the ammonium concentration increases in the polymer. Ionic transport activation energies are calculated and the transport mechanism of ammonium ions in Nafion is discussed. (C) 2012 The Electrochemical Society. [DOI: 10.1149/2.040203jes] All rights reserved.

Study of the Influence of Localized Vibrational Modes in Charge Transport Properties at Nanoscale Systems.

In molecular and atomic devices the interaction between electrons and ionic vibrations has an important role in electronic transport. The electron-phonon coupling can cause the loss of the electron's phase coherence, the opening of new conductance channels and the suppression of purely elastic ones. From the technological viewpoint phonons might restrict the efficiency of electronic devices by energy dissipation, causing heating, power loss and instability. The state of the art in electron transport calculations consists in combining ab initio calculations via Density Functional Theory (DFT) with Non-Equilibrium Green's Function formalism (NEGF). In order to include electron-phonon interactions, one needs in principle to include a self-energy scattering term in the open system Hamiltonian which takes into account the effect of the phonons over the electrons and vice versa. Nevertheless this term could be obtained approximately by perturbative methods. In the First Born Approximation one considers only the first order terms of the electronic Green's function expansion. In the Self-Consistent Born Approximation, the interaction self-energy is calculated with the perturbed electronic Green's function in a self-consistent way. In this work we describe how to incorporate the electron-phonon interaction to the SMEAGOL program (Spin and Molecular Electronics in Atomically Generated Orbital Landscapes), an ab initio code for electronic transport based on the combination of DFT + NEGF. This provides a tool for calculating the transport properties of materials' specific system, particularly in molecular electronics. Preliminary results will be presented, showing the effects produced by considering the electron-phonon interaction in nanoscale devices.

Evaluation of an Experimental Gel Containing Euclea natalensis: An In Vitro Study

Objective. To evaluate the effect of an experimental gel containing Euclea natalensis extract on dentin permeability. Methods. Thirty-six dentin discs, 1-mm-thick. The discs were prepared from the coronal dentin of extracted human third molars that were divided into 3 groups (n = 10). The dentin discs in each group were treated with the groups following experimental materials: (FG): 1.23% fluoride gel, pH 4.1; (EG): Euclea natalensis extract gel, pH 4.1; (CG): control gel, pH 4.1. The gels were applied to the occlusal slide of the dentin under the following conditions: after 37% phosphoric acid and before 6% citric acid. The hydraulic conductance (HC) of each condition was determined four times using a fluid flow apparatus (Flodec). The data were analyzed using Two-way ANOVA and Tukey's test (P < 0.05). Results. The greatest mean reduction in HC was produced in group EG dentin discs (61.2%; P < 0.05). Even after acid challenge with 6% citric acid the great reduction occurred in group EG (66.0%; P < 0.05) than other groups (CG-77.1%, FG-90.8%). Conclusion. E. natalensis gel not only reduced dentin permeability, but also resisted posttreatment citric acid challenge without changing its permeability. Further research has to confirm this promising result in the clinical situation.

Numerical simulation of drop impact and jet buckling problems using the eXtended Pom-Pom model

This work presents numerical simulations of two fluid flow problems involving moving free surfaces: the impacting drop and fluid jet buckling. The viscoelastic model used in these simulations is the eXtended Pom-Pom (XPP) model. To validate the code, numerical predictions of the drop impact problem for Newtonian and Oldroyd-B fluids are presented and compared with other methods. In particular, a benchmark on numerical simulations for a XPP drop impacting on a rigid plate is performed for a wide range of the relevant parameters. Finally, to provide an additional application of free surface flows of XPP fluids, the viscous jet buckling problem is simulated and discussed. (C) 2011 Elsevier B.V. All rights reserved.

Prediction of the Temperature Distribution of Partially Submersed Umbilical Cables

The objective of this work is to predict the temperature distribution of partially submersed umbilical cables under different operating and environmental conditions. The commercial code Fluent (R) was used to simulate the heat transfer and the air fluid flow of part of a vertical umbilical cable near the air-water interface. A free-convective three-dimensional turbulent flow in open-ended vertical annuli was solved. The influence of parameters such as the heat dissipating rate, wind velocity, air temperature and solar radiation was analyzed. The influence of the presence of a radiation shield consisting of a partially submersed cylindrical steel tube was also considered. The air flow and the buoyancy-driven convective heat transfer in the annular region between the steel tube and the umbilical cable were calculated using the standard k-epsilon turbulence model. The radiative heat transfer between the umbilical external surface and the radiation shield was calculated using the Discrete Ordinates model. The results indicate that the influence of a hot environment and intense solar radiation may affect the umbilical cable performance in its dry portion.

An in-depth analysis of the HIV-1/AIDS dynamics by comprehensive mathematical modeling

This work presents major results from a novel dynamic model intended to deterministically represent the complex relation between HIV-1 and the human immune system. The novel structure of the model extends previous work by representing different host anatomic compartments under a more in-depth cellular and molecular immunological phenomenology. Recently identified mechanisms related to HIV-1 infection as well as other well known relevant mechanisms typically ignored in mathematical models of HIV-1 pathogenesis and immunology, such as cell-cell transmission, are also addressed. (C) 2011 Elsevier Ltd. All rights reserved.

BEM modeling of saturated porous media susceptible to damage

This paper deals with the numerical analysis of saturated porous media, taking into account the damage phenomena on the solid skeleton. The porous media is taken into poro-elastic framework, in full-saturated condition, based on Biot's Theory. A scalar damage model is assumed for this analysis. An implicit boundary element method (BEM) formulation, based on time-independent fundamental solutions, is developed and implemented to couple the fluid flow and two-dimensional elastostatic problems. The integration over boundary elements is evaluated using a numerical Gauss procedure. A semi-analytical scheme for the case of triangular domain cells is followed to carry out the relevant domain integrals. The non-linear problem is solved by a Newton-Raphson procedure. Numerical examples are presented, in order to validate the implemented formulation and to illustrate its efficacy. (C) 2011 Elsevier Ltd. All rights reserved.

Linking origin of the electric field-assisted beta-PbF2 crystallization in lead oxyfluoroborate glasses below T (g) to simultaneous cathode/anode-compensated electrochemical reactions

Full validation of the electrochemical mechanisms so far postulated as driving force of electric field-assisted non-spontaneous crystallization development in given glasses has suffered experimental restrictions. In this work, we looked into origin of this phenomenon in lead oxyfluoroborate glasses, resulting in beta-PbF2 growth even below the corresponding glass transition temperatures, through achieving a systematic study of not only Pt,Ag/Glass/Ag,Pt- but also Pt,Ag/Glass/YSZ:PbF2/Ag,Pt-type cells, where YSZ:PbF2 represents a two-phase system (formed by Y2O3-doped ZrO2 and PbF2). It is demonstrated that crystallization induction in these glasses involves Pb2+ ions reduction at the cathode, the phenomenon being, however, confirmed only when the F- ions were simultaneously also able to reach the anode for oxidation, after assuring either a direct glass-anode contact or percolation pathways for free fluoride migration across the YSZ:PbF2 mixtures. A further support of this account is that the electrochemically induced beta-PbF2 phase crystallizes showing ramified-like microstructure morphology that arises, accordingly, from development of electroconvective diffusion processes under electric field action.

Efficiency of extrinsic and intrinsic charge-carrier photogeneration processes obtained from the steady-state photocurrent action spectra of poly(p-phenylene vinylene) derivatives

The efficiency of the charge-carrier photogeneration processes in poly(2,5-bis(3',7'-dimethyl-octyloxy)-1,4-phenylene vinylene) (OC(1)OC10-PPV) has been analyzed by the spectral response of the photocurrent of devices in ITO/polymer/Al structures. The symbatic response of the photocurrent action spectra of the OC1OC10-PPV devices, obtained for light-excitation through the ITO electrode and for forward bias, has been fitted using a phenomenological model which considers that the predominant transport mechanism under external applied electric field is the drift of photogenerated charge-carriers, neglecting charge-carrier diffusion. The proposed model takes into account that charge-carrier photogeneration occurs via intermediate stages of bounded pairs (excitonic states), followed by dissociation processes. Such processes result in two different contributions to the photoconductivity: The first one, associated to direct creation of unbound polaron pairs due to intrinsic photoionization; and the second one is associated to secondary processes like extrinsic photoinjection at the metallic electrodes. The results obtained from the model have shown that the intrinsic component of the photoconductivity at higher excitation energies has a considerably higher efficiency than the extrinsic one, suggesting a dependence on the photon energy for the efficiency of the photogeneration process.

A comparison for a multiscale study of structural lineaments in southern Brazil: LANDSAT-7 ETM+ and shaded relief images from SRTM3-DEM

This paper presents a comparison of descriptive statistics obtained for brittle structural lineaments extracted manually from LANDSAT images and shaded relief images from SRTM 3 DEM at 1:100, 000 and 1:500, 000 scales. The selected area is located in the southern of Brazil and comprises Precambrian rocks and stratigraphic units of the Paraná Basin. The application of this methodology shows that the visual interpretation depends on the kind of remote sensing image. The resulting descriptive statistics obtained for lineaments extracted from the images do not follow the same pattern according to the scale adopted. The main direction obtained for Proterozoic rocks using both image types at a 1:500, 000 scale are close to NS±10, whereas at a 1:100, 000 scale N45E was obtained for shaded relief images from SRTM 3 DEM and N10W for LANDSAT images. The Paleozoic sediments yielded the best results for the different images and scales (N50W). On the other hand, the Mesozoic igneous rocks showed greatest differences, the shaded relief images from SRTM 3 DEM images highlighting NE structures and the LANDSAT images highlighting NW structures. The accumulated frequency demonstrated high similarity between products for each image type no matter the scale, indicating that they can be used in multiscale studies. Conversely, major differences were found when comparing data obtained using shaded relief images from SRTM 3 DEM and Landsat images at a 1:100, 000 scale.