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.     

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.     

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.     

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.     

Robust mixed H_2 /H_\infty active control for offshore steel jacket platform

This paper presents a robust mixed \(H_2 /H_\infty \) control method for wave-excited offshore jacket platforms. Its objective was to design a controller that minimizes the upper bound of the \(H_2 \) performance measure on platform dynamics satisfying some \(H_\infty \) norm bound constraint simultaneously. Based on mixed \(H_2 /H_\infty \) control theory and linear matrix inequality techniques, a novel approach to stabilize offshore platform vibration with constrained \(H_2 /H_\infty \) performances is proposed. Uncertainties of the wave excitation are considered in dynamic performance analysis of offshore platforms. A reduced mode offshore platform structure under wave excitation is analyzed, and simulations are used to verify the effectiveness of the proposed approach. Compared with existing \(H_\infty \) control methods, the proposed approach makes a significant improvement for dynamic performances of offshore platforms under random wave excitation.     

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.     

Global stabilization of inherently non-linear systems using continuously differentiable controllers

This paper concerns the problem of constructing \(C^1\) (continuously differentiable) controllers to stabilize a class of uncertain non-linear systems whose linearization around the origin may contain uncontrollable modes. Based on a new definition of homogeneity with monotone degrees, a polynomial Lyapunov function and a \(C^1\) global stabilizer are constructed recursively. Moreover, several special cases are investigated to show the advantages of the proposed approaches using the generalized homogeneity compared to the existing approaches using the traditional homogeneity.     

On the Characterization of the Navier–Stokes Flows with the Power-Like Energy Decay

We study the energy decay of the turbulent solutions to the Navier–Stokes equations in the whole three-dimensional space. We show as the main result that the solutions with the energy decreasing at the rate \({O(t^{-\alpha}), t \rightarrow \infty, \alpha \in [0, 5/2]}\) , are exactly characterized by their initial conditions belonging into the homogeneous Besov space \({\dot{B}^{-\alpha}_{2, \infty}}\) . Similarly, for a solution u and \({p \in [1, \infty]}\) the integral \({\int_{0}^{\infty} \|t^{\alpha/2} u(t)\|^p \frac{1}{t} dt}\) is finite if and only if the initial condition of u belongs to the homogeneous Besov space \({\dot{B}_{2, p}^{-\alpha}}\) . For the case \({\alpha \in (5/2, 9/2]}\) we present analogical results for some subclasses of turbulent solutions.     

On the neutralization of bacterial spores in post-detonation flows

In multiple operational scenarios, explosive charges are used to neutralize confined or unconfined stores of bacterial spores. The spore destruction is achieved by post-detonation combustion and mixing of hot detonation product gases with the ambient flow and spore clouds. In this work, blast wave interaction with bacterial spore clouds and the effect of post-detonation combustion on spore neutralization are investigated using numerical simulations. Spherical explosive charges (radius, \(R_\mathrm{C}\) = 5.9 cm) comprising of nitromethane are modeled in the vicinity of a spore cloud, and the spore kill in the post-detonation flow is quantified. The effect of the mass of the spores and the initial distance, \(d^0\) , of the spore cloud from the explosive charge on the percentage of spores neutralized is investigated. When the spores are initially placed within a distance of 3.0 \(R_\mathrm{C}\) , within 0.1 ms after detonation of the charge, all the spores are neutralized by the blast wave and the hot detonation product gases. In contrast, almost all the spores survived the explosion when \(d^0\) is greater than 8.0 \(R_\mathrm{C}\) . The percentage of intact spores varied from 0 to 100 for 3.0 \(R_\mathrm{C}\) \( 8.0 \(R_\mathrm{C}\) with spore neutralization dependent on time spent by the spores in the post-detonation mixing/combustion zone.     

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.     

Robust {\varvec{L}}_{\varvec{2}}-gain control for a class of cascade switched nonlinear systems

In this paper, we study the robust finite \(L_2 \) -gain control for a class of cascade switched nonlinear systems with parameter uncertainty. Each subsystem of the switched system under consideration is composed of a zero-input asymptotically stable nonlinear part which is a lower dimension switched system, and of a linearizable part. The uncertainty appears in the control channel of each subsystem. We give sufficient conditions under which the nonlinear feedback controllers are derived to guarantee that the \(L_2 \) -gain of the closed-loop switched system is less than a prespecified value for all admissible uncertainty under arbitrary switching. Moreover, we also develop the \(L_2\) -gain controllers for the switched systems with nonminimum phase case.     

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\) .     

On Formation of a Locally Self-Similar Collapse in the Incompressible Euler Equations

The paper addresses the question of the existence of a locally self-similar blow-up for the incompressible Euler equations. Several exclusion results are proved based on the L ^p -condition for velocity or vorticity and for a range of scaling exponents. In particular, in N dimensions if in self-similar variables ${u \in L^p}$ and ${u \sim \frac{1}{t^{\alpha/(1+\alpha)}}}$ , then the blow-up does not occur, provided ${\alpha > N/2}$ or ${-1 < \alpha \leq N\,/p}$ . This includes the L ³ case natural for the Navier–Stokes equations. For ${\alpha = N\,/2}$ we exclude profiles with asymptotic power bounds of the form ${ |y|^{-N-1+\delta} \lesssim |u(y)| \lesssim |y|^{1-\delta}}$ . Solutions homogeneous near infinity are eliminated, as well, except when homogeneity is scaling invariant.     

Optimal Continuous Dependence Estimates for Fractional Degenerate Parabolic Equations

We derive continuous dependence estimates for weak entropy solutions of degenerate parabolic equations with nonlinear fractional diffusion. The diffusion term involves the fractional Laplace operator, \({\triangle^{\alpha/2}}\) for \({\alpha \in (0,2)}\) . Our results are quantitative and we exhibit an example for which they are optimal. We cover the dependence on the nonlinearities, and for the first time, the Lipschitz dependence on α in the BV-framework. The former estimate (dependence on nonlinearity) is robust in the sense that it is stable in the limits \({\alpha \downarrow 0}\) and \({\alpha \uparrow 2}\) . In the limit \({\alpha \uparrow 2}\) , \({\triangle^{\alpha/2}}\) converges to the usual Laplacian, and we show rigorously that we recover the optimal continuous dependence result of Cockburn and Gripenberg (J Differ Equ 151(2):231–251, 1999) for local degenerate parabolic equations (thus providing an alternative proof).     

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.     

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.     

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.     

Global Compensated Compactness Theorem for General Differential Operators of First Order

Let A ₁(x, D) and A ₂(x, D) be differential operators of the first order acting on l-vector functions ${u= (u_1, \ldots, u_l)}$ in a bounded domain ${\Omega \subset \mathbb{R}^{n}}$ with the smooth boundary ${\partial\Omega}$ . We assume that the H ¹-norm ${\|u\|_{H^{1}(\Omega)}}$ is equivalent to ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_1u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ and ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_2u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ , where B _i  = B _i (x, ν) is the trace operator onto ${\partial\Omega}$ associated with A _i (x, D) for i = 1, 2 which is determined by the Stokes integral formula (ν: unit outer normal to ${\partial\Omega}$ ). Furthermore, we impose on A ₁ and A ₂ a cancellation property such as ${A_1A_2^{\prime}=0}$ and ${A_2A_1^{\prime}=0}$ , where ${A^{\prime}_i}$ is the formal adjoint differential operator of A _i (i = 1, 2). Suppose that ${\{u_m\}_{m=1}^{\infty}}$ and ${\{v_m\}_{m=1}^{\infty}}$ converge to u and v weakly in ${L^2(\Omega)}$ , respectively. Assume also that ${\{A_{1}u_m\}_{m=1}^{\infty}}$ and ${\{A_{2}v_{m}\}_{m=1}^{\infty}}$ are bounded in ${L^{2}(\Omega)}$ . If either ${\{B_{1}u_m\}_{m=1}^{\infty}}$ or ${\{B_{2}v_m\}_{m=1}^{\infty}}$ is bounded in ${H^{\frac{1}{2}}(\partial\Omega)}$ , then it holds that ${\int_{\Omega}u_m\cdot v_m \,{\rm d}x \to \int_{\Omega}u\cdot v \,{\rm d}x}$ . We also discuss a corresponding result on compact Riemannian manifolds with boundary.