A marker-and-cell approach to free surface 2-D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front-tracking method. The velocity field is computed using a finite-difference discretization of a modification of the NavierStokes equations. These equations together with the continuity equation are solved for the two-dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright (c) 2012 John Wiley & Sons, Ltd.

Interpolation estimate for a finite-element space with embedded discontinuities

We consider a recently proposed finite-element space that consists of piecewise affine functions with discontinuities across a smooth given interface Γ (a curve in two dimensions, a surface in three dimensions). Contrary to existing extended finite element methodologies, the space is a variant of the standard conforming Formula space that can be implemented element by element. Further, it neither introduces new unknowns nor deteriorates the sparsity structure. It is proved that, for u arbitrary in Formula, the interpolant Formula defined by this new space satisfies Graphic where h is the mesh size, Formula is the domain, Formula, Formula, Formula and standard notation has been adopted for the function spaces. This result proves the good approximation properties of the finite-element space as compared to any space consisting of functions that are continuous across Γ, which would yield an error in the Formula-norm of order Graphic. These properties make this space especially attractive for approximating the pressure in problems with surface tension or other immersed interfaces that lead to discontinuities in the pressure field. Furthermore, the result still holds for interfaces that end within the domain, as happens for example in cracked domains.

Numerical assessment of stability of interface discontinuous finite element pressure spaces

The stability of two recently developed pressure spaces has been assessed numerically: The space proposed by Ausas et al. [R.F. Ausas, F.S. Sousa, G.C. Buscaglia, An improved finite element space for discontinuous pressures, Comput. Methods Appl. Mech. Engrg. 199 (2010) 1019-1031], which is capable of representing discontinuous pressures, and the space proposed by Coppola-Owen and Codina [A.H. Coppola-Owen, R. Codina, Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions, Int. J. Numer. Methods Fluids, 49 (2005) 1287-1304], which can represent discontinuities in pressure gradients. We assess the stability of these spaces by numerically computing the inf-sup constants of several meshes. The inf-sup constant results as the solution of a generalized eigenvalue problems. Both spaces are in this way confirmed to be stable in their original form. An application of the same numerical assessment tool to the stabilized equal-order P-1/P-1 formulation is then reported. An interesting finding is that the stabilization coefficient can be safely set to zero in an arbitrary band of elements without compromising the formulation's stability. An analogous result is also reported for the mini-element P-1(+)/P-1 when the velocity bubbles are removed in an arbitrary band of elements. (C) 2012 Elsevier B.V. All rights reserved.

Transiting exoplanets from the CoRoT space mission XX. CoRoT-20b: A very high density, high eccentricity transiting giant planet

We report the discovery by the CoRoT space mission of a new giant planet, CoRoT-20b. The planet has a mass of 4.24 +/- 0.23 M-Jup and a radius of 0.84 +/- 0.04 R-Jup. With a mean density of 8.87 +/- 1.10 g cm(-3), it is among the most compact planets known so far. Evolutionary models for the planet suggest a mass of heavy elements of the order of 800 M-circle plus if embedded in a central core, requiring a revision either of the planet formation models or both planet evolution and structure models. We note however that smaller amounts of heavy elements are expected by more realistic models in which they are mixed throughout the envelope. The planet orbits a G-type star with an orbital period of 9.24 days and an eccentricity of 0.56. The star's projected rotational velocity is v sin i = 4.5 +/- 1.0 km s(-1), corresponding to a spin period of 11.5 +/- 3.1 days if its axis of rotation is perpendicular to the orbital plane. In the framework of Darwinian theories and neglecting stellar magnetic breaking, we calculate the tidal evolution of the system and show that CoRoT-20b is presently one of the very few Darwin-stable planets that is evolving toward a triple synchronous state with equality of the orbital, planetary and stellar spin periods.

A new enrichment space for the treatment of discontinuous pressures in multi-fluid flows

In this work, a new enrichment space to accommodate jumps in the pressure field at immersed interfaces in finite element formulations, is proposed. The new enrichment adds two degrees of freedom per element that can be eliminated by means of static condensation. The new space is tested and compared with the classical P1 space and to the space proposed by Ausas et al (Comp. Meth. Appl. Mech. Eng., Vol. 199, 10191031, 2010) in several problems involving jumps in the viscosity and/or the presence of singular forces at interfaces not conforming with the element edges. The combination of this enrichment space with another enrichment that accommodates discontinuities in the pressure gradient has also been explored, exhibiting excellent results in problems involving jumps in the density or the volume forces. Copyright (c) 2011 John Wiley & Sons, Ltd.

Hofmeister effects on the colloidal stability of poly(ethylene glycol)-decorated nanoparticles

The colloidal stability of poly(ethylene glycol)-decorated poly(methyl methacrylate), PMMA/Tween-20, particles was investigated by means of phase separation measurements, in the presence of sodium fluoride (NaF), sodium chloride, sodium bromide, sodium nitrate, or sodium thiocyanate (NaSCN) at 1.0 mol L-1. Following Hofmeister's series, the dispersions of PMMA/Tween-20 destabilized faster in the presence of NaF than with NaSCN. After the phase separation, the systems were homogenized and except for the dispersions in NaF, re-dispersed particles took longer to destabilize, indicating that anions adsorbed on the particles, creating a new surface. Except for F- ions, the adsorption of anions on the polar outmost shell was evidenced by means of tensiometry and small-angle X-ray scattering measurements. Fluoride ions induced the dehydration of the polar shell, without affecting the polar shell electron density, and the formation of very large aggregates. A model was proposed to explain the colloidal behavior in the presence of Hofmeister ions.

Comparison of physical properties of two exogenous surfactants: new parameter.

The pulmonary surfactant has essential physical properties for normal lung function. The most important property is the surface tension. In this work, it was evaluated the surface tension of two commercial exogenous surfactants used in surfactant replacement therapy, poractant alfa (Curosurf, Chiesi Farmaceuticals, Italy) and beractant (Survanta, Abbott Laboratories, USA) using new parameters. A Langmuir film balance (Minitrough, KSV Instruments, Finland) was used to measure surface tension of poractant alfa and beractant samples. For both samples, we prepared a solution of 1 mg/m dissolved in chloroform (100π`), which was applied over a subphase of milli-Q water (175 ml) in the chamber of the balance. The chamber has two moving barriers that can change its surface area between a maximal value of 112.5 cm 2 , and a minimal value of 22.5 cm 2, defining a balance cycle. Each surfactant had its surface tension evaluated during 20 balance cycles for three times. Four quantities were calculated from the experiment: Minimum Surface Tension (MTS), defined as the surface tension at minimal surface area during the first cycle; Mean Work Cycle (MWC), defined as the mean hysteresis area of the measured surface tension curve of the last 16 balance cycles; Critical Active Surface Area in Compression (CASAC) or in Expansion (CASAE), defined as the maximal chamber area where the surfactant is active on the surface in compression or expansion. The t-test was applied to verify for statistical significance of the results. Comproved with the MST is the same reported in literature, the differences between MWC, CASAC, and CASAE were statistically significant (p<0.001). The MWC, CASAC and CASAE were higher for poractant alfa than for beractant. A higher MWC for poractant alfa means higher elastic recoil of the lung in comparison with beractant. Using a different methodology, our results showed that poractant alfa is probably more effective in a surfactant replacement therapy than beractant due the use of poractant alfa in relation to the use of beractant in preterm infants with Respiratory Distress Syndrome (RDS).

Mechanical Measurements Of Two Exogenous Surfactants Using A Langmuir Film Balance:A Comparative Study

Rationale: The primary function of surfactant is to reduce the surface tension at air-liquid interface. In this study, the surface tension behavior of two commercial surfactants, poractant alfa (ChiesiFarmaceuticals,ltaly) and beractant (Abbott Laboratories,USA), were evaluated,using new parameters. Methods: We used a Langmuir film balance (Minitrough,KSV lnstruments,Finland) to measure of surface tension of both poractant alfa and beractant samples. For both samples,we prepared a solution of 1 mg/mdl dissolved in chloroform. The solution (1uL) was applied over a subphase of milli-Q water (175 ml) in the chamber of the balance. The chamber has two moving barriers that can change its surface area between a maximum value of 112.5 cm2 anda minimum value of 22.5 cm2, defining a balance cycle.lhree sample's films were evaluated for each sample, during 20 balance cycles. Here quantify two new variables, which is the mean hysteresis area of the measured curve surface tension of the last 16 balance cycles,defined here as Mean Work Cycle (MWC), and the moment that the surfactant is active in the surface, this measure is defined here as Active Surface Area Critical in the compression (ASACC) and the expansion (ASACE). The test was applied to compare the statistical significance of the results.

Is the meander growth in the Brazil Current system off Southeast Brazil due to baroclinic instability?

Temporally-growing frontal meandering and occasional eddy-shedding is observed in the Brazil Current (BC) as it flows adjacent to the Brazilian Coast. No study of the dynamics of this phenomenon has been conducted to date in the region between 22 degrees S and 25 degrees S. Within this latitude range, the flow over the intermediate continental slope is marked by a current inversion at a depth that is associated with the Intermediate Western Boundary Current (IWBC). A time series analysis of 10-current-meter mooring data was used to describe a mean vertical profile for the BC-IWBC jet and a typical meander vertical structure. The latter was obtained by an empirical orthogonal function (EOF) analysis that showed a single mode explaining 82% of the total variance. This mode structure decayed sharply with depth, revealing that the meandering is much more vigorous within the BC domain than it is in the IWBC region. As the spectral analysis of the mode amplitude time series revealed no significant periods, we searched for dominant wavelengths. This search was done via a spatial EOF analysis on 51 thermal front patterns derived from digitized AVHRR images. Four modes were statistically significant at the 95% confidence level. Modes 3 and 4, which together explained 18% of the total variance, are associated with 266 and 338-km vorticity waves, respectively. With this new information derived from the data, the [Johns, W.E., 1988. One-dimensional baroclinically unstable waves on the Gulf Stream potential vorticity gradient near Cape Hatteras. Dyn. Atmos. Oceans 11, 323-350] one-dimensional quasi-geostrophic model was applied to the interpolated mean BC-IWBC jet. The results indicated that the BC system is indeed baroclinically unstable and that the wavelengths depicted in the thermal front analysis are associated with the most unstable waves produced by the model. Growth rates were about 0.06 (0.05) days(-1) for the 266-km (338-km) wave. Moreover, phase speeds for these waves were low compared to the surface BC velocity and may account for remarks in the literature about growing standing or stationary meanders off southeast Brazil. The theoretical vertical structure modes associated with these waves resembled very closely to the one obtained for the current-meter mooring EOF analysis. We interpret this agreement as a confirmation that baroclinic instability is an important mechanism in meander growth in the BC system. (C) 2008 Elsevier B.V. All rights reserved.

Global interannual trends and amplitude modulations of the sea surface height anomaly from the TOPEX/Jason-1 altimeters

This study uses the global Ocean Topography Experiment (TOPEX)/Jason-1 altimeters` time series to estimate the 13-yr trend in sea surface height anomaly. These trends are estimated at each grid point by two methods: one fits a straight line to the time series and the other is based on the difference between the average height between the two halves of the time series. In both cases the trend shows large regional variability, mostly where the intense western boundary currents turn. The authors hypothesize that the regional variability of the sea surface height trends leads to changes in the local geostrophic transport. This in turn affects the instability-related processes that generate mesoscale eddies and enhances the Rossby wave signals. This hypothesis is verified by estimates of the trend of the amplitude of the filtered sea surface height anomaly that contains the spectral bands associated with Rossby waves and mesoscale eddies. The authors found predominantly positive tendency in the amplitude of Rossby waves and eddies, which suggests that, on average, these events are becoming more energetic. In some regions, the variation in amplitude over 13 yr is comparable to the standard deviation of the data and is statistically significant according to both methods employed in this study. It is plausible that in this case, the energy is transferred from the mean currents to the waves and eddies through barotropic and baroclinic instability processes that are more pronounced in the western boundary current extension regions. If these heat storage patterns and trends are confirmed on longer time series, then it will be justified to argue that the warming trend of the last century provides the energy that amplifies both Rossby waves and mesoscale eddies.

Synthetic Phytochelatin Surface Display in Cupriavidus metallidurans CH34 for Enhanced Metals Bioremediation

This work describes the effects of the cell surface display of a synthetic phytochelatin in the highly metal tolerant bacterium Cupriavidus metallidurans CH34. The EC20sp synthetic phytochelatin gene was fused between the coding sequences of the signal peptide (SS) and of the autotransporter beta-domain of the Neisseria gonorrhoeae IgA protease precursor (IgA beta), which successfully targeted the hybrid protein toward the C. metallidurans outer membrane. The expression of the SS-EC20sp-IgA beta gene fusion was driven by a modified version of the Bacillus subtilis mrgA promoter showing high level basal gene expression that is further enhanced by metal presence in C. metallidurans. The recombinant strain showed increased ability to immobilize Pb2+, Zn2+, Cu2+, Cd2+, Mn2+, and Ni2+ ions from the external medium when compared to the control strain. To ensure plasmid stability and biological containment, the MOB region of the plasmid was replaced by the E. coli hok/sok coding sequence.

Representation and analysis of control surface freeplay nonlinearity

Different representations for a control surface freeplay nonlinearity in a three degree of freedom aeroelastic system are assessed. These are the discontinuous, polynomial and hyperbolic tangent representations. The Duhamel formulation is used to model the aerodynamic loads. Assessment of the validity of these representations is performed through comparison with previous experimental observations. The results show that the instability and nonlinear response characteristics are accurately predicted when using the discontinuous and hyperbolic tangent representations. On the other hand, the polynomial representation fails to predict chaotic motions observed in the experiments. (c) 2012 Elsevier Ltd. All rights reserved.

Impact of the Manaus urban plume on trace gas mixing ratios near the surface in the Amazon Basin: Implications for the NO-NO2-O-3 photostationary state and peroxy radical levels

We measured the mixing ratios of NO, NO2, O-3, and volatile organic carbon as well as the aerosol light-scattering coefficient on a boat platform cruising on rivers downwind of the city of Manaus (Amazonas State, Brazil) in July 2001 (Large-Scale Biosphere-Atmosphere Experiment in Amazonia-Cooperative LBA Airborne Regional Experiment-2001). The dispersion and impact of the Manaus plume was investigated by a combined analysis of ground-based (boat platform) and airborne trace gas and aerosol measurements as well as by meteorological measurements complemented by dispersion calculations (Hybrid Single-Particle Lagrangian Integrated Trajectory model). For the cases with the least anthropogenic influence (including a location in a so far unexplored region similar to 150 km west of Manaus on the Rio Manacapuru), the aerosol scattering coefficient, sigma(s), was below 11 Mm(-1), NOx mixing ratios remained below 0.6 ppb, daytime O-3 mixing ratios were mostly below 20 ppb and maximal isoprene mixing ratios were about 3 ppb in the afternoon. The photostationary state (PSS) was not established for these cases, as indicated by values of the Leighton ratio, Phi, well above unity. Due to the influence of river breeze systems and other thermally driven mesoscale circulations, a change of the synoptic wind direction from east-northeast to south-southeast in the afternoon often caused a substantial increase of ss and trace gas mixing ratios (about threefold for sigma(s), fivefold for NOx, and twofold for O-3), which was associated with the arrival of the Manaus pollution plume at the boat location. The ratio F reached unity within its uncertainty range at NOx mixing ratios of about 3 ppb, indicating "steady-state" conditions in cases when radiation variations, dry deposition, emissions, and reactions mostly involving peroxy radicals (XO2) played a minor role. The median midday/afternoon XO2 mixing ratios estimated using the PSS method range from 90 to 120 parts per trillion (ppt) for the remote cases (sigma(s) < 11 Mm(-1) and NOx < 0.6 ppb), while for the polluted cases our estimates are 15 to 60 ppt. These values are within the range of XO2 estimated by an atmospheric chemistry box model (Chemistry As A Box model Application-Module Efficiently Calculating the Chemistry of the Atmosphere (CAABA/MECCA)-3.0).

Selectivity and Mechanisms Driven by Reaction Dynamics: The Case of the Gas-Phase OH-+CH3ONO2 Reaction

Well-established statistical approaches such as transition-state theory based on high-level calculated potential energy profiles are unable to account for the selectivity observed in the gas-phase OH- + CH3ONO2 reaction. This reaction can undergo bimolecular nucleophilic displacement at either the carbon center (S(N)2@C) or the nitrogen center (S(N)2@N) as well as a proton abstraction followed by dissociation (E(CO)2) pathway. Direct dynamics simulations yield an S(N)2:E(CO)2 product ratio in close agreement with experiment and show that the lack of reactivity at the nitrogen atom is due to the highly negative electrostatic potential generated by the oxygen atoms in the ONO2 group that scatters the incoming OH-. In addition to these dynamical effects, the nonstatistical behavior of these reactions is attributed to the absence of equilibrated reactant complexes and to the large number of recrossings, which might be present in several ion-molecule gas-phase reactions.

Surface modification by metal ion implantation forming metallic nanoparticles in insulating matrix.

There is special interest in the incorporation of metallic nanoparticles in a surrounding dielectric matrix for obtaining composites with desirable characteristics such as for surface plasmon resonance, which can be used in photonics and sensing, and controlled surface electrical conductivity. We investigated nanocomposites produced through metallic ion implantation in insulating substrate, where the implanted metal self-assembles into nanoparticles. During the implantation, the excess of metal atom concentration above the solubility limit leads to nucleation and growth of metal nanoparticles, driven by the temperature and temperature gradients within the implanted sample including the beam-induced thermal characteristics. The nanoparticles nucleate near the maximum of the implantation depth profile (projected range), that can be estimated by computer simulation using the TRIDYN. This is a Monte Carlo simulation program based on the TRIM (Transport and Range of Ions in Matter) code that takes into account compositional changes in the substrate due to two factors: previously implanted dopant atoms, and sputtering of the substrate surface. Our study suggests that the nanoparticles form a bidimentional array buried few nanometers below the substrate surface. More specifically we have studied Au/PMMA (polymethylmethacrylate), Pt/PMMA, Ti/alumina and Au/alumina systems. Transmission electron microscopy of the implanted samples showed the metallic nanoparticles formed in the insulating matrix. The nanocomposites were characterized by measuring the resistivity of the composite layer as function of the dose implanted. These experimental results were compared with a model based on percolation theory, in which electron transport through the composite is explained by conduction through a random resistor network formed by the metallic nanoparticles. Excellent agreement was found between the experimental results and the predictions of the theory. It was possible to conclude, in all cases, that the conductivity process is due only to percolation (when the conducting elements are in geometric contact) and that the contribution from tunneling conduction is negligible.