Transitional and Turbulent Flow in a Bed of Spheres as Measured with Stereoscopic Particle Image Velocimetry

Stereoscopic particle image velocimetry has been used to investigate inertia dominated, transitional and turbulent flow in a randomly packed bed of monosized PMMA spheres. By using an index-matched fluid, the bed is optically transparent and measurements can be performed in an arbitrary position within the porous bed. The velocity field observations are carried out for particle Reynolds numbers, \({Re}_{{p}}\), between 20 and 3220, and the sampling is done at a frequency of 75 Hz. Results show that, in porous media, the dynamics of the flow can vary significantly from pore to pore. At \({Re}_{{p}}\) around 400 the spatially averaged time fluctuations of total velocity reach a maximum and the spatial variation of the time-averaged total velocity, \(u_\mathrm{tot}\) increases up to about the same \({Re}_{{p}}\) and then it decreases. Also in the studied planes, a considerable amount of the fluid moves in the perpendicular directions to the main flow direction and the time-averaged magnitude of the velocity in the main direction, \(u_{x}\), has an averaged minimum of 40% of the magnitude of \(u_\mathrm{tot}\) at \({Re}_{{p}}\) about 400. For \({Re}_{{p}} > 1600\), this ratio is nearly constant and \(u_{x}\) is on average a little bit less than 50% of \(u_\mathrm{tot}\). The importance of the results for longitudinal and transverse dispersion is discussed.     

Approximation of Flat <Emphasis Type="Italic">W</Emphasis><Superscript>2,2</Superscript> Isometric Immersions by Smooth Ones

Let \({S\subset\mathbb{R}^2}\) be a bounded Lipschitz domain and denote by \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\) the set of mappings \({u\in W^{2,2}(S;\mathbb{R}^3)}\) which satisfy \({(\nabla u)^T(\nabla u) = Id}\) almost everywhere. Under an additional regularity condition on the boundary \({\partial S}\) (which is satisfied if \({\partial S}\) is piecewise continuously differentiable), we prove that the strong W ^2,2 closure of \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)\cap C^{\infty}(\overline{S};\mathbb{R}^3)}\) agrees with \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\).     

Insight into Heterogeneity Effects in Methane Hydrate Dissociation via Pore-Scale Modeling

A large number (1253) of high-quality streaming potential coefficient (\(C_\mathrm{sp})\) measurements have been carried out on Berea, Boise, Fontainebleau, and Lochaline sandstones (the latter two including both detrital and authigenic overgrowth forms), as a function of pore fluid salinity (\(C_\mathrm{f})\) and rock microstructure. All samples were saturated with fully equilibrated aqueous solutions of NaCl (10\(^{-5}\) and 4.5 mol/dm\(^{3})\) upon which accurate measurements of their electrical conductivity and pH were taken. These \(C_\mathrm{sp}\) measurements represent about a fivefold increase in streaming potential data available in the literature, are consistent with the pre-existing 266 measurements, and have lower experimental uncertainties. The \(C_\mathrm{sp}\) measurements follow a pH-sensitive power law behaviour with respect to \(C_\mathrm{f}\) at medium salinities (\(C_\mathrm{sp} =-\,1.44\times 10^{-9} C_\mathrm{f}^{-\,1.127} \), units: V/Pa and mol/dm\(^{3})\) and show the effect of rock microstructure on the low salinity \(C_\mathrm{sp}\) clearly, producing a smaller decrease in \(C_\mathrm{sp}\) per decade reduction in \(C_\mathrm{f}\) for samples with (i) lower porosity, (ii) larger cementation exponents, (iii) smaller grain sizes (and hence pore and pore throat sizes), and (iv) larger surface conduction. The \(C_\mathrm{sp}\) measurements include 313 made at \(C_\mathrm{f} > 1\) mol/dm\(^{3}\), which confirm the limiting high salinity \(C_\mathrm{sp}\) behaviour noted by Vinogradov et al., which has been ascribed to the attainment of maximum charge density in the electrical double layer occurring when the Debye length approximates to the size of the hydrated metal ion. The zeta potential (\(\zeta \)) was calculated from each \(C_\mathrm{sp}\) measurement. It was found that \(\zeta \) is highly sensitive to pH but not sensitive to rock microstructure. It exhibits a pH-dependent logarithmic behaviour with respect to \(C_\mathrm{f}\) at low to medium salinities (\(\zeta =0.01133 \log _{10} \left( {C_\mathrm{f} } \right) +0.003505\), units: V and mol/dm\(^{3})\) and a limiting zeta potential (zeta potential offset) at high salinities of \({\zeta }_\mathrm{o} = -\,17.36\pm 5.11\) mV in the pH range 6–8, which is also pH dependent. The sensitivity of both \(C_\mathrm{sp}\) and \(\zeta \) to pH and of \(C_\mathrm{sp}\) to rock microstructure indicates that \(C_\mathrm{sp}\) and \(\zeta \) measurements can only be interpreted together with accurate and equilibrated measurements of pore fluid conductivity and pH and supporting microstructural and surface conduction measurements for each sample.     

Global Well-Posedness for Navier–Stokes Equations with Small Initial Value in $${B^{0}_{n,\infty}(\Omega)}$$

We prove global well-posedness for instationary Navier–Stokes equations with initial data in Besov space \({B^{0}_{n,\infty}(\Omega)}\) in whole and half space, and bounded domains of \({{\mathbb R}^{n}}\), \({n \geq 3}\). To this end, we prove maximal \({L^{\infty}_{\gamma}}\) -regularity of the sectorial operators in some Banach spaces and, in particular, maximal \({L^{\infty}_{\gamma}}\) -regularity of the Stokes operator in little Nikolskii spaces \({b^{s}_{q,\infty}(\Omega)}\), \({s \in (-1, 2)}\), which are of independent significance. Then, based on the maximal regularity results and \({b^{s_{1}}_{q_{1},\infty}-B^{s_{2}}_{q_{2,1}}}\) estimates of the Stokes semigroups, we prove global well-posedness for Navier–Stokes equations under smallness condition on \({\|u_{0}\|_{B^{0}_{n,\infty}(\Omega)}}\) via a fixed point argument using Banach fixed point theorem.     

Displacement and Dissolution Characteristics of CO$$_{}$$/Brine System in Unconsolidated Porous Media

This study investigated the dynamic displacement and dissolution of \(\hbox {CO}_{2}\) in porous media at 313 K and 6/8 MPa. Gaseous (\(\hbox {gCO}_{2}\)) at 6 MPa and supercritical \(\hbox {CO}_{2 }(\hbox {scCO}_{2}) \) at 8 MPa were injected downward into a glass bead pack at different flow rates, following upwards brine injection. The processes occurring during \(\hbox {CO}_{2}\) drainage and brine imbibition were visualized using magnetic resonance imaging. The drainage flow fronts were strongly influenced by the flow rates, resulting in different gas distributions. However, brine imbibition proceeded as a vertical compacted front due to the strong effect of gravity. Additionally, the effects of flow rate on distribution and saturation were analyzed. Then, the front movement of \(\hbox {CO}_{2}\) dissolution was visualized along different paths after imbibition. The determined \(\hbox {CO}_{2}\) concentrations implied that little \(\hbox {scCO}_{2}\) dissolved in brine after imbibition. The dissolution rate was from \(10^{-8}\) to \(10^{-9}\, \hbox {kg}\, \hbox {m}^{-3} \, \hbox {s}^{-1}\) and from \(10^{-6}\) to \(10^{-8}\, \hbox {kg}\, \hbox {m}^{-3} \, \hbox {s}^{-1}\) for \(\hbox {gCO}_{2}\) at 6 MPa and \(\hbox {scCO}_{2 }\) at 8 MPa, respectively. The total time for the \(\hbox {scCO}_{2}\) dissolution was short, indicating fast mass transfer between the \(\hbox {CO}_{2}\) and brine. Injection of \(\hbox {CO}_{2}\) under supercritical conditions resulted in a quick establishment of a steady state with high storage safety.     

Uniqueness of Positive Ground State Solutions of the Logarithmic Schrödinger Equation

We prove the uniqueness of positive ground state solutions of the problem \({ {\frac {d^{2}u}{dr^{2}}} + {\frac {n-1}{r}}{\frac {du}{dr}} + u \ln(|u|) = 0}\), \({u(r) > 0~\forall r \ge 0}\), and \({(u(r),u'(r)) \to (0, 0)}\) as \({r \to \infty}\). This equation is derived from the logarithmic Schrödinger equation \({{\rm i}\psi_{t} = {\Delta} \psi + u \ln \left(|u|^{2}\right)}\), and also from the classical equation \({{\frac {\partial u}{\partial t}} = {\Delta} u +u \left(|u|^{p-1}\right) -u}\). For each \({n \ge 1}\), a positive ground state solution is \({ u_{0}(r) = \exp \left(-{\frac{r^2}{4}} + {\frac{n}{2}}\right),~0 \le r < \infty}\). We combine \({u_{0}(r)}\) with energy estimates and associated Ricatti equation estimates to prove that, for each \({n \in \left[1, 9 \right]}\), \({u_{0}(r)}\) is the only positive ground state. We also investigate the stability of \({u_{0}(r)}\). Several open problems are stated.     

On the Global Regularity of the Two-Dimensional Density Patch for Inhomogeneous Incompressible Viscous Flow

Regarding P.-L. Lions’ open question in Oxford Lecture Series in Mathematics and its Applications, Vol. 3 (1996) concerning the propagation of regularity for the density patch, we establish the global existence of solutions to the two-dimensional inhomogeneous incompressible Navier–Stokes system with initial density given by \({(1 - \eta){\bf 1}_{{\Omega}_{0}} + {\bf 1}_{{\Omega}_{0}^{c}}}\) for some small enough constant \({\eta}\) and some \({W^{k+2,p}}\) domain \({\Omega_{0}}\), with initial vorticity belonging to \({L^{1} \cap L^{p}}\) and with appropriate tangential regularities. Furthermore, we prove that the regularity of the domain \({\Omega_0}\) is preserved by time evolution.     

Linearly Stable Quasi-Periodic Breathers in a Class of Random Hamiltonian Systems

In this paper, we construct linearly stable quasi-periodic breathers for the Hamiltonian systems in the form \({{\rm i} \dot{q}_n+v_n q_n+\delta|q_n|^2q_n+\varepsilon_n \left(q_{n+1}+q_{n-1} \right)=0,\quad n \in \mathbb{Z}}\) where \({\{v_n\}_{n \in \mathbb{Z}}}\) is a family of time independent identically distributed (i.i.d) random variables with common distribution \({g = dv_n, v_n \in [0,1]}\) and \({|\varepsilon_n| \leq \varepsilon e^{-\varrho |n|}}\) with \({\varepsilon,\varrho > 0}\) . We prove that for \({\varepsilon, \delta}\) sufficiently small, the equation admits a family of small-amplitude and linearly stable, time quasi-periodic solutions for most of the parameters \({\{v_n\}_{n \in \mathbb{Z}}}\) .     

On Nonperiodic Euler Flows with Hölder Regularity

In (Isett, Regularity in time along the coarse scale flow for the Euler equations, 2013), the first author proposed a strengthening of Onsager’s conjecture on the failure of energy conservation for incompressible Euler flows with Hölder regularity not exceeding \({1/3}\). This stronger form of the conjecture implies that anomalous dissipation will fail for a generic Euler flow with regularity below the Onsager critical space \({L_t^\infty B_{3,\infty}^{1/3}}\) due to low regularity of the energy profile. This paper is the first and main paper in a series of two, the results of which may be viewed as first steps towards establishing the conjectured failure of energy regularity for generic solutions with Hölder exponent less than \({1/5}\). The main result of the present paper shows that any given smooth Euler flow can be perturbed in \({C^{1/5-\epsilon}_{t,x}}\) on any pre-compact subset of \({\mathbb{R}\times \mathbb{R}^3}\) to violate energy conservation. Furthermore, the perturbed solution is no smoother than \({C^{1/5-\epsilon}_{t,x}}\). As a corollary of this theorem, we show the existence of nonzero \({C^{1/5-\epsilon}_{t,x}}\) solutions to Euler with compact space-time support, generalizing previous work of the first author (Isett, Hölder continuous Euler flows in three dimensions with compact support in time, 2012) to the nonperiodic setting.     

Exponential Decay of the Vorticity in the Steady-State Flow of a Viscous Liquid Past a Rotating Body

We construct a Sobolev homeomorphism in dimension \({n \geqq 4,\,f \in W^{1,1}((0, 1)^n,\mathbb{R}^n)}\) such that \({J_f = {\rm det} Df > 0}\) on a set of positive measure and J _f  < 0 on a set of positive measure. It follows that there are no diffeomorphisms (or piecewise affine homeomorphisms) f _k such that \({f_k\to f}\) in \({W^{1,1}_{\rm loc}}\).     

On the Damped Harmonic Oscillator with Time Dependent Damping Coefficient

Conditions guaranteeing asymptotic stability for the differential equation

$$\begin{aligned} x''+h(t)x'+\omega ^2x=0 \qquad (x\in \mathbb {R}) \end{aligned}$$

are studied, where the damping coefficient \(h:[0,\infty )\rightarrow [0,\infty )\) is a locally integrable function, and the frequency \(\omega >0\) is constant. Our conditions need neither the requirement \(h(t)\le \overline{h}<\infty \) (\(t\in [0,\infty )\); \(\overline{h}\) is constant) (“small damping”), nor \(0< \underline{h}\le h(t)\) (\(t\in [0,\infty )\); \(\underline{h}\) is constant) (“large damping”); in other words, they can be applied to the general case \(0\le h(t)<\infty \) (\(t\in [0,\infty \))). We establish a condition which combines weak integral positivity with Smith’s growth condition

$$\begin{aligned} \int ^\infty _0 \exp [-H(t)]\int _0^t \exp [H(s)]\,\mathrm{{d}}s\,\mathrm{{d}}t=\infty \qquad \left( H(t):=\int _0^t h(\tau )\,\mathrm{{d}}\tau \right) , \end{aligned}$$

so it is able to control both the small and the large values of the damping coefficient simultaneously.

     

Visual Investigation of Improvement in Extra-Heavy Oil Recovery by Borate-Assisted $$\hbox {CO}_{2}$$-Foam Injection

The significant reduction in heavy oil viscosity when mixed with \(\hbox {CO}_{2}\) is well documented. However, for \(\hbox {CO}_{2}\) injection to be an efficient method for improving heavy oil recovery, other mechanisms are required to improve the mobility ratio between the \(\hbox {CO}_{2}\) front and the resident heavy oil. In situ generation of \(\hbox {CO}_{2}\)-foam can improve \(\hbox {CO}_{2}\) injection performance by (a) increasing the effective viscosity of \(\hbox {CO}_{2}\) in the reservoir and (b) increasing the contact area between the heavy oil and injected \(\hbox {CO}_{2}\) and hence improving \(\hbox {CO}_{2}\) dissolution rate. However, in situ generation of stable \(\hbox {CO}_{2}\)-foam capable of travelling from the injection well to the production well is hard to achieve. We have previously published the results of a series of foam stability experiments using alkali and in the presence of heavy crude oil (Farzaneh and Sohrabi 2015). The results showed that stability of \(\hbox {CO}_{2}\)-foam decreased by addition of NaOH, while it increased by addition of \(\hbox {Na}_{2}\hbox {CO}_{3}\). However, the highest increase in \(\hbox {CO}_{2}\)-foam stability was achieved by adding borate to the surfactant solution. Borate is a mild alkaline with an excellent pH buffering ability. The previous study was performed in a foam column in the absence of a porous medium. In this paper, we present the results of a new series of experiments carried out in a high-pressure glass micromodel to visually investigate the performance of borate–surfactant \(\hbox {CO}_{2}\)-foam injection in an extra-heavy crude oil in a transparent porous medium. In the first part of the paper, the pore-scale interactions of \(\hbox {CO}_{2}\)-foam and extra-heavy oil and the mechanisms of oil displacement and hence oil recovery are presented through image analysis of micromodel images. The results show that very high oil recovery was achieved by co-injection of the borate–surfactant solution with \(\hbox {CO}_{2}\), due to in-situ formation of stable foam. Dissolution of \(\hbox {CO}_{2}\) in heavy oil resulted in significant reduction in its viscosity. \(\hbox {CO}_{2}\)-foam significantly increased the contact area between the oil and \(\hbox {CO}_{2}\) significantly and thus the efficiency of the process. The synergy effect between the borate and surfactant resulted in (1) alteration of the wettability of the porous medium towards water wet and (2) significant reduction of the oil–water IFT. As a result, a bank of oil-in-water (O/W) emulsion was formed in the porous medium and moved ahead of the \(\hbox {CO}_{2}\)-foam front. The in-situ generated O/W emulsion has a much lower viscosity than the original oil and plays a major role in the observed additional oil recovery in the range of performed experiments. Borate also made \(\hbox {CO}_{2}\)-foam more stable by changing the system to non-spreading oil and reducing coalescence of the foam bubbles. The results of these visual experiments suggest that borate can be a useful additive for improving heavy oil recovery in the range of the performed tests, by increasing \(\hbox {CO}_{2}\)-foam stability and producing O/W emulsions.     

Asymmetric Domain Walls of Small Angle in Soft Ferromagnetic Films

We focus on a special type of domain wall appearing in the Landau–Lifshitz theory for soft ferromagnetic films. These domain walls are divergence-free \({\mathbb{S}^2}\)-valued transition layers that connect two directions \({m_\theta^\pm \in \mathbb{S}^2}\) (differing by an angle \({2\theta}\)) and minimize the Dirichlet energy. Our main result is the rigorous derivation of the asymptotic structure and energy of such “asymmetric” domain walls in the limit \({\theta \downarrow 0}\). As an application, we deduce that a supercritical bifurcation causes the transition from symmetric to asymmetric walls in the full micromagnetic model.     

Global Strong Well-Posedness of the Three Dimensional Primitive Equations in {L^p}-Spaces

In this article, an \({L^p}\)-approach to the primitive equations is developed. In particular, it is shown that the three dimensional primitive equations admit a unique, global strong solution for all initial data \({a \in [X_p,D(A_p)]_{1/p}}\) provided \({p \in [6/5,\infty)}\). To this end, the hydrostatic Stokes operator \({A_p}\) defined on \({X_p}\), the subspace of \({L^p}\) associated with the hydrostatic Helmholtz projection, is introduced and investigated. Choosing \({p}\) large, one obtains global well-posedness of the primitive equations for strong solutions for initial data \({a}\) having less differentiability properties than \({H^1}\), hereby generalizing in particular a result by Cao and Titi (Ann Math 166:245–267, 2007) to the case of non-smooth initial data.     

Resolvent Estimates for High-Contrast Elliptic Problems with Periodic Coefficients

We study the asymptotic behaviour of the resolvents \({(\mathcal{A}^\varepsilon+I)^{-1}}\) of elliptic second-order differential operators \({{\mathcal{A}}^\varepsilon}\) in \({\mathbb{R}^d}\) with periodic rapidly oscillating coefficients, as the period \({\varepsilon}\) goes to zero. The class of operators covered by our analysis includes both the “classical” case of uniformly elliptic families (where the ellipticity constant does not depend on \({\varepsilon}\)) and the “double-porosity” case of coefficients that take contrasting values of order one and of order \({\varepsilon^2}\) in different parts of the period cell. We provide a construction for the leading order term of the “operator asymptotics” of \({(\mathcal{A}^\varepsilon+I)^{-1}}\) in the sense of operator-norm convergence and prove order \({O(\varepsilon)}\) remainder estimates.     

On the Justification of Viscoelastic Elliptic Membrane Shell Equations

We consider a family of linearly viscoelastic shells with thickness \(2\varepsilon\), clamped along their entire lateral face, all having the same middle surface \(S=\boldsymbol{\theta}(\bar{\omega})\subset \mathbb{R}^{3}\), where \(\omega\subset\mathbb{R}^{2}\) is a bounded and connected open set with a Lipschitz-continuous boundary \(\gamma\). We make an essential geometrical assumption on the middle surface \(S\), which is satisfied if \(\gamma\) and \(\boldsymbol{\theta}\) are smooth enough and \(S\) is uniformly elliptic. We show that, if the applied body force density is \(O(1)\) with respect to \(\varepsilon\) and surface tractions density is \(O(\varepsilon)\), the solution of the scaled variational problem in curvilinear coordinates, \(\boldsymbol{u}( \varepsilon)\), defined over the fixed domain \(\varOmega=\omega\times (-1,1)\) for each \(t\in[0,T]\), converges to a limit \(\boldsymbol{u}\) with \(u_{\alpha}(\varepsilon)\rightarrow u_{\alpha}\) in \(W^{1,2}(0,T,H ^{1}(\varOmega))\) and \(u_{3}(\varepsilon)\rightarrow u_{3}\) in \(W^{1,2}(0,T,L^{2}(\varOmega))\) as \(\varepsilon\to0\). Moreover, we prove that this limit is independent of the transverse variable. Furthermore, the average \(\bar{\boldsymbol{u}}= \frac{1}{2}\int_{-1}^{1} \boldsymbol{u}dx_{3}\), which belongs to the space \(W^{1,2}(0,T, V_{M}( \omega))\), where

$$V_{M}(\omega)=H^{1}_{0}(\omega)\times H^{1}_{0}(\omega)\times L ^{2}(\omega), $$

satisfies what we have identified as (scaled) two-dimensional equations of a viscoelastic membrane elliptic shell, which includes a long-term memory that takes into account previous deformations. We finally provide convergence results which justify those equations.

     

On Regularity Criteria for the 2D Generalized MHD System

This paper deals with the problem of regularity criteria for the 2D generalized MHD system with fractional dissipative terms \({-\Lambda^{2\alpha}u}\) for the velocity field and \({-\Lambda^{2\beta}b}\) for the magnetic field respectively. Various regularity criteria are established to guarantee smoothness of solutions. It turns out that our regularity criteria imply previous global existence results naturally.     

Comparative Study of Five Outcrop Chalks Flooded at Reservoir Conditions: Chemo-mechanical Behaviour and Profiles of Compositional Alteration

This study presents experimental results from a flooding test series performed at reservoir conditions for five high-porosity Cretaceous onshore chalks from Denmark, Belgium and the USA, analogous to North Sea reservoir chalk. The chalks are studied in regard to their chemo-mechanical behaviour when performing tri-axial compaction tests while injecting brines (0.219 mol/L \(\hbox {MgCl}_{2}\) or 0.657 mol/L NaCl) at reservoir conditions for 2–3 months (T = 130 \(^\circ \hbox {C}\); 1 PV/d). Each chalk type was examined in terms of its mineralogical and chemical composition before and after the mechanical flooding tests, using an extensive set of analysis methods, to evaluate the chalk- and brine-dependent chemical alterations. All \(\hbox {MgCl}_{2}\)-flooded cores showed precipitation of Mg-bearing minerals (mainly magnesite). The distribution of newly formed Mg-bearing minerals appears to be chalk-dependent with varying peaks of enrichment. The chalk samples from Aalborg originally contained abundant opal-CT, which was dissolved with both NaCl and \(\hbox {MgCl}_{2}\) and partly re-precipitated as Si-Mg-bearing minerals. The Aalborg core injected with \(\hbox {MgCl}_{2}\) indicated strongly increased specific surface area (from 4.9 \(\hbox {m}^{2}\hbox {/g}\) to within 7–9 \(\hbox {m}^{2}\hbox {/g}\)). Mineral precipitation effects were negligible in chalk samples flooded with NaCl compared to \(\hbox {MgCl}_{2}\). Silicates were the main mineralogical impurity in the studied chalk samples (0.3–6 wt%). The cores with higher \(\hbox {SiO}_{2}\) content showed less deformation when injecting NaCl brine, but more compaction when injecting \(\hbox {MgCl}_{2}\)-brine. The observations were successfully interpreted by mathematical geochemical modelling which suggests that the re-precipitation of Si-bearing minerals leads to enhanced calcite dissolution and mass loss (as seen experimentally) explaining the high compaction seen in \(\hbox {MgCl}_{2}\)-flooded Aalborg chalk. Our work demonstrates that the original mineralogy, together with the newly formed minerals, can control the chemo-mechanical interactions during flooding and should be taken into account when predicting reservoir behaviour from laboratory studies. This study improves the understanding of complex flow reaction mechanisms also relevant for field-scale dynamics seen during brine injection.