Two-Phase Inertial Flow in Homogeneous Porous Media: A Theoretical Derivation of a Macroscopic Model

The purpose of this article is to derive a macroscopic model for a certain class of inertial two-phase, incompressible, Newtonian fluid flow through homogenous porous media. Starting from the continuity and Navier–Stokes equations in each phase β and γ, the method of volume averaging is employed subjected to constraints that are explicitly provided to obtain the macroscopic mass and momentum balance equations. These constraints are on the length- and time-scales, as well as, on some quantities involving capillary, Weber and Reynolds numbers that define the class of two-phase flow under consideration. The resulting macroscopic momentum equation relates the phase-averaged pressure gradient to the filtration or Darcy velocity in a coupled nonlinear form explicitly given by

Transient shear and elongational behavior of blends of PET with a LCP

Blends of polyethylene terephthalate (PET) with a liquid crystalline polymer (LCP) and a compatibilizer were produced by twin screw extrusion and injection molding. Transesterification and compatibilization studies were made in a torque rheometer. The morphology of the injection-molded plaques was studied by scanning electron microscopy. The blends shear growth function was measured in a cone and plate rheometer. The elongational growth function was measured in a modified rotational rheometer. Transesterification was observed in the PET/LCP/compatibilizer 95/5/0 blend. The injection-molded plaques displayed the usual “skin-core” morphology. All the blends were highly shear-thinning, even at low shear rates; thus, a zero-shear viscosity could not be calculated. The compatibilized blend had the highest shear viscosity of all the blends, confirming the strong PET/LCP interphase and the effectiveness of the compatibilizing agent. On the other hand, the 90/10/0 blend had the lowest shear viscosity. All the blends showed strain softening behavior, similar to the PET. The 90/10/0 blend had the highest elongational growth function, while the 95/5/0 had the lowest. The compatibilized blend had an intermediate behavior between both blends.     

Evolution of viscosities and morphology for the two-phase system polyethylene oxide/poly(dimethylsiloxane)

Shear viscosities and oscillatory viscosities were measured for the two-phase system polyethylene oxide/poly(dimethylsiloxane) at 70 °C as a function of composition. This blend exhibits the usual droplet/matrix structures in the vicinity of the pure components and a region of co-continuity within which two droplet/matrix structures coexist. A stepwise reduction in the shear rate, , leads to a rapid increase in viscosity followed by a much slower exponential decay; plots of the corresponding rate constants as a function of composition exhibit two discontinuities marking the boundaries of co-continuity; a similar dependence is obtained for the time independent final viscosities . Keeping the blend composition constant and determining as a function of yields a curve that passes a distinct maximum, where the viscosities are very close to that of the less viscous pure component on both ends of this dependence. The dynamic mechanical measurements of the blends yield at low frequencies storage moduli G′ that are orders of magnitude larger than that of the components because of the deformation of the interfaces. At high frequencies, the loss moduli G″ reflect the increasing alignment of the drops suspended in the matrix phases. The composition dependencies of G′ and of the complex viscosities can again be used to determine the limits of co-continuity.     

Rotating Navier-Stokes Equations in $${\mathbb R}^{3}_{+}$$ with Initial Data Nondecreasing at Infinity: The Ekman Boundary Layer Problem

We prove time local existence and uniqueness of solutions to a boundary layer problem in a rotating frame around the stationary solution called the Ekman spiral. We choose initial data in the vector-valued homogeneous Besov space for 2 <  p <  ∞. Here the L ^p -integrability is imposed in the normal direction, while we may have no decay in tangential components, since the Besov space contains nondecaying functions such as almost periodic functions. A crucial ingredient is theory for vector-valued homogeneous Besov spaces. For instance we provide and apply an operator-valued bounded H ^∞-calculus for the Laplacian in for a general Banach space .     

Fluid flow and heat transfer in microchannels with rectangular cross section

Forced convective heat transfer coefficients and friction factors for flow of water in microchannels with a rectangular cross section were measured. An integrated microsystem consisting of five microchannels on one side and a localized heater and seven polysilicon temperature sensors along the selected channels on the other side was fabricated using a double-polished-prime silicon wafer. For the microchannels tested, the friction factor constant obtained are values between 53.7 and 60.4, which are close to the theoretical value from a correlation for macroscopic dimension, 56.9 for D _h  = 100 μm. The heat transfer coefficients obtained by measuring the wall temperature along the micro channels were linearly dependent on the wall temperature, in turn, the heat transfer mechanism is strongly dependent on the fluid properties such as viscosity. The measured Nusselt number in the laminar flow regime tested could be correlated by which is quite different from the constant value obtained in macrochannels.     

Effects of surfactant and salt concentrations on capillary flow and its entry flow for wormlike micelle solutions

Dynamic viscoelasticities and flow properties were measured for aqueous solutions of cetyltrimethylammonium bromide (CTAB) and sodium salicylate (NaSal) to examine the effects of surfactant (C _D) and salt (C _S). The relaxation time λ of a single mode Maxwell model was obtained, and the relationship between λ and free NaSal concentration was discussed. The relation between λ and was applied to the classification of flow curves, which were obtained using a capillary rheometer. In the flow curves, a shear rate jump occurred at low shear rates for the solutions with low , while bending was seen at high shear rates for all the flow curves. On the other hand, vortex growth at the salient corner in the entrance region of the capillary was also investigated. Four different flow patterns were identified: Newtonian-like flow (A), steady vortex flow (B), periodically oscillated flow (C), and perfectly unstable flow (D). In the steady vortex of the flow pattern B, the vortex length increased with increasing shear rate. In the flow patterns C and D, white turbidity was observed. Furthermore, the relation between λ and was also applied to the discussion on the development of the vortex.     

Resolvent Estimates and Maximal Regularity in Weighted L q -spaces of the Stokes Operator in an Infinite Cylinder

Let be an infinite cylinder of , n ≥ 3, with a bounded cross-section of C ^1,1-class. We study resolvent estimates and maximal regularity of the Stokes operator in for 1 < q, r < ∞ and for arbitrary Muckenhoupt weights ω ∈ A _r with respect to x′ ∈ Σ. The proofs use an operator-valued Fourier multiplier theorem and techniques of unconditional Schauder decompositions based on the -boundedness of the family of solution operators for a system in Σ parametrized by the phase variable of the one-dimensional partial Fourier transform. Supported by the Gottlieb Daimler- und Karl Benz-Stiftung, grant no. S025/02-10/03.     

Full Measure Reducibility for Generic One-parameter Family of Quasi-periodic Linear Systems

Let be the set of m × m matrices A(λ) depending analytically on a parameter λ in a closed interval . Consider one-parameter families of quasi-periodic linear differential equations: , where is analytic and sufficiently small. We prove that there is an open and dense set in , such that for each the equation can be reduced to an equation with constant coefficients by a quasi-periodic linear transformation for almost all in Lebesgue measure sense provided that g is sufficiently small. The result gives an affirmative answer to a conjecture of Eliasson (In: Proceeding of Symposia in Pure Mathematics). Dedicated to Professor Zhifen Zhang on the occasion of her 80th birthday     

The Semigeostrophic Equations Discretized in Reference and Dual Variables

We study the evolution of a system of n particles in . That system is a conservative system with a Hamiltonian of the form , where W ₂ is the Wasserstein distance and μ is a discrete measure concentrated on the set . Typically, μ(0) is a discrete measure approximating an initial L ^∞ density and can be chosen randomly. When d  =  1, our results prove convergence of the discrete system to a variant of the semigeostrophic equations. We obtain that the limiting densities are absolutely continuous with respect to the Lebesgue measure. When converges to a measure concentrated on a special d–dimensional set, we obtain the Vlasov–Monge–Ampère (VMA) system. When, d = 1 the VMA system coincides with the standard Vlasov–Poisson system.     

Non-autonomous 2D Navier–Stokes System with Singularly Oscillating External Force and its Global Attractor

We study the global attractor of the non-autonomous 2D Navier–Stokes (N.–S.) system with singularly oscillating external force of the form . If the functions g ₀(x, t) and g ₁ (z, t) are translation bounded in the corresponding spaces, then it is known that the global attractor is bounded in the space H, however, its norm may be unbounded as since the magnitude of the external force is growing. Assuming that the function g ₁ (z, t) has a divergence representation of the form where the functions (see Section 3), we prove that the global attractors of the N.–S. equations are uniformly bounded with respect to for all . We also consider the “limiting” 2D N.–S. system with external force g ₀(x, t). We have found an estimate for the deviation of a solution of the original N.–S. system from a solution u ⁰(x, t) of the “limiting” N.–S. system with the same initial data. If the function g ₁ (z, t) admits the divergence representation, the functions g ₀(x, t) and g ₁ (z, t) are translation compact in the corresponding spaces, and , then we prove that the global attractors converges to the global attractor of the “limiting” system as in the norm of H. In the last section, we present an estimate for the Hausdorff deviation of from of the form: in the case, when the global attractor is exponential (the Grashof number of the “limiting” 2D N.–S. system is small).      

The Boundary Riemann Solver Coming from the Real Vanishing Viscosity Approximation

We study the limit of the hyperbolic–parabolic approximation

Rheological behavior of poly(methyl methacrylate)/polystyrene (PMMA/PS) blends with the addition of PMMA-<Emphasis Type="BoldItalic">ran</Emphasis>-PS

In this work, the dynamic behavior of poly(methyl methacrylate)/polystyrene blend to which P(S_0.5-ran-MMA_0.5) was added was studied. Several blend (ranging from 5 to 20 wt% of dispersed phase) and copolymer (up to 20 wt% with respect to dispersed phase) concentrations were studied. The rheological behavior of the blends was compared to Bousmina’s (Rheol Acta 38:73–83, 1999) and Palierne’s (Rheol Acta 29:204–214, 1990) generalized models. The relaxation spectra of the blends were also inferred, and the results were analyzed in light of the analysis of Jacobs et al. [J Rheol 43:1495–1509, 1999]. The relaxation spectra of the blends with smaller dispersed phase (below 10 wt%) and larger copolymer concentrations (above 0.4 wt%) showed the presence of four relaxation times, two corresponding to the blend phases, τ _F , corresponding to the relaxation of the shape of the dispersed phase of the blend and that can be attributed to the relaxation of Marangoni stresses tangential to the interface between the dispersed phase and matrix. The experimental values of and were used to infer the interfacial tension (Γ) and the interfacial complex shear modulus (β) for the different blends, Γ decreased with increasing copolymer concentration. β decreased with increasing blend dispersed phase concentration and decreasing copolymer concentration. The predictions of Palierne’s generalized model were found to corroborate the experimental data once the values of Γ and β, found analyzing the relaxation spectra, were used in the calculations. Bousmina’s model was found to corroborate the data only for larger dispersed phase concentration. Paper was presented at the 3rd Annual Rheology Conference, AERC 2006, April 27–29, 2006, Crete, Greece.     

Large amplitude oscillatory elongation flow

A filament stretching rheometer (FSR) was used for measuring the elongation flow with a large amplitude oscillative elongation imposed upon the flow. The large amplitude oscillation imposed upon the elongational flow as a function of the time t was defined as where ε is the Hencky strain, is a constant elongational rate for the base elongational flow, Λ the strain amplitude (Λ ≥ 0), and Ω the strain frequency. A narrow molecular mass distribution linear polystyrene with a molecular weight of 145 kg/mol was subjected to the oscillative flow. The onset of the steady periodic regime is reached at the same Hencky strain as the onset of the steady elongational viscosity ( Λ = 0). The integral molecular stress function formulation within the ‘interchain pressure’ concept agrees qualitatively with the experiments.     

A New Approach to the Fundamental Theorem of Surface Theory

The fundamental theorem of surface theory classically asserts that, if a field of positive-definite symmetric matrices (a _αβ ) of order two and a field of symmetric matrices (b _αβ ) of order two together satisfy the Gauss and Codazzi-Mainardi equations in a simply connected open subset ω of , then there exists an immersion such that these fields are the first and second fundamental forms of the surface , and this surface is unique up to proper isometries in . The main purpose of this paper is to identify new compatibility conditions, expressed again in terms of the functions a _αβ and b _αβ , that likewise lead to a similar existence and uniqueness theorem. These conditions take the form of the matrix equation

CO2 Injection in Geological Formations: Determining Macroscale Coefficients from Pore Scale Processes

Carbon dioxide (CO₂) injections in geological formations are usually performed for enhanced hydrocarbon recovery in oil and gas reservoirs and storage and sequestration in saline aquifers. Once CO₂ is injected into the formation, it propagates in the porous rock by dispersion and convection. Chemical reactions between brine ions and CO₂ molecules and consequent reactions with mineral grains are also important processes. The dynamics of CO₂ molecules in random porous media are modeled with a set of differential equations corresponding to pore scale and continuum macroscale. On the pore scale, convective–dispersive equation is solved considering reactions on the inner boundaries in a unit cell. A unit cell is the smallest portion of a porous media that can reproduce the porous media by repetition. Inner boundaries in a unit cell are the surfaces of the mineral grains. Dispersion process at the pore scale is transformed into continuum macroscale by adopting periodic boundary conditions for contiguous unit cells and applying Taylor-Aris dispersion theory known as macrotransport theory. Using this theory, the discrete porous system changes into a continuum system within which the propagation and interaction of CO₂ molecules with fluid and solid matrix of the porous media are characterized by three position-independent macroscopic coefficients: the mean velocity vector , dispersivity dyadic , and mean volumetric CO₂ depletion coefficient .     

The flow regime in the turbulent near wake of a hemisphere

The work presented is a wind tunnel study of the near wake region behind a hemisphere immersed in three different turbulent boundary layers. In particular, the effect of different boundary layer profiles on the generation and distribution of near wake vorticity and on the mean recirculation region is examined. Visualization of the flow around a hemisphere has been undertaken, using models in a water channel, in order to obtain qualitative information concerning the wake structure.List of symbols C _p pressure coefficient, - D diameter of hemisphere - n vortex shedding frequency - p pressure on model surface - p ₀ static pressure - Re Reynolds number, - St Strouhal number, - U, V, W local mean velocity components - mean freestream velocity inX direction - U ^* shear velocity, - u, v, w velocity fluctuations inX, Y andZ directions - X Cartesian coordinate in longitudinal direction - Y Cartesian coordinate in lateral direction - Z Cartesian coordinate in direction perpendicular to the wall - it* boundary layer displacement thickness, - diameter of model surface roughness - elevation angleI - O boundary layer momentum thickness, - _w wall shearing stress - dynamic viscosity of fluid - density of fluid - streamfunction - _x longitudinal component of vorticity, - _y lateral component of vorticity, - _z vertical component of vorticity, This paper was presented at the Ninth symposium on turbulence, University of Missouri-Rolla, October 1–3, 1984     

On the Non-Isentropic Sound Waves Propagation in a Cylindrical Tube Filled with Saturated Porous Media

A rigid frame, cylindrical capillary theory of sound propagation in porous media that includes the nonlinear effects of the Forchheimer type is laid out by using variational solutions. It is shown that the five main parameters governing the propagation of sound waves in a fluid contained in rigid cylindrical tubes filled with a saturated porous media are: the shear wave number, , the reduced frequency parameter, , the porosity, ε, Darcy number, , and Forchheimer number, . The manner in which the flow influences the attenuation and the phase velocities of the forward and backward propagating non-isentropic acoustic waves is deduced. It is found that the inclusion of the solid matrix increases wave’s attenuations and phase velocities for both forward and backward sound waves, while increasing the porosity and the reduced frequency number decreased attenuation and increased phase velocities. The effect of the steady flow is found to decrease the attenuation and phase velocities for forward sound waves, and enhance them for the backward sound waves. This work is done during a sabbatical leave year granted form the University of Jordan to Dr. Hamzeh Duwairi for the academic year 2007/2008 at the German Jordanian University.     

The MSF model: relation of nonlinear parameters to molecular structure of long-chain branched polymer melts

The elongational viscosity data of model PS combs (Hepperle J, Einfluss der Molekularen Struktur auf Rheologische Eigenschaften von Polystyrol- und Polycarbonatschmelzen. Doctoral Thesis, University Erlangen-Nürnberg, 2003) are reconsidered by including the interchain pressure term of Marrucci and Ianniruberto [Macromolecules 37:3934–3942, 2004] in the Molecular Stress Function model [Wagner et al., J Rheol 47(3):779–793, 2003, Wagner et al., J Rheol 49:1317–1327, 2005d]. Two nonlinear model parameters are needed to describe elongational flow, β and . The parameterβ determines the slope of the elongational viscosity after the inception of strain hardening. It is directly related to the molecular structure of the polymer and represents the ratio of the molar mass of the (branched) polymer to the molar mass of the backbone alone. β follows from the hypothesis of Wagner et al. [J Rheol 47(3):779–793, 2003] that side chains are compressed onto the backbone. We consider also the case that side chains are oriented by deformation, but not stretched, and found little difference in the model predictions. The parameter represents the maximum strain energy stored in the polymeric system and determines the steady-state value of the viscosity in extensional flows. The relation of this energy parameter to the molecular structure is discussed. Good correlations between the energy parameter and different coil contraction ratios, as determined either experimentally or calculated theoretically by considering the topology of the macromolecule, are found. The smaller the relative size of the polymer coil, the larger is the energy parameter and the more strain energy can be stored in the polymeric system. Presented at the 3rd Annual European Rheology Conference, AERC2006, Crete, Greece.