Strichartz Estimates and Moment Bounds for the Relativistic Vlasov–Maxwell System

We consider the relativistic Vlasov–Maxwell system with data of unrestricted size and without compact support in momentum space. In the two-dimensional and the two-and-a-half-dimensional cases, Glassey–Schaeffer proved (Commun Math Phys 185:257–284, 1997; Arch Ration Mech Anal 141:331–354, 1998; Arch Ration Mech Anal. 141:355–374, 1998) that for regular initial data with compact momentum support this system has unique global in time classical solutions. In this work we do not assume compact momentum support for the initial data and instead require only that the data have polynomial decay in momentum space. In the two-dimensional and the two-and-a-half-dimensional cases, we prove the global existence, uniqueness and regularity for solutions arising from this class of initial data. To this end we use Strichartz estimates and prove that suitable moments of the solution remain bounded. Moreover, we obtain a slight improvement of the temporal growth of the \({L^\infty_x}\) norms of the electromagnetic fields compared to Glassey and Schaeffer (Commun Math Phys 185:257–284, 1997; Arch Ration Mech Anal 141:355–374, 1998). In the three-dimensional case, we apply Strichartz estimates and moment bounds to show that a regular solution can be extended as long as \({{\|p_0^{\theta} f \|_{L^{q}_{x}L^1_{p}}}}\) remains bounded for \({\theta > \frac{2}{q}}\), \({2 < q \leqq \infty}\). This improves previous results of Pallard (Indiana Univ Math J 54(5):1395–1409, 2005; Commun Math Sci 13(2):347–354, 2015).     

Approximate Semi-analytical Solution for Injection–Falloff–Production Well Test: An Analytical Tool for the In Situ Estimation of Relative Permeability Curves

An injection–falloff–production test (IFPT) was originally proposed in Chen et al. (in: SPE conference paper, 2006. doi: 10.2118/103271-MS, SPE Reserv Eval Eng 11(1):95–107, 2008) as a well test for the in situ estimation of two-phase relative permeability curves to be used for simulating multiphase flows in porous media. Hence, we develop an approximate semi-analytical solution for the two-phase saturation distribution in an oil–water system during the flowback period of an IFPT according to the mathematical theory of waves. In fact, we show that the weak solution we construct for the saturation equation for the flowback period satisfies the Oleinik entropy condition and hence is unique. In addition, we allow the governing relative permeabilities during the flowback period to be different from the relative permeabilities during injection. Using the saturation solution with the steady-state pressure theory of Thompson and Reynolds, we obtain a solution for the wellbore pressure during the flowback period. By comparing results from our solution with those from a commercial numerical simulator, we show that our approximate semi-analytical solution yields accurate saturation profiles and bottom hole pressures history. The use of very small time steps and a highly refined radial grid is necessary to generate a good solution from a reservoir simulator. The approximate analytical pressure solution developed is used as a forward model to match pressure and water flow rate data from an IFPT in order to estimate reservoir rock absolute permeability and skin factor in conjunction with in situ imbibition and drainage water–oil relative permeabilities.     

On a Nonlinear Model for Tumor Growth in a Cellular Medium

We investigate the dynamics of a nonlinear model for tumor growth within a cellular medium. In this setting the “tumor” is viewed as a multiphase flow consisting of cancerous cells in either proliferating phase or quiescent phase and a collection of cells accounting for the “waste” and/or dead cells in the presence of a nutrient. Here, the tumor is thought of as a growing continuum \(\Omega \) with boundary \(\partial \Omega \) both of which evolve in time. In particular, the evolution of the boundary \(\partial \Omega \) is prescibed by a given velocity \({{{\varvec{V}}}.}\) The key characteristic of the present model is that the total density of cancerous cells is allowed to vary, which is often the case within cellular media. We refer the reader to the articles (Enault in Mathematical study of models of tumor growth, 2010; Li and Lowengrub in J Theor Biol, 343:79–91, 2014) where compressible type tumor growth models are investigated. Global-in-time weak solutions are obtained using an approach based on penalization of the boundary behavior, diffusion, viscosity and pressure in the weak formulation, as well as convergence and compactness arguments in the spirit of Lions (Mathematical topics in fluid dynamics. Compressible models, 1998) [see also Donatelli and Trivisa (J Math Fluid Mech 16: 787–803, 2004), Feireisl (Dynamics of viscous compressible fluids, 2014)].     

On Reducibility of 1d Wave Equation with Quasiperiodic in Time Potentials

In this paper we use a KAM theorem of Grébert and Thomann (Commun Math Phys 307:383–427, 2011) to prove the reducibility of the 1d wave equation with Dirichlet boundery conditions on \([0,\pi ]\) with a quasi-periodic in time potential under some symmetry assumptions. From Mathieu–Hill operator’s known results (Eastham in The spectral theory of periodic differential operators, Hafner, New York, 1974; Magnus and Winkler in Hill’s equation, Wiley-Interscience, London, 1969) and Bourgain’s techniques (Commun Math Phys 204:207–247, 1999), we prove that for any \(\epsilon \) small enough, there exist a \(0<m_{\epsilon }\le 1\) and one solution \(u_{\epsilon }(t,x)\) with

$$\begin{aligned} \Vert u_{\epsilon }(t_n,x)\Vert _{H^1({\mathbb {T}})}\rightarrow \infty , \qquad |t_n|\rightarrow \infty , \end{aligned}$$

where \(u_{\epsilon }(t,x)\) satisfies 1d wave equation

$$\begin{aligned} u_{tt}-u_{xx}+m_{\epsilon }u-\epsilon \cos 2t u=0, \end{aligned}$$

with Dirichlet boundery conditions on \([0,\pi ]\).

     

Bail-Down Test Simulation at Laboratory Scale

This paper presents a comparison of hydraulic oil conductivity obtained from interpreting bail-down test data to values calculated from theory. The bail-down tests were performed at laboratory scale, on a radial portion of a circular domain filled with calibrated sand allowing hydraulic oil conductivity to be calculated using Parker’s theoretical model (Parker et al. in Water Resour Res 23(4):618–624, 1987). The bail-down tests were interpreted using the modified Bouwer and Rice (Huntley in Ground Water 38(1):46–52, 2000) and the modified Cooper methods (Beckett and Lyverse in API Interact LNAPL Guide 2:1–27, 2002). The results show that (1) both interpretation methods from bail-down test data give similar hydraulic oil conductivities, and (2) the hydraulic oil conductivities estimated from bail-down test data agree well with the hydraulic oil conductivity predicted when using the Parker theoretical model. Overall, this paper confirms that the modified Bouwer and Rice (Huntley 2000) and the modified Cooper methods (Beckett and Lyverse 2002) are valid to estimate hydraulic oil conductivity, giving realistic values despite test conditions not meeting all the assumptions and boundary conditions of each analytical solution.     

On the Validation of a Compositional Model for the Simulation of $$\text {CO}_2$$ Injection into Saline Aquifers

In this work we apply a recently proposed Bayesian Markov chain Monte Carlo framework (Akbarabadi et al. in Comput Geosci 19(6):1231–1250, 2015) to quantify uncertainty in the three-dimensional permeability field of a rock core. This process establishes the credibility of a compositional two-phase flow model to describe the displacement of brine by \(\text {CO}_2\) and \(\text {CO}_2\) storage in saline aquifers. We investigate the predictive capabilities of the compositional model in the context of an unsteady-state \(\text {CO}_2\)-brine drainage experiment at the laboratory scale, performed at field-scale aquifer conditions. We employ forward models consisting of a system of discretized partial differential equations along with relative permeability curves obtained by a curve fitting of experimental measurements. We consider a forward model to be validated when: (1) numerical simulations reveal that the Bayesian framework has accurately characterized the core’s permeability and (2) Monte Carlo predictions show excellent agreement between measured and simulated data. A large set of numerical studies with an accurate compositional simulator shows that forward models have been successfully validated. For such models, our numerical results show that we are able to capture all the dominant features and general trends of the \(\text {CO}_2\) saturation fields observed in the core. Our study is consistent with the design and findings of real experiments. Fluid properties, relative permeability data, measured porosity field, physical dimensions, and thermodynamic conditions are the same as those reported in Akbarabadi and Piri (Adv Water Resour 52:190–206, 2013). However, the measured saturation data are from flow experiments different from those reported in Akbarabadi and Piri (2013), and will be presented here.     

Topological and Measure-Theoretical Entropies of Nonautonomous Dynamical Systems

We consider the dynamics of a nonautonomous dynamical system determined by a sequence of continuous self-maps \(f_n:X \rightarrow X,\) where \( n \in {\mathbb {N}},\) defined on a compact metric space X. Applying the theory of the Carathéodory structures, elaborated by Pesin (Dimension Theory in Dynamical Systems. Chicago Lectures in Mathematics. The University of Chicago Press, Chicago, 1997), we construct a Carathéodory structure whose capacity coincides with the topological entropy of the considered system. Generalizing the notion of local measure entropy, introduced by Brin and Katok (in: Palis (ed) Geometric Dynamics, Lecture Notes in Mathematics. Springer, Berlin 1983) for a single map, to a nonautonomous dynamical system we provide a lower and upper estimations of the topological entropy by local measure entropies. The theorems of the paper generalize results of Kawan (Nonautonomous Stoch Dyn Syst 1:26–52, 2013) and of Feng and Huang (J Funct Anal 263:2228–2254, 2012). Also, we construct a new entropy-like invariant such the entropy of a sequence \(\{f_n:X \rightarrow X\}_{n=1}^{\infty }\) of Lipschitz continuous maps with the same Lipschitz constant \(L >1,\) restricted to a subset \(Y\subset X,\) is upper bounded by Hausdorff dimension of Y multiplied by the logarithm of the Lipschitz constant L. This gives a generalizations of results of Dai et al. (Sci China Ser A 41:1068–1075, 1998) and Misiurewicz (Discret Contin Dyn Syst 10:827–833, 2004).     

Two-subsystem thermodynamics for the mechanics of aging amorphous solids

The effect of physical aging on the mechanics of amorphous solids as well as mechanical rejuvenation is modeled with nonequilibrium thermodynamics, using the concept of two thermal subsystems, namely a kinetic one and a configurational one. Earlier work (Semkiv and Hütter in J Non-Equilib Thermodyn 41(2):79–88, 2016) is extended to account for a fully general coupling of the two thermal subsystems. This coupling gives rise to hypoelastic-type contributions in the expression for the Cauchy stress tensor, that reduces to the more common hyperelastic case for sufficiently long aging. The general model, particularly the reversible and irreversible couplings between the thermal subsystems, is compared in detail with models in the literature (Boyce et al. in Mech Mater 7:15–33, 1988; Buckley et al. in J Mech Phys Solids 52:2355–2377, 2004; Klompen et al. in Macromolecules 38:6997–7008, 2005; Kamrin and Bouchbinder in J Mech Phys Solids 73:269–288 2014; Xiao and Nguyen in J Mech Phys Solids 82:62–81, 2015). It is found that only for the case of Kamrin and Bouchbinder (J Mech Phys Solids 73:269–288, 2014) there is a nontrivial coupling between the thermal subsystems in the reversible dynamics, for which the Jacobi identity is automatically satisfied. Moreover, in their work as well as in Boyce et al. (Mech Mater 7:15–33, 1988), viscoplastic deformation is driven by the deviatoric part of the Cauchy stress tensor, while for Buckley et al. (J Mech Phys Solids 52:2355–2377, 2004) and Xiao and Nguyen (J Mech Phys Solids 82:62–81, 2015) this is not the case.     

Statistical and Numerical Considerations of Backus-Average Product Approximation

In this paper, we examine the applicability of the approximation, \(\overline{f g}\approx \overline{f}\,\overline{g}\), within Backus (J. Geophys. Res. 67(11):4427–4440, 1962) averaging. This approximation is a crucial step in the method proposed by Backus (J. Geophys. Res. 67(11):4427–4440, 1962), which is widely used in studying wave propagation in layered Hookean solids. According to this approximation, the average of the product of a rapidly varying function and a slowly varying function is approximately equal to the product of the averages of those two functions.Considering that the rapidly varying function represents the mechanical properties of layers, we express it as a step function. The slowly varying function is continuous, since it represents the components of the stress or strain tensors. In this paper, beyond the upper bound of the error for that approximation, which is formulated by Bos et al. (J. Elast. 127:179–196, 2017), we provide a statistical analysis of the approximation by allowing the function values to be sampled from general distributions.Even though, according to the upper bound, Backus (J. Geophys. Res. 67(11):4427–4440, 1962) averaging might not appear as a viable approach, we show that—for cases representative of physical scenarios modelled by such an averaging—the approximation is typically quite good. We identify the cases for which there can be a deterioration in its efficacy.In particular, we examine a special case for which the approximation results in spurious values. However, such a case—though physically realizable—is not likely to appear in seismology, where Backus (J. Geophys. Res. 67(11):4427–4440, 1962) averaging is commonly used. Yet, such values might occur in material sciences, in general, for which Backus (J. Geophys. Res. 67(11):4427–4440, 1962) averaging is also considered.     

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.     

The peculiar elongational viscosity of concentrated solutions of monodisperse PMMA in oligomeric MMA

Concentrated solutions of nearly monodisperse poly(methyl methacrylate), PMMA-270k and PMMA-86k, in oligo(methyl methacrylate), MMA o-4k and MMA o-2k, investigated by Wingstrand et al. (Phys Rev Lett 115:078302, 2015) and Wingstrand (2015) do not follow the linear-viscoelastic scaling relations of monodisperse polystyrenes (PS) dissolved in oligomeric styrene (Wagner in Rheol Acta 53:765–777, 2014a, in Non-Newtonian Fluid Mech 222:121–131, 2014b; Wagner et al. in J Rheol 59:1113–1130, 2015). Rather, PMMA-270k shows an attractive interaction with MMA, in contrast to the interaction of PMMA-86k and MMA. This different behavior can be traced back to different tacticities of the two polymers. The attractive interaction of PMMA-270k with o-4k creates pseudo entanglements, which increase the interchain tube pressure, and therefore, the solution PMMA-270k/o-4k shows, as reported by Wingstrand et al. (Phys Rev Lett 115:078302, 2015), qualitatively a similar scaling of the elongational viscosity with \( {\left(\dot{\varepsilon}{\tau}_R\right)}^{-1/2} \) as observed for polystyrene melts. For the solution PMMA-270/o-2k, this effect is only seen at the highest elongation rates investigated. The elongational viscosity of PMMA-86k dissolved in oligomeric MMA is determined by the Rouse time of the melt, as in the case of polystyrene solutions.

Open image in new window

Graphical abstract ?

     

Thermodynamic model formulation for viscoplastic solids as general equations for non-equilibrium reversible–irreversible coupling

Thermodynamic models for viscoplastic solids are often formulated in the context of continuum thermodynamics and the dissipation principle. The purpose of the current work is to show that models for such material behavior can also be formulated in the form of a General Equation for Non-Equilibrium Reversible–Irreversible Coupling (GENERIC), see, e.g., Grmela and Öttinger (Phys Rev E, 56:6620–6632, 1997), Öttinger and Grmela (Phys Rev E, 56:6633–6655, 1997), Grmela (J Non-Newtonian Fluid Mech, 165:980–986, 2010). A GENERIC combines Hamiltonian-dynamics-based modeling of time-reversible processes with Onsager–Casimir-based modeling of time-irreversible processes. The result is a model for the approach of non-equilibrium systems to thermodynamic equilibrium. Originally developed to model complex fluids, it has recently been applied to anisotropic inelastic solids in Eulerian (Hütter and Tervoort, in J Non-Newtonian Fluid Mech, 152:45–52, 2008; Hütter and Tervoort, in J Non-Newtonian Fluid Mech, 152:53–65, 2008; Hütter and Tervoort, in Adv Appl Mech, 42:254–317, 2008) and Lagrangian (Hütter and Svendsen, in J Elast 104:357–368, 2011) settings, as well as to damage mechanics. For simplicity, attention is focused in the current work on the case of thermoelastic viscoplasticity. Central to this formulation is a GENERIC-based form for the viscoplastic flow rule. A detailed comparison with the formulation based on continuum thermodynamics and the dissipation principle is given.     

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.     

On the Generalization of Reissner Plate Theory to Laminated Plates,Part I: Theory

This is the first part of a two-part paper presenting the generalization of Reissner thick plate theory (Reissner in J. Math. Phys. 23:184–191, 1944) to laminated plates and its relation with the Bending-Gradient theory (Lebée and Sab in Int. J. Solids Struct. 48(20):2878–2888, 2011 and in Int. J. Solids Struct. 48(20):2889–2901, 2011). The original thick and homogeneous plate theory derived by Reissner (J. Math. Phys. 23:184–191, 1944) is based on the derivation of a statically compatible stress field and the application of the principle of minimum of complementary energy. The static variables of this model are the bending moment and the shear force. In the present paper, the rigorous extension of this theory to laminated plates is presented and leads to a new plate theory called Generalized-Reissner theory which involves the bending moment, its first and second gradients as static variables. When the plate is homogeneous or functionally graded, the original theory from Reissner is retrieved. In the second paper (Lebée and Sab, 2015), the Bending-Gradient theory is obtained from the Generalized-Reissner theory and comparison with an exact solution for the cylindrical bending of laminated plates is presented.     

Second Order Corrections to Mean Field Evolution for Weakly Interacting Bosons in The Case of Three-body Interactions

In this paper, we consider the Hamiltonian evolution of N weakly interacting bosons. Assuming triple collisions, its mean field approximation is given by a quintic Hartree equation. We construct a second order correction to the mean field approximation using a kernel k(t, x, y) and derive an evolution equation for k. We show global existence for the resulting evolution equation for the correction and establish an a priori estimate comparing the approximation to the exact Hamiltonian evolution. Our error estimate is global and uniform in time. Comparing with the work of Rodnianski and Schlein (Commun Math Phys 291:31–61, 2009), and Grillakis, Machedon and Margetis (Commun Math Phys 294:273–301, 2010; Adv Math 288:1788–1815, 2011), where the error estimate grows in time, our approximation tracks the exact dynamics for all time with an error of the order \({O(1/\sqrt{N}).}\)     

On the Regularity Set and Angular Integrability for the Navier–Stokes Equation

We investigate the size of the regular set for suitable weak solutions of the Navier–Stokes equation, in the sense of Caffarelli–Kohn–Nirenberg (Commun Pure Appl Math 35:771–831, 1982). We consider initial data in weighted Lebesgue spaces with mixed radial-angular integrability, and we prove that the regular set increases if the data have higher angular integrability, invading the whole half space \({\{t > 0\}}\) in an appropriate limit. In particular, we obtain that if the \({L^{2}}\) norm with weight \({|x|^{-\frac12}}\) of the data tends to 0, the regular set invades \({\{t > 0\}}\); this result improves Theorem D of Caffarelli et al. (Commun Pure Appl Math 35:771–831, 1982).     

On the Generalization of Reissner Plate Theory to Laminated Plates,Part II: Comparison with the Bending-Gradient Theory

In the first part of this two-part paper (Lebée and Sab in On the generalization of Reissner plate theory to laminated plates, Part I: theory, doi: 10.1007/s10659-016-9581-6, 2015), the original thick plate theory derived by Reissner (J. Math. Phys. 23:184–191, 1944) was rigorously extended to the case of laminated plates. This led to a new plate theory called Generalized-Reissner theory which involves the bending moment, its first and second gradients as static variables. In this second paper, the Bending-Gradient theory (Lebée and Sab in Int. J. Solids Struct. 48(20):2878–2888, 2011 and 2889–2901, 2011) is obtained from the Generalized-Reissner theory and several projections as a Reissner–Mindlin theory are introduced. A comparison with an exact solution for the cylindrical bending of laminated plates is presented. It is observed that the Generalized-Reissner theory converges faster than the Kirchhoff theory for thin plates in terms of deflection. The Bending-Gradient theory does not converge faster but improves considerably the error estimate.     

The Effects of Brine Species on the Formation of Residual Water in a \hbox {CO}_{2}–Brine System

Geological storage of \(\hbox {CO}_{2}\) in deep saline aquifers is achieved by injecting \(\hbox {CO}_{2}\) into the aquifers and displacing the brine. Although most of the brine is displaced, some residual groundwater remains in the rock pores. We conducted experiments to investigate factors that influence how much of this residual water remains after \(\hbox {CO}_{2}\) is injected. A rock sample was saturated with brines of two different salts. Supercritical \(\hbox {CO}_{2}\) was injected into the samples at aquifer temperature and pressure, and the displaced water and water–gas mixtures were collected and measured. The results show that deionized water drains more completely than either of the two brines, and NaCl brine drains more completely than \(\hbox {CaCl}_{2}\) brine. The ranking of the irreducible water saturation at the end of the experiment is deionized \(\hbox {water}<\hbox {NaCl brine } <\hbox {CaCl}_{2}\) brine. The process of drainage can be divided into three stages according to the drainage flow rates; the Pushing Drainage, Portable Drainage, and Dissolved Drainage stages. This paper proposed a capillary model which is used to interpret the mechanisms that characterize these three stages.     

Partial Regularity of the Suitable Weak Solutions to the Multi-dimensional Incompressible Boussinesq Equations

This work is concerned with the partial regularity of the suitable weak solutions to the Boussinesq equations in \(\mathbb {R}^{n}\) where \(n=3,\,4\). By means of the De Giorgi iteration method developed in Vasseur (Nonlinear Differ Equ Appl 14(5–6):753–785, 2007), Wang, Wu (J Differ Equ 256(3):1224–1249, 2014), we obtain that \(n-2\) dimensional parabolic Hausdorff measure of the possible singular points set of the suitable weak solutions to this system is zero. Particularly, we obtain some interior regularity criteria only in terms of the scaled mixed norm of velocity for the suitable weak solutions to the Boussinesq equations, which implies that the potential singular points may only stem from the velocity field.     

Stability of Concatenated Traveling Waves

We consider a reaction–diffusion equation in one space dimension whose initial condition is approximately a sequence of widely separated traveling waves with increasing velocity, each of which is individually asymptotically stable. We show that the sequence of traveling waves is itself asymptotically stable: as \(t\rightarrow \infty \), the solution approaches the concatenated wave pattern, with different shifts of each wave allowed. Essentially the same result was previously proved by Wright (J Dyn Differ Equ 21:315–328, 2009) and Selle (Decomposition and stability of multifronts and multipulses, 2009), who regarded the concatenated wave pattern as a sum of traveling waves. In contrast to their work, we regard the pattern as a sequence of traveling waves restricted to subintervals of \(\mathbb {R}\) and separated at any finite time by small jump discontinuities. Our proof uses spatial dynamics and Laplace transform.