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,

Nanostructured organic layers via polymer demixing for interface-enhanced photovoltaic cells

Significant progress is being made in the photovoltaic energy conversion using organic semiconducting materials. One of the focuses of attention is the morphology of the donor-acceptor heterojunction at the nanometer scale, to ensure efficient charge generation and loss-free charge transport at the same time. Here, we present a method for the controlled, sequential design of a bilayer polymer cell architecture that consists of a large interface area with connecting paths to the respective electrodes for both materials. We used the surface-directed demixing of a donor conjugated/guest polymer blend during spin coating to produce a nanostructured interface, which was, after removal of the guest with a selective solvent, covered with an acceptor layer. With use of a donor poly(p-phenylenevinylene) derivative and the acceptor C-60 fullerene, this resulted in much-improved device performance, with external power efficiencies more than 3 times higher than those reported for that particular material combination so far.

Numerical assessment of mass conservation on a MAC-type method for viscoelastic free surface flows

A numerical study of mass conservation of MAC-type methods is presented, for viscoelastic free-surface flows. We use an implicit formulation which allows for greater time steps, and therefore time marching schemes for advecting the free surface marker particles have to be accurate in order to preserve the good mass conservation properties of this methodology. We then present an improvement by using a Runge-Kutta scheme coupled with a local linear extrapolation on the free surface. A thorough study of the viscoelastic impacting drop problem, for both Oldroyd-B and XPP fluid models, is presented, investigating the influence of timestep, grid spacing and other model parameters to the overall mass conservation of the method. Furthermore, an unsteady fountain flow is also simulated to illustrate the low mass conservation error obtained.

Desenvolvimento de interface amigável de rotina em linguagem c++ para o processamento aritmético e estatístico de dados censitários

A routine was developed in C++ for the processing of social and environmental census data acquired by the Brazilian Institute of Geography and Statistics (IBGE). The routine employs a simple graphical environment. The data generated are presented in a tabular format, which facilitates a broad and objective view of the values, and provides a convenient means of querying the database. The source code used to develop the routine permits updates and changes, as required by the user. Statistical and mathematical analysis enables the generation of social and environmental indicators, together with quantitative and qualitative classification of the socio-environmental quality of the region analyzed. As an example, the routine was applied using census data for the city of Sorocaba (São Paulo State, Brazil), including conditions of household occupation, water supply, sanitation, level of education, income, and other factors. It is envisaged that the proposed analytical model will assist professionals from different fields of research and teaching to develop urban planning and management strategies.

Estudo da interface titânio/porcelana, com três porcelanas de ultra baixa fusão, empregando-se microscopia eletrônica de varredura e dilatometria

Pós-graduação em Ciências Odontológicas - FOAR

Avaliação da interface dentina/agentes protetores em restaurações com resina composta, através de microscopia eletrônica de varredura

Pós-graduação em Ciências Odontológicas - FOAR

Characterization of the Al-3wt.%Si alloy in unsteady-state horizontal directional solidification

The main purpose of this paper is to investigate both the columnar to equiaxed transition and primary dendritic arm spacings of Al-3wt.%Si alloy during the horizontal directional solidification. The transient heat transfer coefficient at the metal-mold interface is calculated based on comparisons between the experimental thermal profiles in castings and the simulations provided by a finite difference heat flow program. Simulated curve of the interfacial heat transfer coefficient was used in another numerical solidification model to determine theoretical values of tip growth rates, cooling rates and thermal gradients that are associated with both columnar to equiaxed transition and primary dendritic arm spacings. A good agreement was observed between the experimental values of these thermal variables and those numerically simulated for the alloy examined. A comparative analysis is carried out between the experimental data of this work and theoretical models from the literature that have been proposed to predict the primary dendritic spacings. In this context, this study may contribute to the understanding of how to manage solidification operational parameters aiming at designing the microstructure of Al-Si alloys.

Application of the PTT Model to Axisymmetric Free Surface Flows

This work is concerned with numerical simulation of axisymmetric viscoelastic free surface flows using the Phan-Thien-Tanner (PTT) constitutive equation. A finite difference technique for solving the governing equations for unsteady incompressible flows written in Cylindrical coordinates on a staggered grid is described. 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 applied. The numerical method is verified by comparing numerical predictions of fully developed flow in a pipe with the corresponding analytic solutions. To demonstrate that the numerical method can simulate axisymmetric free surface flows governed by the PTT model, numerical results of the flow evolution of a drop impacting on a rigid dry plate are presented. In these simulations, the rheological effects of the parameters epsilon and xi are investigated.

Serological evidence of Infection by Leishmania (Leishmania) infantum (Synonym: Leishmania (Leishmania) chagasi) in free-ranging wild mammals in a nonendemic region of the state of São Paulo, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Avaliação do selamento de ápices radiculares retrobturados ou recobertos com diferentes materiais mediante infiltração marginal por corante e análise da interface em microscopia eletrônica de varredura

The purpose of this study was to evaluate the apical sealing ability and the marginal adaptation of five dental materials used in retrofillings or applied to the bevelled root surface. One hundred and forty extracted single rooted human teeth were used, divided into seven groups of twenty each. ln the first, second, third and fourth groups, the teeth were apicoectomized, submitted to cavity preparations and retrofilled with one of the following materiais: zinc free silver amalgam, a dentin bonding system plus composite resin, a glass ionomer cement ora compomer. In the fifth, sixth and seventh groups, the teeth were apicoectomized and capped on the bevelled root surface with one of the following materiais: a dentin bonding system plus composite resin, a glass ionomer cement or a compomer. Two specimens of each experimental group were evaluated for the marginal adaptation using scanning electron microscopy. The remaining 126 specimens were immersed in 2% methylene blue dye, stored for one week at 37ºC and the infiltration was evaluated with a stereomicroscope. The results showed that the glass ionomer cement presented the lowest values of marginal infiltration when used as retrofilling material, with a significant statistical difference when compared with the others tested materials. When used as apical capping, the glass ionomer cement and the compomer were equivalent and significantly better than the dentin bonding system plus composite resin. Using scanning electron microscopy, all the materials showed some slight adjustment problem. ln the retrofilling, the smallest marginal gaps were observed with the compomer and the dentin bonding system plus composite resin, while the largest were observed with the glass ionomer cement and zinc free silver amalgam. ln the apical capping, the smallest marginal gaps were observed with the compomer and the dentin bonding system plus composite resin and ...

Interface tailoring for adhesion enhancement of diamond-like carbon thin films

We have explored the suitability and characteristics of interface tailoring as a tool for enhancing the adhesion of hydrogen-free diamond-like carbon (DLC) thin films to silicon substrates. DLC films were deposited on silicon with and without application of an initial high energy carbon ion bombardment phase that formed a broad Si-C interface of gradually changing Si:C composition. The interface depth profile was calculated using the TRIDYN simulation program, revealing a gradient of carbon concentration including a region with the stoichiometry of silicon carbide. DLC films on silicon, with and without interface tailoring, were characterized using Raman spectroscopy, scanning electron microscopy, atomic force microscopy and scratch tests. The Raman spectroscopy results indicated sp3-type carbon bonding content of up to 80%. Formation of a broadened Si:C interface as formed here significantly enhances the adhesion of DLC films to the underlying silicon substrate. (C) 2012 Elsevier B.V. All rights reserved.

Microemulsion: prediction of the phase diagram with a modified Helfrich free energy

Microemulsions are thermodynamically stable, macroscopically homogeneous but microscopically heterogeneous, mixtures of water and oil stabilised by surfactant molecules. They have unique properties like ultralow interfacial tension, large interfacial area and the ability to solubilise other immiscible liquids. Depending on the temperature and concentration, non-ionic surfactants self assemble to micelles, flat lamellar, hexagonal and sponge like bicontinuous morphologies. Microemulsions have three different macroscopic phases (a) 1phase- microemulsion (isotropic), (b) 2phase-microemulsion coexisting with either expelled water or oil and (c) 3phase- microemulsion coexisting with expelled water and oil.rnrnOne of the most important fundamental questions in this field is the relation between the properties of the surfactant monolayer at water-oil interface and those of microemulsion. This monolayer forms an extended interface whose local curvature determines the structure of the microemulsion. The main part of my thesis deals with the quantitative measurements of the temperature induced phase transitions of water-oil-nonionic microemulsions and their interpretation using the temperature dependent spontaneous curvature [c0(T)] of the surfactant monolayer. In a 1phase- region, conservation of the components determines the droplet (domain) size (R) whereas in 2phase-region, it is determined by the temperature dependence of c0(T). The Helfrich bending free energy density includes the dependence of the droplet size on c0(T) as

Computational simulations of latent heat thermal energy storage systems - with innovative and first-principles based simulation for the underlying unsteady melting (and solidification) processes

This thesis develops an effective modeling and simulation procedure for a specific thermal energy storage system commonly used and recommended for various applications (such as an auxiliary energy storage system for solar heating based Rankine cycle power plant). This thermal energy storage system transfers heat from a hot fluid (termed as heat transfer fluid - HTF) flowing in a tube to the surrounding phase change material (PCM). Through unsteady melting or freezing process, the PCM absorbs or releases thermal energy in the form of latent heat. Both scientific and engineering information is obtained by the proposed first-principle based modeling and simulation procedure. On the scientific side, the approach accurately tracks the moving melt-front (modeled as a sharp liquid-solid interface) and provides all necessary information about the time-varying heat-flow rates, temperature profiles, stored thermal energy, etc. On the engineering side, the proposed approach is unique in its ability to accurately solve – both individually and collectively – all the conjugate unsteady heat transfer problems for each of the components of the thermal storage system. This yields critical system level information on the various time-varying effectiveness and efficiency parameters for the thermal storage system.

Ragu: a free tool for the analysis of EEG and MEG event-related scalp field data using global randomization statistics

We present a program (Ragu; Randomization Graphical User interface) for statistical analyses of multichannel event-related EEG and MEG experiments. Based on measures of scalp field differences including all sensors, and using powerful, assumption-free randomization statistics, the program yields robust, physiologically meaningful conclusions based on the entire, untransformed, and unbiased set of measurements. Ragu accommodates up to two within-subject factors and one between-subject factor with multiple levels each. Significance is computed as function of time and can be controlled for type II errors with overall analyses. Results are displayed in an intuitive visual interface that allows further exploration of the findings. A sample analysis of an ERP experiment illustrates the different possibilities offered by Ragu. The aim of Ragu is to maximize statistical power while minimizing the need for a-priori choices of models and parameters (like inverse models or sensors of interest) that interact with and bias statistics.

Uso de DES y RANS no estacionario para la simulación del flujo alrededor de una carena Wigley - Use of Unsteady RANS and DES to simulate flows around a Wigley Hull)

Computer Fluid Dynamics tools have already become a valuable instrument for Naval Architects during the ship design process, thanks to their accuracy and the available computer power. Unfortunately, the development of RANSE codes, generally used when viscous effects play a major role in the flow, has not reached a mature stage, being the accuracy of the turbulence models and the free surface representation the most important sources of uncertainty. Another level of uncertainty is added when the simulations are carried out for unsteady flows, as those generally studied in seakeeping and maneuvering analysis and URANS equations solvers are used. Present work shows the applicability and the benefits derived from the use of new approaches for the turbulence modeling (Detached Eddy Simulation) and the free surface representation (Level Set) on the URANS equations solver CFDSHIP-Iowa. Compared to URANS, DES is expected to predict much broader frequency contents and behave better in flows where boundary layer separation plays a major role. Level Set methods are able to capture very complex free surface geometries, including breaking and overturning waves. The performance of these improvements is tested in set of fairly complex flows, generated by a Wigley hull at pure drift motion, with drift angle ranging from 10 to 60 degrees and at several Froude numbers to study the impact of its variation. Quantitative verification and validation are performed with the obtained results to guarantee their accuracy. The results show the capability of the CFDSHIP-Iowa code to carry out time-accurate simulations of complex flows of extreme unsteady ship maneuvers. The Level Set method is able to capture very complex geometries of the free surface and the use of DES in unsteady simulations highly improves the results obtained. Vortical structures and instabilities as a function of the drift angle and Fr are qualitatively identified. Overall analysis of the flow pattern shows a strong correlation between the vortical structures and free surface wave pattern. Karman-like vortex shedding is identified and the scaled St agrees well with the universal St value. Tip vortices are identified and the associated helical instabilities are analyzed. St using the hull length decreases with the increase of the distance along the vortex core (x), which is similar to results from other simulations. However, St scaled using distance along the vortex cores shows strong oscillations compared to almost constants for those previous simulations. The difference may be caused by the effect of the free-surface, grid resolution, and interaction between the tip vortex and other vortical structures, which needs further investigations. This study is exploratory in the sense that finer grids are desirable and experimental data is lacking for large α, especially for the local flow. More recently, high performance computational capability of CFDSHIP-Iowa V4 has been improved such that large scale computations are possible. DES for DTMB 5415 with bilge keels at α = 20º were conducted using three grids with 10M, 48M and 250M points. DES analysis for flows around KVLCC2 at α = 30º is analyzed using a 13M grid and compared with the results of DES on the 1.6M grid by. Both studies are consistent with what was concluded on grid resolution herein since dominant frequencies for shear-layer, Karman-like, horse-shoe and helical instabilities only show marginal variation on grid refinement. The penalties of using coarse grids are smaller frequency amplitude and less resolved TKE. Therefore finer grids should be used to improve V&V for resolving most of the active turbulent scales for all different Fr and α, which hopefully can be compared with additional EFD data for large α when it becomes available.