Dynamical Parallelepipeds in Minimal Systems

For a topological dynamical system $(X,T)$ ( X , T ) and $d\in \mathbb N $ d ∈ N , the associated dynamical parallelepiped $\mathbf{Q}^{[d]}$ Q [ d ] was defined by Host–Kra–Maass. For a minimal distal system it was shown by them that the relation $\sim _{d-1}$ ～ d ? 1 defined on $\mathbf{Q}^{[d-1]}$ Q [ d ? 1 ] is an equivalence relation; the closing parallelepiped property holds, and for each $x\in X$ x ∈ X the collection of points in $\mathbf{Q}^{[d]}$ Q [ d ] with first coordinate $x$ x is a minimal subset under the face transformations. We give examples showing that the results do not extend to general minimal systems.     

A New Insight into the Convective Boundary Condition

The steady mixed convection boundary layer flows over a vertical surface adjacent to a Darcy porous medium and subject respectively to (i) a prescribed constant wall temperature, (ii) a prescribed variable heat flux, $q_\mathrm{w} =q_0 x^{-1/2}$ q w = q 0 x ? 1 / 2 , and (iii) a convective boundary condition are compared to each other in this article. It is shown that, in the characteristic plane spanned by the dimensionless flow velocity at the wall ${f}^{\prime }(0)\equiv \lambda $ f ′ ( 0 ) ≡ λ and the dimensionless wall shear stress $f^{\prime \prime }(0)\equiv S$ f ′ ′ ( 0 ) ≡ S , every solution $(\lambda , S)$ ( λ , S ) of one of these three flow problems at the same time is also a solution of the other two ones. There also turns out that with respect to the governing mixed convection and surface heat transfer parameters $\varepsilon $ ε and $\gamma $ γ , every solution $(\lambda , S)$ ( λ , S ) of the flow problem (iii) is infinitely degenerate. Specifically, to the very same flow solution $(\lambda , S)$ ( λ , S ) there corresponds a whole continuous set of values of $\varepsilon $ ε and $\gamma $ γ which satisfy the equation $S=-\gamma (1+\varepsilon -\lambda )$ S = ? γ ( 1 + ε ? λ ) . For the temperature solutions, however, the infinite degeneracy of the velocity solutions becomes lifted. These and further outstanding features of the convective problem (iii) are discussed in the article in some detail.     

Mixed Convection Boundary-Layer Flow on a Vertical Surface in a Porous Medium with a Constant Convective Boundary Condition

The mixed convection boundary-layer flow on a vertical surface heated convectively is considered when a constant surface heat transfer parameter is assumed. The problem is seen to be chararterized by a mixed convection parameter $\gamma $ γ . The flow and heat transfer near the leading edge correspond to forced convection solution and numerical solutions are obtained to determine how the solution then develops. The solution at large distances is obtained and this identifies a critical value $\gamma _c$ γ c of the parameter $\gamma $ γ . For $\gamma > \gamma _c$ γ > γ c a solution at large distances is possible and this is approached in the numerical integrations. For $\gamma <\gamma _c$ γ < γ c the numerical solution breaks down at a finite distance along the surface with a singularity, the nature of which is discussed.     

Invariant Parallels, Invariant Meridians and Limit Cycles of Polynomial Vector Fields on Some 2-Dimensional Algebraic Tori in \mathbb R ^3

We consider the polynomial vector fields of arbitrary degree in $\mathbb R ^3$ R 3 having the 2-dimensional algebraic torus $$\begin{aligned} \mathbb T ^2(l,m,n)=\{(x,y,z)\in \mathbb R ^3 : (x^{2l}+y^{2m}-r^2)^2+z^{2n}-1=0\}, \end{aligned}$$ T 2 ( l , m , n ) = { ( x , y , z ) ∈ R 3 : ( x 2 l + y 2 m - r 2 ) 2 + z 2 n - 1 = 0 } , where $l,m$ l , m , and $n$ n positive integers, and $r\in (1,\infty )$ r ∈ ( 1 , ∞ ) , invariant by their flow. We study the possible configurations of invariant meridians and parallels that these vector fields can exhibit on $\mathbb T ^2(l,m,n)$ T 2 ( l , m , n ) . Furthermore, we analyze when these invariant meridians or parallels are limit cycles.     

Local Dynamics at Focal Points Coupled with Elliptic Directions

We investigate local dynamics of a 4-dimensional system with small and slowly varying time periodic forcing. By assuming the unperturbed system is autonomous and has a fixed point with eigenvalues $(0,0,i,-i)$ ( 0 , 0 , i , ? i ) , we study homoclinic, subharmonic solutions and Hopf bifurcation in a $O(\epsilon )$ O ( ? ) neighborhood of the fixed point, where $\epsilon $ ? is the perturbation parameter.     

Basic Reproduction Ratios for Almost Periodic Compartmental Epidemic Models

The theory of the basic reproduction ratio $R_{0}$ R 0 and its computation formulae for almost periodic compartmental epidemic models are established. It is shown that the disease-free almost periodic solution is stable if $R_{0}<1$ R 0 < 1 , and unstable if $R_{0}>1$ R 0 > 1 . We also apply the developed theory to a patchy model with almost periodic population dispersal and disease transmission coefficients to obtain a threshold type result for uniform persistence and global extinction of the disease.     

Combined Effect of Stress, Pore Pressure and Temperature on Methane Permeability in Anthracite Coal: An Experimental Study

The effects of particle-size distribution on the longitudinal dispersion coefficient ( $D_{\mathrm{L}})$ D L ) in packed beds of spherical particles are studied by simulating a tracer column experiment. The packed-bed models consist of uniform and different-sized spherical particles with a ratio of maximum to minimum particle diameter in the range of 1–4. The modified version of Euclidian Voronoi diagrams is used to discretize the system of particles into cells that each contains one sphere. The local flow distribution is derived with the use of Laurent series. The flow pattern at low particle Reynolds number is then obtained by minimization of dissipation rate of energy for the dual stream function. The value of $D_{\mathrm{L}}$ D L is obtained by comparing the effluent curve from large discrete systems of spherical particles to the solution of the one-dimensional advection–dispersion equation. Main results are that at Peclet numbers above 1, increasing the width of the particle-size distribution increases the values of $D_{\mathrm{L}}$ D L in the packed bed. At Peclet numbers below 1, increasing the width of the particle-size distribution slightly lowers $D_{\mathrm{L}}$ D L .     

Inverse Approach in Ordinary Differential Equations: Applications to Lagrangian and Hamiltonian Mechanics

This paper is on the so called inverse problem of ordinary differential equations, i.e. the problem of determining the differential system satisfying a set of given properties. More precisely we characterize under very general assumptions the ordinary differential equations in \(\mathbb {R}^N\) which have a given set of either \(M\) partial integrals, or \(M first integral, or \(M partial and first integrals. Moreover, for such systems we determine the necessary and sufficient conditions for the existence of \(N-1\) independent first integrals. We give two relevant applications of the solutions of these inverse problem to constrained Lagrangian and Hamiltonian systems respectively. Additionally we provide the general solution of the inverse problem in dynamics.     

Linear Instability of a Horizontal Thermal Boundary Layer Formed by Vertical Throughflow in a Porous Medium: The Effect of Local Thermal Nonequilibrium

In this paper we investigate the onset of convection in a saturated porous medium where uniform suction into a horizontal and uniformly hot bounding surface induces a stationary thermal boundary layer. Particular attention is paid to how the well-known linear stability characteristics of this boundary layer are modified by the presence of local thermal nonequilibrium effects. The basic conduction state is determined and it is found that the boundary layer forms two distinct regions when the porosity is small or when the conductivity of the fluid is small compared with that of the solid. A linearised stability analysis is performed which results in an ordinary differential eigenvalue problem for the critical Darcy–Rayleigh number as a function of the wave number and the two nondimensional parameters, $H$ H and $\gamma $ γ , which are associated with local thermal nonequilibrium. This eigenvalue problem is solved numerically by first approximating the equations by fourth order compact finite differences, and then the critical Rayleigh number is computed iteratively using the inverse power method and minimised over the wavenumber. The variation of the critical Rayleigh number and wavenumber with $H$ H and $\gamma $ γ is presented. One of the unusual effects of local thermal nonequilibrium is that there exists a parameter regime within which the neutral curve is bimodal.     

Fracture Propagation Driven by Fluid Outflow from a Low-Permeability Aquifer

Deep saline aquifers are promising geological reservoirs for $\mathrm {CO}_2$ CO 2 sequestration if they do not leak. The absence of leakage is provided by the caprock integrity. However, $\mathrm {CO}_2$ CO 2 injection operations may change the geomechanical stresses and cause fracturing of the caprock. We present a model for the propagation of a fracture in the caprock driven by the outflow of fluid from a low-permeability aquifer. We show that to describe the fracture propagation, it is necessary to solve the pressure diffusion problem in the aquifer. We solve the problem numerically for the two-dimensional domain and show that, after a relatively short time, the solution is close to that of one-dimensional problem, which can be solved analytically. We use the relations derived in the hydraulic fracture literature to relate the width of the fracture to its length and the flux into it, which allows us to obtain an analytical expression for the fracture length as a function of time. Using these results we predict the propagation of a hypothetical fracture at the In Salah $\mathrm {CO}_2$ CO 2 injection site to be as fast as a typical hydraulic fracture. We also show that the hydrostatic and geostatic effects cause the increase of the driving force for the fracture propagation and, therefore, our solution serves as an estimate from below. Numerical estimates show that if a fracture appears, it is likely that it will become a pathway for $\mathrm {CO}_2$ CO 2 leakage.     

Conjugate Natural Convection in a Porous Enclosure Sandwiched by Finite Walls Under the Influence of Non-uniform Heat Generation and Radiation

Conjugate natural convection in a square porous enclosure sandwiched by finite walls under the influence of non-uniform heat generation and radiation is studied numerically in the present article. The horizontal heating is considered, where the vertical walls heated isothermally at different temperatures, while the horizontal walls are kept adiabatic. The Darcy model is used in the mathematical formulation for the porous layer and finite difference method is applied to solve the dimensionless governing equations. The governing parameters considered are the ratio of wall thickness to its width $(0.02 \le D \le 0.3)$ ( 0.02 ≤ D ≤ 0.3 ) , the wall to porous thermal conductivity ratio $(0.1 \le k_\mathrm{r} \le 10.0)$ ( 0.1 ≤ k r ≤ 10.0 ) , the internal heating $(0 \le \gamma \le 5)$ ( 0 ≤ γ ≤ 5 ) , and the local heating exponent parameters $(1 \le \lambda \le 20)$ ( 1 ≤ λ ≤ 20 ) . It is found that the average Nusselt number on the hot and cold interfaces increases with increasing the radiation intensity. Very high non-uniformity heating does not affect the average Nusselt number at very thick walls.     

Some Aspects of Stability for Semigroup Actions and Control Systems

The present paper introduces both the notions of Lagrange and Poisson stabilities for semigroup actions. Let \(S\) be a semigroup acting on a topological space \(X\) with mapping \(\sigma :S\times X\rightarrow X\) , and let \(\mathcal {F}\) be a family of subsets of \(S\) . For \(x\in X\) the motion \(\sigma _{x}:S\rightarrow X\) is said to be forward Lagrange stable if the orbit \(Sx\) has compact closure in \(X\) . The point \(x\) is forward \(\mathcal {F}\) -Poisson stable if and only if it belongs to the limit set \(\omega \left( x,\mathcal {F}\right) \) . The concept of prolongational limit set is also introduced and used to describe nonwandering points. It is shown that a point \(x\) is \( \mathcal {F}\) -nonwandering if and only if \(x\) lies in its forward \(\mathcal {F} \) -prolongational limit set \(J\left( x,\mathcal {F}\right) \) . The paper contains applications to control systems.     

A KAM Theorem for Higher Dimensional Nonlinear Schrödinger Equations

We prove an infinite dimensional KAM theorem. As an application, we use the theorem to study the higher dimensional nonlinear Schrödinger equation $$\begin{aligned} iu_t-\triangle u +M_\xi u+f(|u|^2)u=0, \quad t\in \mathbb{R }, x\in \mathbb{T }^d \end{aligned}$$ with periodic boundary conditions, where $M_\xi $ is a real Fourier multiplier and $f(|u|^2)$ is a real analytic function near $u=0$ with $f(0)=0$ . We obtain for the equation a Whitney smooth family of real-analytic small-amplitude linearly-stable quasi-periodic solutions with a nice linear normal form.     

The Stabilizing Effect of Spacetime Expansion on Relativistic Fluids With Sharp Results for the Radiation Equation of State

In this article, we study the 1 + 3-dimensional relativistic Euler equations on a pre-specified conformally flat expanding spacetime background with spatial slices that are diffeomorphic to ${\mathbb{R}^3.}$ R 3 . We assume that the fluid verifies the equation of state ${p = c^{2}_{s} \rho,}$ p = c s 2 ρ , where ${0 \leq c_{s} \leq \sqrt{1/3}}$ 0 ≤ c s ≤ 1 / 3 is the speed of sound. We also assume that the reciprocal of the scale factor associated with the expanding spacetime metric verifies a c _s ?dependent time-integrability condition. Under these assumptions, we use the vector field energy method to prove that an explicit family of physically motivated, spatially homogeneous, and spatially isotropic fluid solutions are globally future-stable under small perturbations of their initial conditions. The explicit solutions corresponding to each scale factor are analogs of the well-known spatially flat Friedmann–Lemaître–Robertson–Walker family. Our nonlinear analysis, which exploits dissipative terms generated by the expansion, shows that the perturbed solutions exist for all future times and remain close to the explicit solutions. This work is an extension of previous results, which showed that an analogous stability result holds when the spacetime is exponentially expanding. In the case of the radiation equation of state p = (1/3)ρ, we also show that if the time-integrability condition for the reciprocal of the scale factor fails to hold, then the explicit fluid solutions are unstable. More precisely, we show the existence of an open family of initial data such that (i) it contains arbitrarily small smooth perturbations of the explicit solutions’ data and (ii) the corresponding perturbed solutions necessarily form shocks in finite time. The shock formation proof is based on the conformal invariance of the relativistic Euler equations when ${c^{2}_{s} = 1/3,}$ c s 2 = 1 / 3 , which allows for a reduction to a well-known result of Christodoulou.     

Stable disarrangement phases of elastic aggregates: a setting for the emergence of no-tension materials with non-linear response in compression

A recent field theory of elastic bodies undergoing non-smooth submacroscopic geometrical changes (disarrangements) provides a setting in which, for a given homogeneous macroscopic deformation \(F\) of the body, there are typically a number of different states \(G\) of smooth, submacroscopic deformation (disarrangement phases) available to the body. A tensorial consistency relation and the inequality \(\det G\le \det F\) that guarantees that \(F\) accommodates \(G\) determine the totality of disarrangement phases \(G\) corresponding to \(F\) , and it is natural to seek for a given \(F\) those disarrangement phases that minimize the Helmholtz free energy (stable disarrangement phases). We introduce these concepts in the particular context of continuous bodies comprised of many small elastic bodies (elastic aggregates) and in the context where disarrangements do not contribute to the Helmholtz free energy (purely dissipative disarrangements). In this setting, the Helmholtz free energy response \(G\longmapsto \varPsi (G)\) of the pieces of the aggregate determines the totality of disarrangement phases corresponding to \(F\) , which necessarily includes the phase \(G=F\) (compact phase) in which every piece of the aggregate undergoes the given macroscopic deformation \(F\) . When the response function \(\varPsi \) is isotropic and smooth, and when \(\varPsi \) possesses standard semiconvexity and growth properties, the body also admits phases of the form \(G=\zeta _{\min }R\) (loose phases) with \(R\) an arbitrary rotation, provided that \(\zeta _{\min }R \) satisfies the accommodation inequality \(\zeta _{\min }^{3}\le \det F\) . Loose phases, when available, achieve the global minimum \(\varPsi (\zeta _{\min }R)\) of the free energy and consequently are stable and stress-free. When \( \varPsi (G)\) has the specific form \(\varPsi _{\alpha \beta }(G)=(\alpha /2)(\det G)^{-2}+(\beta /2)tr(GG^{T})\) , with \(\alpha \) , \(\beta \) given elastic constants, we determine all of the disarrangement phases corresponding to \(F\) . These include not only the compact and loose phases, but also disarrangement phases \(G\) in which the stress \(D\varPsi (G)\) is uniaxial or planar. Our main result (“stability implies no-tension”) is the assertion that every stable disarrangement phase for \(\varPsi _{\alpha \beta }\) cannot support tensile tractions, and our treatment of elastic aggregates thus provides a natural setting for the emergence of no-tension materials whose response in compression is non-linear. Existing treatments of no-tension materials assume at the outset that the body cannot support tension and that the response in compression is linear.     

Modeling and analysis of nonholonomic dynamic systems with a class of rheonomous affine constraints

This paper is devoted to modeling and theoretical analysis of dynamic control systems subject to a class of rheonomous affine constraints, which are called $A$ -rheonomous affine constraints. We first define $A$ -rheonomous affine constraints and explain their geometric representation. Next, a necessary and sufficient condition for complete nonholonomicity of $A$ -rheonomous affine constraints is shown. Then, we derive nonholonomic dynamic systems with $A$ -rheonomous affine constraints (NDSARAC), which are included in the class of nonlinear control systems. We also analyze linear approximated systems and accessibility for the NDSARAC. Finally, the results are applied to some physical examples in order to check the application potentiality.     

Analysis of Two Ultrafiltration Fouling Models and Estimation of Model Parameters as a Function of Operational Conditions

This work analyses the measure of fit of experimental data of permeate flux decline with time for ultrafiltration experiments performed with polyethylene glycol aqueous solutions to two different ultrafiltration models. A feed solution of 5 kg/m $^{3}$ 3 of polyethylene glycol and a monotubular ceramic membrane of ${\mathrm{ZrO}}_{2}$ ZrO 2 – ${\mathrm{TiO}}_{2}$ TiO 2 were used in the experiments. The first model considered was developed by Ho and Zydney and it considers two different fouling mechanisms: pore blocking and gel layer formation. The second model was proposed by Yee et al. It is an exponential model that considers three stages: concentration polarization, molecule deposition on the membrane surface and long-term fouling. The results show that both models give very accurate predictions for the severe fouling conditions (high transmembrane pressures and low crossflow velocities). However, both models cannot explain the experimental results obtained for all the experimental conditions tested. An equation for Ho and Zydney’s model parameters as a function of operating conditions was obtained by means of multiple regression analysis.     

A fast algorithm to calculate the critical coupling strength for synchronization in a chain of Kuramoto oscillators

Synchronization in a one-dimensional chain of Kuramoto oscillators with periodic boundary conditions is studied. An algorithm to rapidly calculate the critical coupling strength \(K_c\) for complete frequency synchronization is presented according to the mathematical constraint conditions and the periodic boundary conditions. By this new algorithm, we have checked the relation between \(\langle K_c\rangle \) and \(N\) , which is \(\langle K_c\rangle \sim \sqrt{N}\) , not only for small \(N\) , but also for large \(N\) . We also investigate the heavy-tailed distribution of \(K_c\) for random intrinsic frequencies, which is obtained by showing that the synchronization problem is equivalent to a discretization of Brownian motion. This theoretical result was checked by generating a large sample of \(K_c\) for large \(N\) from our algorithm to get the empirical density of \(K_c\) . Finally, we derive the permutation for the maximum coupling strength and its exact expression, which grows linearly with \(N\) and would provide the theoretical support for engineering applications.     

Semianalytical Solution for \text{ CO}_{2} Plume Shape and Pressure Evolution During \text{ CO}_{2} Injection in Deep Saline Formations

The injection of supercritical carbon dioxide ( $\text{ CO}_{2})$ in deep saline aquifers leads to the formation of a $\text{ CO}_{2}$ rich phase plume that tends to float over the resident brine. As pressure builds up, $\text{ CO}_{2}$ density will increase because of its high compressibility. Current analytical solutions do not account for $\text{ CO}_{2}$ compressibility and consider a volumetric injection rate that is uniformly distributed along the whole thickness of the aquifer, which is unrealistic. Furthermore, the slope of the $\text{ CO}_{2}$ pressure with respect to the logarithm of distance obtained from these solutions differs from that of numerical solutions. We develop a semianalytical solution for the $\text{ CO}_{2}$ plume geometry and fluid pressure evolution, accounting for $\text{ CO}_{2}$ compressibility and buoyancy effects in the injection well, so $\text{ CO}_{2}$ is not uniformly injected along the aquifer thickness. We formulate the problem in terms of a $\text{ CO}_{2}$ potential that facilitates solution in horizontal layers, with which we discretize the aquifer. Capillary pressure is considered at the interface between the $\text{ CO}_{2}$ rich phase and the aqueous phase. When a prescribed $\text{ CO}_{2}$ mass flow rate is injected, $\text{ CO}_{2}$ advances initially through the top portion of the aquifer. As $\text{ CO}_{2}$ is being injected, the $\text{ CO}_{2}$ plume advances not only laterally, but also vertically downwards. However, the $\text{ CO}_{2}$ plume does not necessarily occupy the whole thickness of the aquifer. We found that even in the cases in which the $\text{ CO}_{2}$ plume reaches the bottom of the aquifer, most of the injected $\text{ CO}_{2}$ enters the aquifer through the layers at the top. Both $\text{ CO}_{2}$ plume position and fluid pressure compare well with numerical simulations. This solution permits quick evaluations of the $\text{ CO}_{2}$ plume position and fluid pressure distribution when injecting supercritical $\text{ CO}_{2}$ in a deep saline aquifer.     

Drainage of Capillary-Trapped Oil by an Immiscible Gas: Impact of Transient and Steady-State Water Displacement on Three-Phase Oil Permeability

In a previous paper (Dehghanpour et al., Phys Rev E 83:065302, 2011a), we showed that relative permeability of mobilized oil, $k_\mathrm{ro}$ , measured during tertiary gravity drainage, is significantly higher than that of the same oil saturation in other tests where oil is initially a continuous phase. We also showed that tertiary $k_\mathrm{ro}$ strongly correlates to both water saturation, $S_\mathrm{w}$ , water flux (water relative permeability), $k_\mathrm{rw}$ , and the change in water saturation with time, $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ . To develop a model and understanding of the enhanced oil transport, identifying which of these parameters ( $S_\mathrm{w},\,k_{\mathrm{rw}}$ , or $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ ) plays the controlling role is necessary, but in the previous experiments these could not be deconvolved. To answer the remaining question, we conduct specific three-phase displacement experiments in which $k_{\mathrm{rw}}$ is controlled by applying a fixed water influx, and $S_\mathrm{w}$ develops naturally. We obtain $k_{\mathrm{ro}}$ by using the saturation data measured in time and space. The results suggest that steady-state water influx, in contrast to transient water displacement, does not enhance $k_{\mathrm{ro}}$ . Instead, reducing water influx rate results in excess oil flow. Furthermore, according to our pore scale hydraulic conductivity calculations, viscous coupling and fluid positioning do not sufficiently explain the observed correlation between $k_{\mathrm{ro}}$ and $S_{\mathrm{w}}$ . We conclude that tertiary $k_{\mathrm{ro}}$ is controlled by the oil mobilization rate, which in turn is linked to the rate of water saturation decrease with time, $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ . Finally, we develop a simple model which relates tertiary $k_{\mathrm{ro}}$ to transient two-phase gas/water relative permeability.