Application of a bounded upwinding scheme to complex fluid dynamics problems

This paper is concerned with an overview of upwinding schemes, and further nonlinear applications of a recently introduced high resolution upwind differencing scheme, namely the ADBQUICKEST [V.G. Ferreira, F.A. Kurokawa, R.A.B. Queiroz, M.K. Kaibara, C.M. Oishi, J.A.Cuminato, A.F. Castelo, M.F. Tomé, S. McKee, assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems, International Journal for Numerical Methods in Fluids 60 (2009) 1-26]. The ADBQUICKEST scheme is a new TVD version of the QUICKEST [B.P. Leonard, A stable and accurate convective modeling procedure based on quadratic upstream interpolation, Computer Methods in Applied Mechanics and Engineering 19 (1979) 59-98] for solving nonlinear balance laws. The scheme is based on the concept of NV and TVD formalisms and satisfies a convective boundedness criterion. The accuracy of the scheme is compared with other popularly used convective upwinding schemes (see, for example, Roe (1985) [19], Van Leer (1974) [18] and Arora & Roe (1997) [17]) for solving nonlinear conservation laws (for example, Buckley-Leverett, shallow water and Euler equations). The ADBQUICKEST scheme is then used to solve six types of fluid flow problems of increasing complexity: namely, 2D aerosol filtration by fibrous filters; axisymmetric flow in a tubular membrane; 2D two-phase flow in a fluidized bed; 2D compressible Orszag-Tang MHD vortex; axisymmetric jet onto a flat surface at low Reynolds number and full 3D incompressible flows involving moving free surfaces. The numerical simulations indicate that this convective upwinding scheme is a good generic alternative for solving complex fluid dynamics problems. © 2012.

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.

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.

On a linear steady-state system of equations in microcontinuum fluid mechanics

Galerkin representations and integral representations are obtained for the linearized system of coupled differential equations governing steady incompressible flow of a micropolar fluid. The special case of 2-dimensional Stokes flows is then examined and further representation formulae as well as asymptotic expressions, are generated for both the microrotation and velocity vectors. With the aid of these formulae, the Stokes Paradox for micropolar fluids is established.

Phase relationships and magnetic properties of TbxDy1-xFe1.95 alloys

The effect of Tb/Dy ratio on the structural and magnetic properties of (Tb,Dy)Fe-2 class of alloys has been investigated using nine alloys of TbxDy1-xFe1.95 (x = 0-1) covering the entire range. Our results indicate that the three phases viz. (Tb,Dy)Fe-2 (major phase), (Tb,Dy)Fe-3 and(Tb,Dy)-solid solution (minor phases) coexist in all the alloys. The volume fraction of pro-peritectic (Tb,Dy)Fe-3 phase however, has a minimum at x = 0.4 and a maximum at x = 0.6 compositions. The volume fraction of this phase decreases upon heat treatment at 850 degrees C and 1000 degrees C. A Widmanstatten type precipitate of (Tb,Dy)Fe-3 was observed for Dy-rich compositions (0 <= x <= 0.5). The microstructural investigations indicate that the ternary phase equilibria of Tb-Dy-Fe are sensitive to Tb/Dy ratio including the expansion of (Tb,Dy)Fe-2 phase field which is in contrast to the pseudo-binary assumption that is followed in available literature to date. The lattice parameter, Curie temperature and coercivity are found to increase with Tb addition. Split of (440) peak of (Tb,Dy)Fe-2 observed in x >= 0.3 alloys indicate, a spin reorientation transition from 100] to 111] occurs with Tb addition. (C) 2012 Elsevier B. V. All rights reserved.

Competing structural phase transition scenarios in the giant tetragonality ferroelectric BiFeO3-PbTiO3: Isostructural vs multiphase transition

A systematic x-ray and neutron powder diffraction study of the giant tetragonality multiferroic (1-x) BiFeO3-(x) PbTiO3 have revealed that the compositions close to the morphotropic phase boundary present two different structural phase transition scenarios on cooling from the cubic phase: (i) cubic -> tetragonal (T-2) + tetragonal (T-1) -> tetragonal (T-1) and (ii) cubic -> tetragonal (T-2) + tetragonal (T-1) + rhombohedral (R3c) -> tetragonal (T-1) + rhombohedral (R3c). The comparatively larger tetragonality (c/a - 1) of the T-1 phase as compared to the coexisting isostructural T-2 phase is shown to be a result of significantly greater degree of overlap of the Pb/Bi-6s and Ti/Fe-3d with the O-2p orbitals as compared to that in the T-2 phase. The formation/suppression of the minor metastable rhombohedral phase seems to be governed by subtle play of local kinetic factors. In the scenario when the minor rhombohedral (R) phase is formed along with the tetragonal phases it is able to accommodate the large transformation stress in the system due to formation of the tetragonal phases, and prevent the solid from disintegration into powder after sintering. When the metastable rhombohedral phase is not formed, the large transformation strain ruptures the grain boundaries leading to fragmentation of the dense solid to powder. (C) 2013 American Institute of Physics. http://dx.doi.org/10.1063/1.4792215]

Low temperature sulfurization of electrodeposited Cu(In,Al)Se-2 thin films

Sulfurization of Cu(In,Al)Se-2 films is carried out in an indigenously made set up at moderately low temperature. The films are sulfurized for different time durations of 15, 30, 45 and 60 min at 150 degrees C. InSe and Cu2S phases occurred in the films during the initial stage of sulfurization along with Cu(In,Al)(Se,S)(2) phase. The compositional analysis shows that the sulfur incorporation is saturated after 30 min. Crystallinity increased with the increase in sulfurization time. The band gap of the Cu(In,Al)Se-2 film increased up to 1.35 eV with the addition of sulfur. Single phase Cu(In,Al)(Se,S)(2) with high crystallinity is obtained after 60 min of sulfurization. (C) 2014 Elsevier B.V. All rights reserved.

Computation of Viscous Incompressible Flow using Pressure Correction Method on Unstructured Chimera Grid

In this paper, a pressure correction algorithm for computing incompressible flows is modified and implemented on unstructured Chimera grid. Schwarz method is used to couple the solutions of different sub-domains. A new interpolation to ensure consistency between primary variables and auxiliary variables is proposed. Other important issues such as global mass conservation and order of accuracy in the interpolations are also discussed. Two numerical simulations are successfully performed. They include one steady case, the lid-driven cavity and one unsteady case, the flow around a circular cylinder. The results demonstrate a very good performance of the proposed scheme on unstructured Chimera grids. It prevents the decoupling of pressure field in the overlapping region and requires only little modification to the existing unstructured Navier–Stokes (NS) solver. The numerical experiments show the reliability and potential of this method in applying to practical problems.

Study of drag reduction by gas injection for power-law fluid flow in horizontal stratified and slug flow regimes

In this work, the drag reduction by gas injection for power-law fluid flow in stratified and slug flow regimes has been studied. Experimentswere conducted to measure the pressure gradient within air/CMC solutions in a horizontal Plexiglas pipe that had a diameter of 50mm and a length of 30 m. The drag reduction ratio in stratified flow regime was predicted using the two-fluid model. The results showed that the drag reduction should occur over the large range of the liquid holdup when the flow behaviour index remained at the low value. Furthermore, for turbulent gas-laminar liquid stratified flow, the drag reduction by gas injection for Newtonian fluid was more effective than that for shear-shinning fluid, when the dimensionless liquid height remained in the area of high value. The pressure gradient model for a gas/Newtonian liquid slug flow was extended to liquids possessing the Ostwald–de Waele power law model. The proposed model was validated against 340 experimental data point over a wide range of operating conditions, fluid characteristics and pipe diameters. The dimensionless pressure drop predicted was well inside the 20% deviation region for most of the experimental data. These results substantiated the general validity of the model presented for gas/non-Newtonian two-phase slug flows.

Wave pattern characteristics of a two-phase nozzle flow by shock propagation

In this paper, the wave pattern characteristics of shock-induced two-phase nozzle Hows with the quiescent or moving dusty gas ahead of the incident-shock front has been investigated by using high-resolution numerical method. As compared with the corresponding results in single-phase nozzle flows of the pure gas, obvious differences between these two kinds of flows can be obtained.

Numerical simulations of unsteady wet steam flows with non-equilibrium condensation in the nozzle and the steam turbine

Numerical techniques for non-equilibrium condensing flows are presented. Conservation equations for homogeneous gas-liquid two-phase compressible flows are solved by using a finite volume method based on an approximate Riemann solver. The phase change consists of the homogeneous nucleation and growth of existing droplets. Nucleation is computed with the classical Volmer-Frenkel model, corrected for the influence of the droplet temperature being higher than the steam temperature due to latent heat release. For droplet growth, two types of heat transfer model between droplets and the surrounding steam are used: a free molecular flow model and a semi-empirical two-layer model which is deemed to be valid over a wide range of Knudsen number. The computed pressure distribution and Sauter mean droplet diameters in a convergent-divergent (Laval) nozzle are compared with experimental data. Both droplet growth models capture qualitatively the pressure increases due to sudden heat release by the non-equilibrium condensation. However the agreement between computed and experimental pressure distributions is better for the two-layer model. The droplet diameter calculated by this model also agrees well with the experimental value, whereas that predicted by the free molecular model is too small. Condensing flows in a steam turbine cascade are calculated at different Mach numbers and inlet superheat conditions and are compared with experiments. Static pressure traverses downstream from the blade and pressure distributions on the blade surface agree well with experimental results in all cases. Once again, droplet diameters computed with the two-layer model give best agreement with the experiments. Droplet sizes are found to vary across the blade pitch due to the significant variation in expansion rate. Flow patterns including oblique shock waves and condensation-induced pressure increases are also presented and are similar to those shown in the experimental Schlieren photographs. Finally, calculations are presented for periodically unsteady condensing flows in a low expansion rate, convergent-divergent (Laval) nozzle. Depending on the inlet stagnation subcooling, two types of self-excited oscillations appear: a symmetric mode at lower inlet subcooling and an asymmetric mode at higher subcooling. Plots of oscillation frequency versus inlet sub-cooling exhibit a hysteresis loop, in accord with observations made by other researchers for moist air flow. Copyright © 2006 by ASME.

Phase Diagram of [Amim]Cl plus Salt Aqueous Biphasic Systems and Its Application for [Amim]Cl Recovery

In this study, binodal curves and tie line data of [Amim]Cl + salt (K3PO4, K2HPO4, K2CO3) + water aqueous biphasic systems (ABS) were measured and correlated satisfactorily with the Merchuk equation and Othmer-Tobias and Bancroft equations, respectively. [Amim]Cl could be recovered from aqueous solutions using the ABS, and the recovery efficiency could reach 96.80%. The recovery efficiency was influenced by the concentrations of the salts and their Homeister series: K3PO4 > K2HPO4 > K2CO3. Our method provides a new and effective route for the recovery of hydrophilic IL using [Amim]Cl + salt + water ABS from aqueous solutions.

Defect evolution and hydrodynamic effects in lamellar ordering process of two-dimensional quenched block copolymers

The effects of hydrodynamic interactions on the lamellar ordering process for two-dimensional quenched block copolymers in the presence of extended defects and the topological defect evolutions in lamellar ordering process are numerically investigated by means of a model based on lattice Boltzmann method and self-consistent field theory. By observing the evolution of the average size of domains, it is found that the domain growth is faster with stronger hydrodynamic effects. The morphological patterns formed also appear different. To study the defect evolution, a defect density is defined and is used to explore the defect evolutions in lamellar ordering process. Our simulation results show that the hydrodynamics effects can reduce the density of defects. With our model, the relations between the Flory-Huggins interaction parameter chi, the length of the polymer chains N, and the defect evolutions are studied.

Phase diagram data for several salt plus salt aqueous biphasic systems at 298.15 K

Phase diagrams corresponding to aqueous biphasic systems of salt (the organic ionic liquid of salts [C(4)mim]Cl, [C(6)mim]Cl, and [C(8)mim]Cl) + salt (K3PO4, K2CO3) + water were determined at 298.15 K. The binodal curve was fitted to the Merchuk equation. Tie lines assigned from mass phase ratios according to the lever arm rule were satisfactorily described using the Othmer-Tobias and Bancroft equations.