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.     

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.     

Random parameters induce chaos in power systems

We study how random parameter (namely, noise-perturbed parameter) effects the dynamical behaviors of power systems using random Melnikov technique and numerical simulation. The studied model is described by the classical single-machine-infinite-bus systems which operate in a stable periodic regime far away from chaotic behavior with deterministic parameter. It is found that when the parameter perturbations are weak, chaos is absent in power systems. With the intensity of random parameter \(\rho \) increasing, power systems become unstable and fall into chaos as \(\rho \) is further increased. These phenomena imply that random parameter can induce and enhance chaos in power systems. Our results may provide a useful tip for understanding power systems’ security operation.     

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.     

Expansivity and Cone-fields in Metric Spaces

Due to the results of Lewowicz and Tolosa expansivity can be characterized with the aid of Lyapunov function. In this paper we study a similar problem for uniform expansivity and show that it can be described using generalized cone-fields on metric spaces. We say that a function \(f:X\rightarrow X\) is uniformly expansive on a set \(\varLambda \subset X\) if there exist \(\varepsilon >0\) and \(\alpha \in (0,1)\) such that for any two orbits \(\hbox {x}:\{-N,\ldots ,N\} \rightarrow \varLambda \) , \(\hbox {v}:\{-N,\ldots ,N\} \rightarrow X\) of \(f\) we have $$\begin{aligned} \sup _{-N\le n\le N}d(\hbox {x}_n,\hbox {v}_n) \le \varepsilon \implies d(\hbox {x}_0,\hbox {v}_0) \le \alpha \sup _{-N\le n\le N}d(\hbox {x}_n,\hbox {v}_n). \end{aligned}$$ It occurs that a function is uniformly expansive iff there exists a generalized cone-field on \(X\) such that \(f\) is cone-hyperbolic.     

Gas Sorption and the Consequent Volumetric and Permeability Change of Coal II: Numerical Modeling

The permeability of coalbed methane reservoirs may evolve during the recovery of methane and injection of gas, due to the change of effective stress and gas adsorption and desorption. Experimental and numerical studies were conducted to investigate the sorption-induced permeability change of coal. This paper presents the numerical modeling part of the work. It was found that adsorption of pure gases on coal was well represented by parametric adsorption isotherm models in the literature. Based on the experimental data of this study, adsorption of pure \(\hbox {N}_2\) was modeled using the Langmuir equation, and adsorption of pure \(\hbox {CO}_2\) was well represented by the N-Layer BET equation. For the modeling of CO \(_2\) & N \(_2\) binary mixture adsorption, the ideal adsorbed solution (IAS) model and the real adsorbed solution (RAS) model were used. The IAS model estimated the total amount of mixture adsorption and the composition of the adsorbed phase based on the pure adsorption isotherms. The estimated total adsorption and adsorbed-phase composition were very different from the experimental results, indicating nonideality of the CO \(_2\) –N \(_2\) –Coal-adsorption system. The measured sorption-induced strain was linearly proportional to the total amount of adsorption despite the species of the adsorbed gas. Permeability reduction followed a linear correlation with the volumetric strain with the adsorption of pure \(\hbox {N}_2\) and the tested CO \(_2\) & N \(_2\) binary mixtures, and an exponential correlation with the adsorption of pure \(\hbox {CO}_2\) .     

Investigation of Self-Similar Interface Evolution in Carbon Dioxide Sequestration in Saline Aquifers

In this work, we investigate numerically the injection of supercritical carbon dioxide into a deep saline reservoir from a single well. We analyze systematically the sharp-interface evolution in different flow regimes. The flow regimes can be parameterized by two dimensionless numbers, the gravity number, \(\Gamma \) and the mobility ratio, \(\lambda \) . Numerical simulations are performed using the volume of fluid method, and the results are compared with the solutions of the self-similarity equation established in previous works, which describes the evolution of the sharp interface. We show that these theoretical solutions are in very good agreement with the results from the numerical simulations presented over the different flow regimes, thereby showing that the theoretical and simulation models predict consistently the spreading and migration of the created \(\hbox {CO}_{2}\) plume under complex flow behavior in porous media. Furthermore, we compare the numerical results with known analytic approximations in order to assess their applicability and accuracy over the investigated parametric space. The present study indicates that the self-similar solutions parameterized by the dimensionless numbers \(\lambda , \Gamma \) are significant for examining effectively injection scenarios, as these numbers control the shape of the interface and migration of the \(\hbox {CO}_{2}\) plume. This finding is essential in assessing the storage capacity of saline aquifers.     

Output feedback \varvec{\mathcal {H}}_{{\varvec{\infty }}} control of stochastic nonlinear time-delay systems with state and disturbance-dependent noises

This paper is concerned with the output feedback \(\mathcal {H}_\infty \) control problem for a class of stochastic nonlinear systems with time-varying state delays; the system dynamics is governed by the stochastic time-delay It \(\hat{o}\) -type differential equation with state and disturbance contaminated by white noises. The design of the output feedback \(\mathcal {H}_\infty \) control is based on the stochastic dissipative theory. By establishing the stochastic dissipation of the closed-loop system, the delay-dependent and delay-independent approaches are proposed for designing the output feedback \(\mathcal {H}_\infty \) controller. It is shown that the output feedback \(\mathcal {H}_\infty \) control problem for the stochastic nonlinear time-delay systems can be solved by two delay-involved Hamilton–Jacobi inequalities. A numerical example is provided to illustrate the effectiveness of the proposed methods.     

Elucidating the Role of Interfacial Tension for Hydrological Properties of Two-Phase Flow in Natural Sandstone by an Improved Lattice Boltzmann Method

We investigated the interfacial tension (IFT) effect on fluid flow characteristics inside micro-scale, porous media by a highly efficient multi-phase lattice Boltzmann method using a graphics processing unit. IFT is one of the most important parameters for carbon capture and storage and enhanced oil recovery. Rock pores of Berea sandstone were reconstructed from micro-CT scanned images, and multi-phase flows were simulated for the digital rock model at extremely high resolution (3.2  \(\upmu \) m). Under different IFT conditions, numerical analyses were carried out first to investigate the variation in relative permeability, and then to clarify evolution of the saturation distribution of injected fluid. We confirmed that the relative permeability decreases with increasing IFT due to growing capillary trapping intensity. It was also observed that with certain pressure gradient \(\Delta P\) two crucial IFT values, \(\sigma _{1}\) and \(\sigma _{2}\) , exist, creating three zones in which the displacement process has totally different characteristics. When \(\sigma _{1}< \sigma < \sigma _{2}\) , the capillary fingering patterns are observed, while for \(\sigma < \sigma _{1}\) viscous fingering is dominant and most of the passable pore spaces were invaded. When \(\sigma > \sigma _{2}\) the invading fluid failed to break through. The pore-throat-size distribution estimated from these crucial IFT values ( \(\sigma _{1 }\) and \(\sigma _{2})\) agrees with that derived from mercury porosimetry measurements of Berea sandstone. This study demonstrates that the proposed numerical method is an efficient tool for investigating hydrological properties from pore structures.     

Intermittent control on switched networks via {\varvec{\omega }}-matrix measure method

This paper is concerned with global exponential synchronization problem for a class of switched delay networks with interval parameters uncertainty, different from the most existing results, without constructing complex Lyapunov–Krasovskii functions; \(\omega \) -matrix measure method is firstly introduced to switched interval networks, combining Halanay inequality technique, designing proper intermittent and non-intermittent control strategy; some easy-to-verify synchronization criteria are given to ensure the global exponential synchronization of switched interval networks under arbitrary switching rule and for admissible interval uncertainties. Moreover, as an application, the proposed scheme can be applied to chaotic neural networks. Finally, numerical simulations are provided to illustrate the effectiveness of the theoretical results and show the obtained results via employing \(\omega \) -measure are superior to previous results by using \(1\) -measure.     

First-principles study of full a-dislocations in pure magnesium

Full a-dislocations on the (0001) basal plane, \((10\bar 10)\) prismatic plane, and \((10\bar 11)\) and \((10\bar 12)\) pyramidal planes in pure magnesium are investigated by using the Peierls-Nabarro model combined with generalized stacking fault (GSF) energies from first-principles calculations. The results show that the \(\left( {10\bar 11} \right)\left\langle {11\bar 20} \right\rangle\) and \(\left( {10\bar 12} \right)\left\langle {11\bar 20} \right\rangle\) slip modes have nearly the same GSF energy barriers, which are obviously larger than the GSF energy barriers of the \(\left( {0001} \right)\left\langle {11\bar 20} \right\rangle\) and \(\left( {10\bar 10} \right)\left\langle {11\bar 20} \right\rangle\) slip modes. For both edge and screw full dislocations, the maximum dislocation densities, Peierls energies, and stresses of dislocations on the \((10\bar 10)\) , (0001), \((10\bar 11)\) , and \((10\bar 12)\) planes eventually increase. Moreover, the Peierls energies and the stresses of screw full dislocations are always lower than those of edge full dislocations for all slip systems. Dislocations on the \((10\bar 11)\) and \((10\bar 12)\) pyramidal planes possess smaller core energies, while the \((10\bar 10)\) prismatic plane has the largest ones, implying that the formation of full dislocations on the \((10\bar 10)\) plane is more difficult.     

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.     

Fractional models of anomalous relaxation based on the Kilbas and Saigo function

We revisit the Kilbas and Saigo functions of the Mittag-Leffler type of a real variable \(t\) , with two independent real order-parameters. These functions, subjected to the requirement to be completely monotone for \(t>0\) , can provide suitable models for the responses and for the corresponding spectral distributions in anomalous (non–Debye) relaxation processes, found e.g. in dielectrics. Our analysis includes as particular cases the classical models referred to as Cole–Cole (the one-parameter Mittag-Leffler function) and to as Kohlrausch (the stretched exponential function). After some remarks on the Kilbas and Saigo functions, we discuss a class of fractional differential equations of order \(\alpha \in (0,1]\) with a characteristic coefficient varying in time according to a power law of exponent \(\beta \) , whose solutions will be presented in terms of these functions. We show 2D plots of the solutions and, for a few of them, the corresponding spectral distributions, keeping fixed one of the two order-parameters. The numerical results confirm the complete monotonicity of the solutions via the non-negativity of the spectral distributions, provided that the parameters satisfy the additional condition \(0<\alpha +\beta \le 1\) , assumed by us.     

On solutions of two coupled fractional time derivative Hirota equations

We consider the well-known nonlinear Hirota equation (NLH) with fractional time derivative and derive its periodic wave solution and approximate analytic solitary wave solution using the homotopy analysis method (HAM). We also apply HAM to two coupled time fractional NLHs and construct their periodic wave solution and approximate solitary wave solution. We observe that the obtained periodic wave solution in both cases can be written in terms of the Mittag–Leffler function when the convergence control parameter \({c}_0\) equals \(-1\) . Convergence of the obtained solution is discussed. The derived approximate analytic solution and the effect of time-fractional order \(\alpha \) are shown graphically.     

Miscible and Immiscible Foam Injection for Mobility Control and EOR in Fractured Oil-Wet Carbonate Rocks

Foam injection is a proven enhanced oil recovery (EOR) technique for heterogeneous reservoirs, but is less studied for EOR in fractured systems. We experimentally investigated tertiary \(\text {CO}_{2}\) injections, and \(\text {N}_{2}\) - and \(\text {CO}_{2}\) -foam injections for enhanced oil recovery in fractured, oil-wet limestone core plugs. Miscible \(\text {CO}_{2}\) and \(\text {CO}_{2}\) -foam was compared with immiscible \(\text {CO}_{2}\) - and \(\text {N}_{2}\) -foam as tertiary recovery techniques, subsequent to waterfloods, in fractured rocks with different wettability preferences. At water-wet conditions waterfloods produced approximately 40 % OOIP, by spontaneous imbibition. Waterflood oil recovery at oil-wet conditions was below 20 % OOIP, due to suppressed imbibition where water predominantly flowed through the fractures, unable to mobilize the oil trapped in the matrix. Tertiary, supercritical \(\text {CO}_{2}\) -mobilized oil trapped in the matrix, particularly at weakly oil-wet conditions, by diffusion. Recovery by diffusion was high due to small core samples, high initial oil saturation and a continuous oil phase at oil-wet conditions. Both immiscible \(\text {CO}_{2}\) - and \(\text {N}_{2}\) -foams and miscible, supercritical \(\text {CO}_{2}\) -foam demonstrated high ultimate oil recoveries, but immiscible foam was less efficient (30 pore volumes injected) compared to miscible foam (2 pore volumes injected) to reach ultimate recovery. This is explained by the capillary threshold pressure preventing the injected \(\text {N}_{2}\) gas from entering the matrix, verified by computed X-ray tomography, and the mobilized oil was displaced by the aqueous surfactant in the foam. At miscible conditions, there exists no capillary entry pressure between the oil-saturated matrix and the injected \(\text {CO}_{2}\) , allowing foam to invade the matrix for efficient oil recovery.     

Tomographic PIV investigation of roughness-induced transition in a hypersonic boundary layer

The disagreement between free surface scalar experiments and the two-dimensional (2D) transport equation is discussed. An effective diffusivity coefficient, \(\kappa _{{\rm eff}}\) , is introduced and defined as the quotient between variance decay and mean gradient square. In all the experiments performed, \(\kappa _{{\rm eff}}\) is significantly larger than the scalar diffusivity, \(\kappa \) . Three mechanisms are identified as responsible for the differences between the quasi two-dimensional (Q2D) experiments and the 2D behaviour of a diffusive scalar. These are the vertical velocity gradients, the free surface divergence and the gravity currents induced by the scalar. These mechanisms, which affect the diffusive term in the 2D transport equation for large Péclet number ( \(Pe\gg 1\) ), are evaluated for steady and time-dependant laminar flows driven by electromagnetic body forces.     

Estimability Analysis for Optimization of Hysteretic Soil Hydraulic Parameters Using Data of a Field Irrigation Experiment

To improve the quality of parameter optimization, estimability analysis has been proposed as the first step before inverse modeling. When using field data of irrigation experiments for the determination of soil hydraulic parameters, wetting and drying processes may complicate optimization. The objectives of this study were to compare estimability analysis and inverse optimization of the soil hydraulic parameters in the models with and without considering hysteresis of the soil water retention function. Soil water pressure head data of a field irrigation experiment were used. The one-dimensional vertical water movement in variably saturated soil was described with the Richards equation using the HYDRUS-1D code. Estimability of the unimodal van Genuchten–Mualem hydraulic model parameters as well as of the hysteretic parameter model of Parker and Lenhard was classified according to a sensitivity coefficient matrix. The matrix was obtained by sequentially calculating effects of initial parameter variations on changes in the simulated pressure head values. Optimization was carried out by means of the Levenberg-Marquardt method implemented in the HYDRUS-1D code. The parameters \(\alpha , K_{s}, \theta _{s}\) , and \(n\) in the nonhysteretic model were found sensitive and parameter \(\theta _{s}\) strongly correlated with parameter \(n\) . When assuming hysteresis, the estimability was decreased with soil depth for \(K_{s}\) and \(\alpha ^{d}\) , and increased for \(\theta _{s}\) and n. Among the shape parameters, \(\alpha ^{w}\) was the most estimable. The hysteretic model could approximate the pressure heads in the soil by considering parameters from wetting and drying periods separately as initial estimates. The inverse optimization could be carried out more efficiently with most estimable parameters. Despite the remaining weaknesses of the local optimization algorithm and the inflexibility of the unimodal van Genuchten model, the results suggested that estimability analysis could be considered as a guidance to better define the optimization scenarios and then improved the determination of soil hydraulic parameters.     

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.     

Flow Dynamics of \hbox {CO}_{2}/brine at the Interface Between the Storage Formation and Sealing Units in a Multi-layered Model

Pressure distribution and \(\hbox {CO}_{2}\) plume migration are two major interests in \(\hbox {CO}_{2}\) geologic storage as they determine the injectivity and storage capacity. In this study, we adopted a three-layer model comprising a storage formation and the over- and underlying seals and determined three distinct flow regions based on the vertical flux exchange of \(\hbox {CO}_{2}\) and native brine. Regions 1 and 2 showed \(\hbox {CO}_{2}\) flowing from the storage formation to adjacent seals with counter-flowing brine. The characteristics of these fluxes in Region 1 were governed by permeability change due to salt precipitation whereas buoyancy force controlled the flux pattern in Region 2. Region 3 showed brine flowing from storage formation toward the over- and underlying seals, which enabled the displaced brine to escape from the storage formation and make room for \(\hbox {CO}_{2}\) to store as well as reduce the pressure build-up. In the multi-layered model, the counter-flowing brine in flow Region 1 resulted in localized salt precipitation at the upper and lower boundary of storage formation. We assessed the bottom-hole pressure and \(\hbox {CO}_{2}\) mass in caprock with respect to reservoir size. While the formation thickness influenced the bottom-hole pressure in the early stage of injection, the horizontal extension of the reservoir was more influential to pressure build-up during the injection period, and to the stabilized pressure during the post-injection period. The \(\hbox {CO}_{2}\) mass in caprock gently increased during the injection period as well as during the post-injection period and reached about 4–5 % of injected \(\hbox {CO}_{2}\) . The percentage of escaped brine from the storage formation ranged from 80–100 % of the \(\hbox {CO}_{2}\) mass stored in the storage formation depending on the reservoir scale.     

Integrability with symbolic computation on the Bogoyavlensky–Konoplechenko model: Bell-polynomial manipulation,bilinear representation,and Wronskian solution

With symbolic computation, this paper investigates some integrable properties of a two-dimensional generalization of the Korteweg-de Vries equation, i.e., the Bogoyavlensky–Konoplechenko model, which can govern the interaction of a Riemann wave propagating along the \(y\) -axis and a long wave propagating along the \(x\) -axis. Within the framework of Bell-polynomial manipulations, Bell-polynomial expressions are firstly given, which then are cast into bilinear forms. The \(N\) -soliton solutions in the form of an \(N\) th-order polynomial in the \(N\) exponentials and in terms of the Wronskian determinant are, respectively, constructed with the Hirota bilinear method and Wronskian technique. Bilinear Bäcklund transformation is also derived with the achievement of a family of explicit solutions.