An improved finite element space for discontinuous pressures

We consider incompressible Stokes flow with an internal interface at which the pressure is discontinuous, as happens for example in problems involving surface tension. We assume that the mesh does not follow the interface, which makes classical interpolation spaces to yield suboptimal convergence rates (typically, the interpolation error in the L(2)(Omega)-norm is of order h(1/2)). We propose a modification of the P(1)-conforming space that accommodates discontinuities at the interface without introducing additional degrees of freedom or modifying the sparsity pattern of the linear system. The unknowns are the pressure values at the vertices of the mesh and the basis functions are computed locally at each element, so that the implementation of the proposed space into existing codes is straightforward. With this modification, numerical tests show that the interpolation order improves to O(h(3/2)). The new pressure space is implemented for the stable P(1)(+)/P(1) mini-element discretization, and for the stabilized equal-order P(1)/P(1) discretization. Assessment is carried out for Poiseuille flow with a forcing surface and for a static bubble. In all cases the proposed pressure space leads to improved convergence orders and to more accurate results than the standard P(1) space. In addition, two Navier-Stokes simulations with moving interfaces (Rayleigh-Taylor instability and merging bubbles) are reported to show that the proposed space is robust enough to carry out realistic simulations. (c) 2009 Elsevier B.V. All rights reserved.

Ab initio calculation of the BCC Fe-Al-Mo (Iron-Aluminum-Molybdenum) phase diagram: Implications for the nature of the tau(2) phase

The metastable phase diagram of the BCC-based ordering equilibria in the Fe-Al-Mo system has been calculated via a truncated cluster expansion, through the combination of Full-Potential-Linear augmented Plane Wave (FP-LAPW) electronic structure calculations and of Cluster Variation Method (CVM) thermodynamic calculations in the irregular tetrahedron approximation. Four isothermal sections at 1750 K, 2000 K, 2250 K and 2500 K are calculated and correlated with recently published experimental data on the system. The results confirm that the critical temperature for the order-disorder equilibrium between Fe(3)Al-D0(3) and FeAl-B2 is increased by Mo additions, while the critical temperature for the FeAl-B2/A2 equilibrium is kept approximately invariant with increasing Mo contents. The stabilization of the Al-rich A2 phase in equilibrium with overstoichiometric B2-(Fe,Mo)Al is also consistent with the attribution of the A2 structure to the tau(2) phase, stable at high temperatures in overstoichiometric B2-FeAl. (C) 2009 Elsevier Ltd. All rights reserved.

Bootstrap percolation on homogeneous trees has 2 phase transitions

We study the threshold theta bootstrap percolation model on the homogeneous tree with degree b + 1, 2 <= theta <= b, and initial density p. It is known that there exists a nontrivial critical value for p, which we call p(f), such that a) for p > p(f), the final bootstrapped configuration is fully occupied for almost every initial configuration, and b) if p < p(f) , then for almost every initial configuration, the final bootstrapped configuration has density of occupied vertices less than 1. In this paper, we establish the existence of a distinct critical value for p, p(c), such that 0 < p(c) < p(f), with the following properties: 1) if p <= p(c), then for almost every initial configuration there is no infinite cluster of occupied vertices in the final bootstrapped configuration; 2) if p > p(c), then for almost every initial configuration there are infinite clusters of occupied vertices in the final bootstrapped configuration. Moreover, we show that 3) for p < p(c), the distribution of the occupied cluster size in the final bootstrapped configuration has an exponential tail; 4) at p = p(c), the expected occupied cluster size in the final bootstrapped configuration is infinite; 5) the probability of percolation of occupied vertices in the final bootstrapped configuration is continuous on [0, p(f)] and analytic on (p(c), p(f) ), admitting an analytic continuation from the right at p (c) and, only in the case theta = b, also from the left at p(f).

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

Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures

We design and investigate a sequential discontinuous Galerkin method to approximate two-phase immiscible incompressible flows in heterogeneous porous media with discontinuous capillary pressures. The nonlinear interface conditions are enforced weakly through an adequate design of the penalties on interelement jumps of the pressure and the saturation. An accurate reconstruction of the total velocity is considered in the Raviart-Thomas(-Nedelec) finite element spaces, together with diffusivity-dependent weighted averages to cope with degeneracies in the saturation equation and with media heterogeneities. The proposed method is assessed on one-dimensional test cases exhibiting rough solutions, degeneracies, and capillary barriers. Stable and accurate solutions are obtained without limiters. (C) 2010 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.

Phase behavior of the orange essential oil/sodium bis(2-ethylhexyl)sulfosuccinate/water system

The ternary phase diagram for the orange essential oil (OEO)/sodium bis(2-ethylhexyl)sulfosuccinate (AOT)/water system was constructed at 25 degrees C. It indicates a large single phase region, comprising an isotropic water-in-oil (W/O) microemulsion (ME) phase (L(2)), a liquid crystal (LC) (lamellar or hexagonal) and a large unstable emulsion phase that separates in two phases of normal and reverse micelles (L(1) and L(2)). In this communication the properties of the ME are investigated by viscosity, electric conductivity and small angle X-ray scattering (SAXS) indicating that the isotropic ME phase exhibits different behaviors depending on composition. At low water content low viscous ""dry"" surfactant structures are formed, whereas at higher water content higher viscous water droplets are formed. The experimental data allow the determination of the transition from ""dry"" to the water droplet structures within the L(2) phase. SAXS analyses have also been performed for selected LC samples. (C) 2009 Elsevier B.V. All rights reserved.

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.

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids. (C) 2010 Elsevier B.V. All rights reserved.

An implicit technique for solving 3D low Reynolds number moving free surface flows

This paper describes the development of an implicit finite difference method for solving transient three-dimensional incompressible free surface flows. To reduce the CPU time of explicit low-Reynolds number calculations, we have combined a projection method with an implicit technique for treating the pressure on the free surface. The projection method is employed to uncouple the velocity and the pressure fields, allowing each variable to be solved separately. We employ the normal stress condition on the free surface to derive an implicit technique for calculating the pressure at the free surface. Numerical results demonstrate that this modification is essential for the construction of methods that are more stable than those provided by discretizing the free surface explicitly. In addition, we show that the proposed method can be applied to viscoelastic fluids. Numerical results include the simulation of jet buckling and extrudate swell for Reynolds numbers in the range [0.01, 0.5]. (C) 2008 Elsevier Inc. All rights reserved.

NUCLEAR SPIN 3/2 ELECTRIC QUADRUPOLE RELAXATION AS A QUANTUM COMPUTATION PROCESS

In this work we applied a quantum circuit treatment to describe the nuclear spin relaxation. From the Redfield theory, we obtain a description of the quadrupolar relaxation as a computational process in a spin 3/2 system, through a model in which the environment is comprised by five qubits and three different quantum noise channels. The interaction between the environment and the spin 3/2 nuclei is described by a quantum circuit fully compatible with the Redfield theory of relaxation. Theoretical predictions are compared to experimental data, a short review of quantum channels and relaxation in NMR qubits is also present.

Non-stationary frequency analysis of extreme daily rainfall in Sao Paulo, Brazil

This work is an assessment of frequency of extreme values (EVs) of daily rainfall in the city of Sao Paulo. Brazil, over the period 1933-2005, based on the peaks-over-threshold (POT) and Generalized Pareto Distribution (GPD) approach. Usually. a GPD model is fitted to a sample of POT Values Selected With a constant threshold. However. in this work we use time-dependent thresholds, composed of relatively large p quantities (for example p of 0.97) of daily rainfall amounts computed from all available data. Samples of POT values were extracted with several Values of p. Four different GPD models (GPD-1, GPD-2, GPD-3. and GDP-4) were fitted to each one of these samples by the maximum likelihood (ML) method. The shape parameter was assumed constant for the four models, but time-varying covariates were incorporated into scale parameter of GPD-2. GPD-3, and GPD-4, describing annual cycle in GPD-2. linear trend in GPD-3, and both annual cycle and linear trend in GPD-4. The GPD-1 with constant scale and shape parameters is the simplest model. For identification of the best model among the four models WC used rescaled Akaike Information Criterion (AIC) with second-order bias correction. This criterion isolates GPD-3 as the best model, i.e. the one with positive linear trend in the scale parameter. The slope of this trend is significant compared to the null hypothesis of no trend, for about 98% confidence level. The non-parametric Mann-Kendall test also showed presence of positive trend in the annual frequency of excess over high thresholds. with p-value being virtually zero. Therefore. there is strong evidence that high quantiles of daily rainfall in the city of Sao Paulo have been increasing in magnitude and frequency over time. For example. 0.99 quantiles of daily rainfall amount have increased by about 40 mm between 1933 and 2005. Copyright (C) 2008 Royal Meteorological Society

First-order-reversal-curve analysis of Pr-Fe-B-based exchange spring magnets

Ribbons of nominal composition (Pr(9.5)Fe(84.5)B(6))(0.96)Cr(0.01)(TiC)(0.03) were produced by arc-melting and melt-spinning the alloys on a Cu wheel. X-ray diffraction (XRD) reveals two main phases, one based upon alpha-Fe and the other upon Pr(2)Fe(14)B. The ribbons show exchange spring behavior with H (c) = 12.5 kOe and (BH)(max) = 13.6 MGOe when these two phases are well coupled. Transmission electron microscopy revealed the coupled behavior is observed when the microstructure consists predominantly of alpha-Fe grains (diameter similar to 100 nm.) surrounded by hard material containing Pr(2)Fe(14)B. The microstructure is discussed in terms of a calculation by Skomski and Coey. A first-order-reversal-curve (FORC) analysis was performed for both a well-coupled sample and a poorly coupled sample. The FORC diagrams show two strong peaks for both the poorly coupled sample and for the well-coupled material. In both cases, the localization of the FORC probability suggests magnetizing interactions between particles. Switching field distributions were calculated and are consistent with the sample microstructure.

First-order-reversal-curve analysis of Pr-Fe-B-based nanocomposites

Ribbons of nominal composition (Pr(9.5)Fe(84.5)B(6))(0.96)Cr(0.01)(TiC)(0.03) were produced by arc-melting and melt-spinning the alloys on a Cu wheel. X-ray diffraction reveals two main phases, one based upon alpha-Fe and the other upon Pr(2)Fe(14)B. The ribbons show exchange spring behavior with H(c)=12.5 kOe and (BH)(max)= 13.6 MGOe when these two phases are well coupled. Transmission electron microscopy revealed that the coupled behavior is observed when the microstructure consists predominantly of alpha-Fe grains(diameter similar to 100 nm.) surrounded by hard material containing Pr(2)Fe(14)B. A first-order-reversal-curve (FORC) analysis was performed for both a well-coupled sample and a partially-coupled sample. The FORC diagrams show two strong peaks for both the partially-coupled sample and for the well coupled material. In both cases, the localization of the FORC probability suggests demagnetizing interactions between particles. Switching field distributions were calculated and are consistent with the sample microstructure. (C) 2009 Elsevier B.V. All rights reserved.

Oxide surface modification: Synthesis and characterization of zirconia-coated alumina

Four aluminas were used as Supports for impregnation with a zirconium oxide with the aim to achieve a coating, without phase separation, between Support and modifier. The Supports were impregnated with different concentrations Of zirconium aqueous resin, obtained through the polymeric precursor method. After impregnation the samples were calcined and then characterized by XRD, which led to identification of crystalline zirconia in different concentrations from each support used. Using a simple geometric model the maximum amount Of Surface modifier Oxide required for the complete coating of a support with a layer of unit cells was estimated. According to this estimate, only the support should be identified below the limit proposed and crystalline zirconium oxide Should be identified above this limit when a complete coating is reached. The results obtained From XRD agree with the estimated values and to confirm the coating, the samples were also characterized by EDS/STEM, HRTEM, XPS, and XAS. The results showed that the zirconium oxide oil the Surface of alumina Support reached the coating in the limit of 15 Zr nm(-2), without the formation of the ZrO(2) phase. (c) 2009 Elsevier Inc. All rights reserved.