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.

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.

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.

The T-2 phase in the Nb-Si-B system studied by ab initio calculations and synchrotron X-ray diffraction

The solid solution based on Nb5Si3 (Cr5B3 structure type, D8(l), tl32, 14/mcm, No140, a=6.5767 angstrom, c=11.8967 angstrom) in the Nb-Si-B system was studied from the structural and thermodynamic point of view both experimentally and by ab initio calculations. Rietveld refinement of powder X-ray synchrotron data allowed to determine the boron to silicon substitution mechanism and the structural parameters. Ab initio calculations of different ordered compounds and selected disordered alloys allowed to obtain in addition to the enthalpy of formation of the solution, substitution mechanism and structural parameters which are in excellent agreement with the experimental data. The stability of the phase is discussed. (C) 2012 Elsevier Inc. All rights reserved.

Robust tori-like Lagrangian coherent structures

In general the term "Lagrangian coherent structure" (LCS) is used to make reference about structures whose properties are similar to a time-dependent analog of stable and unstable manifolds from a hyperbolic fixed point in Hamiltonian systems. Recently, the term LCS was used to describe a different type of structure, whose properties are similar to those of invariant tori in certain classes of two-dimensional incompressible flows. A new kind of LCS was obtained. It consists of barriers, called robust tori that block the trajectories in certain regions of the phase space. We used the Double-Gyre Flow system as the model. In this system, the robust tori play the role of a skeleton for the dynamics and block, horizontally, vortices that come from different parts of the phase space. (C) 2012 Elsevier B.V. All rights reserved.

URANS calculations for smooth circular cylinder flow in a wide range of Reynolds Numbers: solution verification and validation

The flow around circular smooth fixed cylinder in a large range of Reynolds numbers is considered in this paper. In order to investigate this canonical case, we perform CFD calculations and apply verification & validation (V&V) procedures to draw conclusions regarding numerical error and, afterwards, assess the modeling errors and capabilities of this (U)RANS method to solve the problem. Eight Reynolds numbers between Re = 10 and Re 5 x 10(5) will be presented with, at least, four geometrically similar grids and five discretization in time for each case (when unsteady), together with strict control of iterative and round-off errors, allowing a consistent verification analysis with uncertainty estimation. Two-dimensional RANS, steady or unsteady, laminar or turbulent calculations are performed. The original 1994 k - omega SST turbulence model by Menter is used to model turbulence. The validation procedure is performed by comparing the numerical results with an extensive set of experimental results compiled from the literature. [DOI: 10.1115/1.4007571]

Saturated flow boiling heat transfer and critical heat flux in small horizontal flattened tubes

This paper presents experimental results for flow boiling heat transfer coefficient and critical heat flux (CHF) in small flattened tubes. The tested flattened tubes have the same equivalent internal diameter of 2.2 mm, but different aspect height/width ratios (H/W) of 1/4, 1/2, 2 and 4. The experimental data were compared against results for circular tubes using R134a and R245fa as working fluids at a nominal saturation temperature of 31 degrees C. For mass velocities higher than 200 kg/m(2)s, the flattened and circular tubes presented similar heat transfer coefficients. Such a behavior is related to the fact that stratification effects are negligible under conditions of higher mass velocities. Heat transfer correlations from the literature, usually developed using only circular-channel experimental data, predicted the flattened tube results for mass velocities higher than 200 kg/m(2)s with mean absolute error lower than 20% using the equivalent diameter to account for the geometry effect. Similarly, the critical heat flux results were found to be independent of the tube aspect ratio when the same equivalent length was kept. Equivalent length is a new parameter which takes into account the channel heat transfer area. The CHF correlations for round tubes predicted the flattened tube data relatively well when using the equivalent diameter and length. Furthermore, a new proposed CHF correlation predicted the present flattened tube data with a mean absolute error of 5%. (C) 2012 Elsevier Ltd. All rights reserved.

Geometrical and kinematic properties of interfacial waves in stratified oil-water flow in inclined pipe

The stratified oil-water flow pattern is common in the petroleum industry, especially in offshore directional wells and pipelines. Previous studies have shown that the phenomenon of flow pattern transition in stratified flow can be related to the interfacial wave structure (problem of hydrodynamic instability). The study of the wavy stratified flow pattern requires the characterization of the interfacial wave properties, i.e., average shape, celerity and geometric properties (amplitude and wavelength) as a function of holdup, inclination angle and phases' relative velocity. However, the data available in the literature on wavy stratified flow is scanty, especially in inclined pipes and when oil is viscous. This paper presents new geometric and kinematic interfacial wave properties as a function of a proposed two-phase Froude number in the wavy-stratified liquid-liquid flow. The experimental work was conducted in a glass test line of 12 m and 0.026 m id., oil (density and viscosity of 828 kg/m(3) and 0.3 Pa s at 20 degrees C, respectively) and water as the working fluids at several inclinations from horizontal (-20 degrees, -10 degrees, 0 degrees, 10 degrees, 20 degrees). The results suggest a physical relation between wave shape and the hydrodynamic stability of the stratified liquid-liquid flow pattern. (C) 2011 Elsevier Inc. All rights reserved.

Finite element methods for the Stokes system based on a Zienkiewicz type N-simplex

Hermite interpolation is increasingly showing to be a powerful numerical solution tool, as applied to different kinds of second order boundary value problems. In this work we present two Hermite finite element methods to solve viscous incompressible flows problems, in both two- and three-dimension space. In the two-dimensional case we use the Zienkiewicz triangle to represent the velocity field, and in the three-dimensional case an extension of this element to tetrahedra, still called a Zienkiewicz element. Taking as a model the Stokes system, the pressure is approximated with continuous functions, either piecewise linear or piecewise quadratic, according to the version of the Zienkiewicz element in use, that is, with either incomplete or complete cubics. The methods employ both the standard Galerkin or the Petrov–Galerkin formulation first proposed in Hughes et al. (1986) [18], based on the addition of a balance of force term. A priori error analyses point to optimal convergence rates for the PG approach, and for the Galerkin formulation too, at least in some particular cases. From the point of view of both accuracy and the global number of degrees of freedom, the new methods are shown to have a favorable cost-benefit ratio, as compared to velocity Lagrange finite elements of the same order, especially if the Galerkin approach is employed.

Gas-phase reactivity of protonated 2-oxazoline derivatives: mass spectrometry and computational studies

RATIONALE: Oxazolines have attracted the attention of researchers worldwide due to their versatility as carboxylic acid protecting groups, chiral auxiliaries, and ligands for asymmetric catalysis. Electrospray ionization tandem mass spectrometric (ESI-MS/MS) analysis of five 2-oxazoline derivatives has been conducted, in order to understand the influence of the side chain on the gas-phase dissociation of these protonated compounds under collision-induced dissociation (CID) conditions. METHODS: Mass spectrometric analyses were conducted in a quadrupole time-of-flight (Q-TOF) spectrometer fitted with electrospray ionization source. Protonation sites have been proposed on the basis of the gas-phase basicity, proton affinity, atomic charges, and a molecular electrostatic potential map obtained on the basis of the quantum chemistry calculations at the B3LYP/6-31 + G(d, p) and G2(MP2) levels. RESULTS: Analysis of the atomic charges, gas-phase basicity and proton affinities values indicates that the nitrogen atom is a possible proton acceptor site. On the basis of these results, two main fragmentation processes have been suggested: one taking place via neutral elimination of the oxazoline moiety (99 u) and another occurring by sequential elimination of neutral fragments with 72 u and 27 u. These processes should lead to formation of R+. CONCLUSIONS: The ESI-MS/MS experiments have shown that the side chain could affect the dissociation mechanism of protonated 2-oxazoline derivatives. For the compound that exhibits a hydroxyl at the lateral chain, water loss has been suggested to happen through an E2-type elimination, in an exothermic step. Copyright (C) 2012 John Wiley & Sons, Ltd.

Evaluation of the Effects of a Preoperative 2-Hour Fast With Maltodextrine and Glutamine on Insulin Resistance, Acute-Phase Response, Nitrogen Balance, and Serum Glutathione After Laparoscopic Cholecystectomy: A Controlled Randomized Trial

Background: Prolonged preoperative fasting increases insulin resistance (IR). The authors investigated whether an abbreviated preoperative fast with glutamine (GLN) plus a carbohydrate (CHO)-based beverage would improve the organic response after surgery. Methods: Forty-eight female patients (19-62 years) were randomized to either standard fasting (control group) or to fasting with 1 of 3 different beverages before video-cholecystectomy. Beverages were consumed 8 hours (400 mL; placebo group: water; GLN group: water with 50 g maltodextrine plus 40 g GLN; and CHO group: water with 50 g maltodextrine) and 2 hours (200 mL; placebo: water; GLN: water with 25 g maltodextrine plus 10 g GLN; and CHO: water with 25 g maltodextrine) before anesthesia. Blood samples were collected pre- and postoperatively. Results: The mean (SEM) postoperative homeostasis model assessment-insulin resistance was greater (P < .05) in control patients (4.3 [1.3]) than in the other groups (placebo, 1.6 [0.3]; CHO, 2.3 [0.4]; and GLN, 1.5 [0.1]). Glutathione was significantly higher (P < .01) in the GLN group than in both CHO and control groups. Interleukin-6 increased in all groups except the GLN group. The C-reactive protein/albumin ratio was higher (P < .05) in controls than in CHO and GLN groups. The nitrogen balance was less negative in GLN (-2.5 [0.8] gN) than in both placebo (-9.0 [2] gN; P = .001) and control (-6.6 [0.4] gN; P = .04) groups. Conclusions Preoperative intake of a GLN-enriched CHO beverage appears to improve IR and antioxidant defenses and decreases the inflammatory response after video-cholecystectomy. (JPEN J Parenter Enteral Nutr. 2012; 36: 43-52)

Fractionation of lemon essential oil by solvent extraction: Phase equilibrium for model systems at T=298.2 K

The fractioning of lemon essential oil can be performed by liquid-liquid extraction using hydrous ethanol as a solvent. A quaternary mixture composed of limonene, gamma-terpinene, beta-pinene, and citral was used to simulate lemon essential oil. In this paper, we present (liquid + liquid) equilibrium data that were experimentally determined for systems containing essential oil compounds, ethanol, and water at T = 298.2 K. The experimental data were correlated using the NRTL and UNIQUAC models, and the mean deviations between calculated and experimental data were less than 0.0053 in all systems, indicating the accuracy of these molecular models in describing our systems. The results show that as the water content in the solvent phase increased, the values of the distribution coefficients decreased, regardless of the type of compound studied. However, the oxygenated compound always showed the highest distribution coefficient among the components of the essential oil, thus making deterpenation of the lemon essential oil a feasible process. (C) 2012 Elsevier Ltd. All rights reserved.

Gas-phase reactivity of 2-hydroxy-1,4-naphthoquinones: a computational and mass spectrometry study of lapachol congeners

In order to understand the influence of alkyl side chains on the gas-phase reactivity of 1,4-naphthoquinone derivatives, some 2-hydroxy-1,4-naphthoquinone derivatives have been prepared and studied by electrospray ionization tandem mass spectrometry in combination with computational quantum chemistry calculations. Protonation and deprotonation sites were suggested on the basis of gas-phase basicity, proton affinity, gas-phase acidity (?Gacid), atomic charges and frontier orbital analyses. The nature of the intramolecular interaction as well as of the hydrogen bond in the systems was investigated by the atoms-in-molecules theory and the natural bond orbital analysis. The results were compared with data published for lapachol (2-hydroxy-3-(3-methyl-2-butenyl)-1,4-naphthoquinone). For the protonated molecules, water elimination was verified to occur at lower proportion when compared with side chain elimination, as evidenced in earlier studies on lapachol. The side chain at position C(3) was found to play important roles in the fragmentation mechanisms of these compounds. Copyright (c) 2012 John Wiley & Sons, Ltd.

Electron collisions with the HCOOH...(H2O)n (n=1, 2 and 3) complexes in liquid phase: the influence on the π* shape resonance of formic acid

In this conference we report cross sections for elastic collisions of low-energy electrons with the HCOOH…(H2O)n complexes, with n = 1, 2 and 3. The scattering cross sections were computed with the Schwinger multichannel method [K. Takatsuka and V. McKoy, Phys. Rev. A 24 , 2473 (1981); Phys. Rev. A 30 , 1734 (1984)] with pseudopotentials [M. H. F. Bettega, L. G. Ferreira, and M. A. P. Lima, Phys. Rev. A 47, 1111 (1993)] in the static-exchange and static-exchange plus polarization approximations, for energies from 0.5 eV to 6 eV. We considered some diÆerent hydrogen-bonded structures for the complexes that were generated with classical Monte Carlo simulations [K. Coutinho and S. Canuto, J. Chem. Phys. 113, 9132, (2000)]. The aim of this work is to investigate the effect of the surrounding water molecules on the π* shape resonance of the solute. Previous theoretical and experimental studies carried out in the gas phase reported a π* state for HCOOH at around 1.9 eV. For the n = 1 case and for all complexes, the stabilization of the resonance was observed (it appears at lower energy compared to the value obtained in the gas phase), as reported previously for the CH2O…H2O complexes [T. C. Freitas, M. A. P. Lima, S. Canuto, and M. H. F. Bettega, Phys. Rev. A 80, 062710 (2009)]. This result indicates that the presence of the solvent may affect the processes related to the π* state, such as the molecular dissociation by electron impact. For the n = 2 case we have observed both stabilization and destabilization of the π* resonance, that is associated with the hydrogen bond donor or acceptor role of the water molecules in the complexes. For the n = 3 case, preliminary static-exchange results show the stabilization of the π* state. We propose an explanation of the stabilization/destabilization of the π* state in terms of the polarization of the solute due to the surrounding water molecules and the net charge in the solute.

Two-Phase Frictional Pressure Drop and Flow Boiling Heat Transfer for R245fa in a 2.32-mm Tube

Experimental two-phase frictional pressure drop and flow boiling heat transfer results are presented for a horizontal 2.32-mm ID stainless-steel tube using R245fa as working fluid. The frictional pressure drop data was obtained under adiabatic and diabatic conditions. Experiments were performed for mass velocities ranging from 100 to 700 kg m−2 s−1 , heat flux from 0 to 55 kW m−2 , exit saturation temperatures of 31 and 41◦C, and vapor qualities from 0.10 to 0.99. Pressures drop gradients and heat transfer coefficients ranging from 1 to 70 kPa m−1 and from 1 to 7 kW m−2 K−1 were measured. It was found that the heat transfer coefficient is a strong function of the heat flux, mass velocity, and vapor quality. Five frictional pressure drop predictive methods were compared against the experimental database. The Cioncolini et al. (2009) method was found to work the best. Six flow boiling heat transfer predictive methods were also compared against the present database. Liu and Winterton (1991), Zhang et al. (2004), and Saitoh et al. (2007) were ranked as the best methods. They predicted the experimental flow boiling heat transfer data with an average error around 19%.