Regularizing Piecewise Smooth Differential Systems: Co-Dimension 2 Discontinuity Surface

In this paper, we are concerned with numerical solution of piecewise smooth initial value problems. Of specific interest is the case when the discontinuities occur on a smooth manifold of co-dimension 2, intersection of two co-dimension 1 singularity surfaces, and which is nodally attractive for nearby dynamics. In this case of a co-dimension 2 attracting sliding surface, we will give some results relative to two prototypical time and space regularizations. We will show that, unlike the case of co-dimension 1 discontinuity surface, in the case of co-dimension 2 discontinuity surface the behavior of the regularized problems is strikingly different. On the one hand, the time regularization approach will not select a unique sliding mode on the discontinuity surface, thus maintaining the general ambiguity of how to select a Filippov vector field in this case. On the other hand, the proposed space regularization approach is not ambiguous, and there will always be a unique solution associated to the regularized vector field, which will remain close to the original co-dimension 2 surface. We will further clarify the limiting behavior (as the regularization parameter goes to 0) of the proposed space regularization to the solution associated to the sliding vector field of Dieci and Lopez (Numer Math 117:779–811, 2011). Numerical examples will be given to illustrate the different cases and to provide some preliminary exploration in the case of co-dimension 3 discontinuity surface.     

The Sliding Bifurcations in Planar Piecewise Smooth Differential Systems

In this paper we are mainly interested in the bifurcation phenomena for a class of planar piecewise smooth differential systems, where a new phenomenon, i.e. sliding heteroclinic bifurcation, is found. Furthermore we will show that the involved systems can present many interesting bifurcation phenomena, such as the (sliding) heteroclinic bifurcation, sliding (homoclinic) cycle bifurcation and semistable limit cycle bifurcation and so on. The system can have two hyperbolic limit cycles, which are bifurcated in one way from a semistable limit cycle, and in another way from a heteroclinic cycle and a sliding cycle. In the proof of our main results, we will use the geometric singular perturbation theory to analyze the dynamics near the sliding region.     

Periodic Orbits in Hamiltonian Systems with Involutory Symmetries

We study the existence of families of periodic solutions in a neighbourhood of a symmetric equilibrium point in two classes of Hamiltonian systems with involutory symmetry. In both classes, the involution reverses the sign of the Hamiltonian function, and the system is in 1:1 resonance. In the first class we study a Hamiltonian system with a reversing involution R acting symplectically. We first recover a result of Buzzi and Lamb showing that the equilibrium point is contained in a three dimensional conical subspace which consists of a two parameter family of periodic solutions with symmetry R, and furthermore that there may or may not exist two families of non-symmetric periodic solutions, depending on the coefficients of the Hamiltonian (correcting a minor error in their paper). In the second problem we study an equivariant Hamiltonian system with a symmetry S that acts anti-symplectically. Generically, there is no S-symmetric solution in a neighbourhood of the equilibrium point. Moreover, we prove the existence of at least 2 and at most 12 families of non-symmetric periodic solutions. We conclude with a brief study of systems with both forms of symmetry, showing they have very similar structure to the system with symmetry R.     

Orbits Homoclinic to Exponentially Small Periodic Orbits for a Class of Reversible Systems. Application to Water Waves

In this paper a class of reversible analytic vector fields is investigated near an equilibrium. For these vector fields, the part of the spectrum of the differential at the equilibrium which lies near the imaginary axis comes from the perturbation of a double eigenvalue 0 and two simple eigenvalues , . In the first part of this paper, we study the 4-dimensional problem. The existence of a family of solutions homoclinic to periodic orbits of size less than μ^N for any fixed N, where μ is the bifurcation parameter, is known for vector fields. Using the analyticity of the vector field, we prove here the existence of solutions homoclinic to a periodic orbit the size of which is exponentially small ( of order . This result receives its significance from the still unsolved question of whether there exist solutions that are homoclinic to the equilibrium or whether the amplitudes of the oscillations at infinity have a positive infimum. In the second part of this paper we prove that the exponential estimates still hold in infinite dimensions. This result cannot be simply obtained from the study of the 4-dimensional analysis by a center-manifold reduction since this result is based on analyticity of the vector field. One example of such a vector field in infinite dimensions occurs when describing the irrotational flow of an inviscid fluid layer under the influence of gravity and small surface tension (Bond number ) for a Froude number F close to 1. In this context a homoclinic solution to a periodic orbit is called a generalized solitary wave. Our work shows that there exist generalized solitary waves with exponentially small oscillations at infinity. More precisely, we prove that for each F close enough to 1, there exist two reversible solutions homoclinic to a periodic orbit, the size of which is less than , l being any number between 0 and π and satisfying . (Accepted October 2, 1995)     

Planar Vector Fields with a Given Set of Orbits

We determine all the \({\mathcal{C}^1}\) planar vector fields with a given set of orbits of the form y ? y(x) = 0 satisfying convenient assumptions. The case when these orbits are branches of an algebraic curve is also study. We show that if a quadratic vector field admits a unique irreducible invariant algebraic curve \({g(x, y) = \sum_{j=0}^S a_j(x) y^{S-j}= 0}\) with S branches with respect to the variable y, then the degree of the polynomial g is at most 4S.     

Numerical Continuation of Homoclinic Orbits to Non-Hyperbolic Equilibria in Planar Systems

In this paper we develop numerical algorithms for thecontinuation of degenerate homoclinic orbits to non-hyperbolicequilibria in planar systems. The first situation corresponds to asaddle-node equilibrium (a zero eigenvalue) and the second one is theso-called cuspidal loop (double-zero eigenvalue). The methods proposedmay deal with codimension-two and -three homoclinic connections.Application of the algorithms to several examples supports its validityand demonstrates its usefulness.     

Averaging Theory of Arbitrary Order for Piecewise Smooth Differential Systems and Its Application

The averaging theory for studying periodic orbits of smooth differential systems has a long history. Whereas the averaging theory for piecewise smooth differential systems appeared only in recent years, where the unperturbed systems are smooth. When the unperturbed systems are only piecewise smooth, there is not an existing averaging theory to study existence of periodic orbits of their perturbed systems. Here we establish such a theory for one dimensional perturbed piecewise smooth periodic differential equations. Then we show how to transform planar perturbed piecewise smooth differential systems to one dimensional piecewise smooth periodic differential equations when the unperturbed planar piecewise smooth differential systems have a family of periodic orbits. Finally as application of our theory we study limit cycle bifurcation of planar piecewise differential systems which are perturbation of a \(\Sigma \)-center.     

The Full Study of Planar Quadratic Differential Systems Possessing a Line of Singularities at Infinity

In this article we make a full study of the class of non-degenerate real planar quadratic differential systems having all points at infinity (in the Poincaré compactification) as singularities. We prove that all such systems have invariant affine lines of total multiplicity 3, give all their configurations of invariant lines and show that all these systems are integrable via the method of Darboux having cubic polynomials as inverse integrating factors. After constructing the topologically distinct phase portraits in this class we give invariant necessary and sufficient conditions in terms of the 12 coefficients of the systems for the realization of each one of them and give representatives of the orbits under the action of the affine group and time rescaling. We construct the moduli space of this class for this action and give the corresponding bifurcation diagram. Dedicated to Professor Zhifen Zhang on the occasion of her 80th birthday     

Continuation of the Periodic Orbits for the Differential Equation with Discontinuous Right Hand Side

In this paper the discontinuous system with one parameter perturbation is considered. Here the unperturbed system is supposed to have at least either one periodic orbit or a limit cycle. The aim is to prove the continuation of the periodic orbits under perturbation by means of the bifurcation map and the zeroes of this map imply the periodic orbits for the perturbed system. The tools for this problem are jumps of fundamental matrix solutions and the Poincare map for discontinuous systems. Therefore, we develop the Diliberto theorem and variation lemma for the system with discontinuous right hand side. At the end, as application of our method, the effect of discontinuous damping on Van der pol equation, and the effect of small force on the discontinuous linear oscillator with add a ·sgn(x) are considered.     

Isoenergetic Families of Planar Orbits Generated by Homogeneous Potentials

In the light of the inverse problem of dynamics, we study in some detail the monoparametric isoenergetic families of planar orbits f (x, y) = c, created by homogeneous potentials V (x, y). For any preassigned family of orbits and for any degree of homogeneity m of the potential, we offer the criteria which the family has to satisfy so that it can result from such a potential. When the criteria are fulfilled, Szebehelys first order partial differential equation for the unknown potential V (x, y) is substantially simplified and, in most of the cases, can be solved to completion and uniqueness.     

Unique Periodic Orbits for Delayed Positive Feedback and the Global Attractor

The delay differential equation, (t)=–x(t)+f(x(t–1)), with >0 and a real function f satisfying f(0)=0 and f>0 models a system governed by delayed positive feedback and instantaneous damping. Recently the geometric, topological, and dynamical properties of a three-dimensional compact invariant set were described in the phase space C=C([–1, 0], ) of initial data for solutions of the equation. In this paper, for a set of and f which include examples from neural network theory, we show that this three-dimensional set is the global attractor, i.e., the compact invariant set which attracts all bounded subsets of C. The proof involves, among others, results on uniqueness and absence of periodic orbits.     

Parabolic Systems with Nowhere Smooth Solutions

We construct smooth 2×2 parabolic systems with smooth initial data and C^α right-hand side which admit solutions that are nowhere C¹. The elliptic part is in variational form and the corresponding energy ϕ is strongly quasiconvex and in particular satisfies a uniform Legendre-Hadamard (or strong ellipticity) condition.     

Global Existence of Smooth Solutions for Partially Dissipative Hyperbolic Systems with a Convex Entropy

We consider the Cauchy problem for a general one-dimensional n×n hyperbolic symmetrizable system of balance laws. It is well known that, in many physical examples, for instance for the isentropic Euler system with damping, the dissipation due to the source term may prevent the shock formation, at least for smooth and small initial data. Our main goal is to find a set of general and realistic sufficient conditions to guarantee the global existence of smooth solutions, and possibly to investigate their asymptotic behavior. For systems which are entropy dissipative, a quite natural generalization of the Kawashima condition for hyperbolic-parabolic systems can be given. In this paper, we first propose a general framework for this kind of problem, by using the so-called entropy variables. Then we go on to prove some general statements about the global existence of smooth solutions, under different sets of conditions. In particular, the present approach is suitable for dealing with most of the physical examples of systems with a relaxation extension. Our main tools will be some refined energy estimates and the use of a suitable version of the Kawashima condition.     

An Analysis of Journal Orbits for Nonlinear Dynamic Bearing Systems

An analysis of journal centre orbits is presented in this paper based on a non-isothermal non-Newtonian fluid model for dynamically loaded bearing systems. A spectral element approach is used to solve a full set of coupled equations (kinematics and constitutive) governing the flow of the lubricant, and an operator-splitting spectral element technique is used to evaluate the dynamic energy equation. The motion of the journal is calculated on the basis of Newtonian mechanics incorporated with a simple cavitation model. The stability of the journal orbits is investigated under a wide range of the rotation speeds of journal. The unstable orbits arise as a sub-harmonic motion when the journal rotation speed is increased beyond a critical value. The influences of the oscillation speeds of the applied loads on the journal orbits are examined. The numerical simulations demonstrate that both the rotation speed of the journal and the oscillation speed of the applied load play an important role in determining the pattern of the journal orbits. The effects of square-wave and rotating applied loads on the journal orbits are also investigated. Received 22 April 1998 and accepted 26 May 1999     

Homogenization of Periodic Systems with Large Potentials

We consider the homogenization of a system of second-order equations with a large potential in a periodic medium. Denoting by the period, the potential is scaled as ^–2. Under a generic assumption on the spectral properties of the associated cell problem, we prove that the solution can be approximately factorized as the product of a fast oscillating cell eigenfunction and of a slowly varying solution of a scalar second-order equation. This result applies to various types of equations such as parabolic, hyperbolic or eigenvalue problems, as well as fourth-order plate equation. We also prove that, for well-prepared initial data concentrating at the bottom of a Bloch band, the resulting homogenized tensor depends on the chosen Bloch band. Our method is based on a combination of classical homogenization techniques (two-scale convergence and suitable oscillating test functions) and of Bloch waves decomposition.     

Binary Decompositions for Planar N-Body Problems and Symmetric Periodic Solutions

A binary decomposition for a system of N masses is a way of treating the system as binaries with the total action exactly the same as that of the original system. By considering binary decompositions, we are able to provide effective lower-bound estimates for the action of collision paths in several spaces of symmetric loops. As applications, we use our estimates to prove the existence of some new classes of symmetric periodic solutions for the N-body problem.     

Representation of a Smooth Isometric Deformation of a Planar Material Region into a Curved Surface

We consider the problem of characterizing the smooth, isometric deformations of a planar material region identified with an open, connected subset \({\mathcal{D}}\) of two-dimensional Euclidean point space \(\mathbb{E}^{2}\) into a surface \({\mathcal{S}}\) in three-dimensional Euclidean point space \(\mathbb{E}^{3}\). To be isometric, such a deformation must preserve the length of every possible arc of material points on \({\mathcal{D}}\). Characterizing the curves of zero principal curvature of \({\mathcal{S}}\) is of major importance. After establishing this characterization, we introduce a special curvilinear coordinate system in \(\mathbb{E}^{2}\), based upon an à priori chosen pre-image form of the curves of zero principal curvature in \({\mathcal{D}}\), and use that coordinate system to construct the most general isometric deformation of \({\mathcal{D}}\) to a smooth surface \({\mathcal{S}}\). A necessary and sufficient condition for the deformation to be isometric is noted and alternative representations are given. Expressions for the curvature tensor and potentially nonvanishing principal curvature of \({\mathcal{S}}\) are derived. A general cylindrical deformation is developed and two examples of circular cylindrical and spiral cylindrical form are constructed. A strategy for determining any smooth isometric deformation is outlined and that strategy is employed to determine the general isometric deformation of a rectangular material strip to a ribbon on a conical surface. Finally, it is shown that the representation established here is equivalent to an alternative previously established by Chen, Fosdick and Fried (J. Elast. 119:335–350, 2015).     

Limit Cycles of Piecewise Smooth Differential Equations on Two Dimensional Torus

In this paper we study the limit cycles of some classes of piecewise smooth vector fields defined in the two dimensional torus. The piecewise smooth vector fields that we consider are composed by linear, Ricatti with constant coefficients and perturbations of these one, which are given in (3). Considering these piecewise smooth vector fields we characterize the global dynamics, studying the upper bound of number of limit cycles, the existence of non-trivial recurrence and a continuum of periodic orbits. We also present a family of piecewise smooth vector fields that posses a finite number of fold points and, for this family we prove that for any 2k number of limit cycles there exists a piecewise smooth vector fields in this family that presents k number of limit cycles and prove that some classes of piecewise smooth vector fields presents a non-trivial recurrence or a continuum of periodic orbits.