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.

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.

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

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

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.

CFD modeling of two-phase boiling flows in the slug flow regime with an interface capturing technique

The objective of this thesis was to improve the commercial CFD software Ansys Fluent to obtain a tool able to perform accurate simulations of flow boiling in the slug flow regime. The achievement of a reliable numerical framework allows a better understanding of the bubble and flow dynamics induced by the evaporation and makes possible the prediction of the wall heat transfer trends. In order to save computational time, the flow is modeled with an axisymmetrical formulation. Vapor and liquid phases are treated as incompressible and in laminar flow. By means of a single fluid approach, the flow equations are written as for a single phase flow, but discontinuities at the interface and interfacial effects need to be accounted for and discretized properly. Ansys Fluent provides a Volume Of Fluid technique to advect the interface and to map the discontinuous fluid properties throughout the flow domain. The interfacial effects are dominant in the boiling slug flow and the accuracy of their estimation is fundamental for the reliability of the solver. Self-implemented functions, developed ad-hoc, are introduced within the numerical code to compute the surface tension force and the rates of mass and energy exchange at the interface related to the evaporation. Several validation benchmarks assess the better performances of the improved software. Various adiabatic configurations are simulated in order to test the capability of the numerical framework in modeling actual flows and the comparison with experimental results is very positive. The simulation of a single evaporating bubble underlines the dominant effect on the global heat transfer rate of the local transient heat convection in the liquid after the bubble transit. The simulation of multiple evaporating bubbles flowing in sequence shows that their mutual influence can strongly enhance the heat transfer coefficient, up to twice the single phase flow value.

Occlusal outcomes and efficiency of 1-and 2-phase protocols in the treatment of Class II Division 1 malocclusion

Introduction: The purpose of this study was to compare the occlusal outcomes and the efficiency of 1-phase and 2-phase treatment protocols in Class II Division 1 malocclusions. Treatment efficiency was defined as a change in the occlusal characteristics in a shorter treatment time. Methods: Class II Division 1 subjects ( n = 139) were divided into 2 groups according to the treatment protocol for Class II correction. Group 1 comprised 78 patients treated with a 1-phase treatment protocol at initial and final mean ages of 12.51 and 14.68 years. Group 2 comprised 61 patients treated with a 2-phase treatment protocol at initial and final mean ages of 11.21 and 14.70 years. Lateral cephalometric radiographs were taken at the pretreatment stage to evaluate morphological differences in the groups. The initial and final study models of the patients were evaluated by using the peer assessment rating index. Chi-square tests were used to test for differences between the 2 groups for categorical variables. Variables regarding occlusal results were compared by using independent t tests. A linear regression analysis was completed, with total treatment time as the dependent variable, to identify clinical factors that predict treatment length for patients with Class II malocclusions. Results: Similar occlusal outcomes were obtained between the 1-phase and the 2-phase treatment protocols, but the duration of treatment was significantly shorter in the 1-phase treatment protocol group. Conclusions: Treatment of Class II Division 1 malocclusions is more efficient with the 1-phase than the 2-phase treatment protocol.

Investigation into Improved Pavement Curing Materials and Techniques: Part 2 - Phase III, March 2003

Appropriate curing is important for concrete to obtain the designed properties. This research was conducted to evaluate the curing effects of different curing materials and methods on pavement properties. At present the sprayed curing compound is a common used method for pavement and other concrete structure construction. Three curing compounds were selected for testing. Two different application rates were employed for the white-pigmented liquid curing compounds. The concrete properties of temperature, moisture content, conductivity, and permeability were examined at several test locations. It was found, in this project, that the concrete properties varied with the depth. Of the tests conducted (maturity, sorptivity, permeability, and conductivity), conductivity appears to be the best method to evaluate the curing effects in the field and bears potential for field application. The results indicated that currently approved curing materials in Iowa, when spread uniformly in a single or double application, provide adequate curing protection and meet the goals of the Iowa Department of Transportation. Experimental curing methods can be compared to this method through the use of conductivity testing to determine their application in the field.

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).

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.

A fast, high resolution, second-order central scheme for incompressible flows

A high resolution, second-order central difference method for incompressible flows is presented. The method is based on a recent second-order extension of the classic Lax–Friedrichs scheme introduced for hyperbolic conservation laws (Nessyahu H. & Tadmor E. (1990) J. Comp. Physics. 87, 408-463; Jiang G.-S. & Tadmor E. (1996) UCLA CAM Report 96-36, SIAM J. Sci. Comput., in press) and augmented by a new discrete Hodge projection. The projection is exact, yet the discrete Laplacian operator retains a compact stencil. The scheme is fast, easy to implement, and readily generalizable. Its performance was tested on the standard periodic double shear-layer problem; no spurious vorticity patterns appear when the flow is underresolved. A short discussion of numerical boundary conditions is also given, along with a numerical example.

Inhibition of S/G(2) phase CDK4 reduces mitotic fidelity

Cyclin-dependent kinase 4 (CDK4)/cyclin D has a key role in regulating progression through late G(1) into S phase of the cell cycle. CDK4-cyclin D complexes then persist through the latter phases of the cell cycle, although little is known about their potential roles. We have developed small molecule inhibitors that are highly selective for CDK4 and have used these to define a role for CDK4-cyclin D in G(2) phase. The addition of the CDK4 inhibitor or small interfering RNA knockdown of cyclin D3, the cyclin D partner, delayed progression through G(2) phase and mitosis. The G(2) phase delay was independent of ATM/ATR and p38 MAPK but associated with elevated Wee1. The mitotic delay was because of failure of chromosomes to migrate to the metaphase plate. However, cells eventually exited mitosis, with a resultant increase in cells with multiple or micronuclei. Inhibiting CDK4 delayed the expression of the chromosomal passenger proteins survivin and borealin, although this was unlikely to account for the mitotic phenotype. These data provide evidence for a novel function for CDK4-cyclin D3 activity in S and G(2) phase that is critical for G(2)/M progression and the fidelity of mitosis.