1. Introduction
There are hundreds of papers, books, theses and lecture notes that deal with the analysis of linear oscillatory systems, which are affected by a small nonlinear perturbation. The pioneering works of Poincaré [1] and Lindstedt [2] , led to the development of the Methods of Normal Forms [3] - [10] , Averaging [3] [11] [12] [13] and Multiple Time Scales [14] [15] [16] [17] [18] . Typically nonlinear oscillatory systems are the first to be studied in these references. A variety of new perturbative procedures that apply to problems, in which the nonlinearity is not small relative to the linear term, have emerged over the years. Some of these are reviewed in [19] [20] [21] [22] [23] . The asymptotic behavior of conservative oscillatory systems with polynomial nonlinearities as the highest power in the polynomial tends to infinity has been studied in the literature using a variety of expansion methods [24] - [35] . Tools for the analysis of integrals of the solutions of oscillatory systems with polynomial nonlinearity, the equation of which is
(1)
have been developed in terms of the inverse of the incomplete Beta function (Ateb(h) functions) [24] [34] [35] . Of particular relevance to the present paper is the expansion in a series of saw-tooth functions to describe the dynamics of both vibro-impact and other strongly nonlinear systems [26] [27] .
In the procedures exploited in these references the amplitude, x(t), and the period, T, are approximated by a power expansion in a small parameter through a finite order. The approximation for T leads to the emergence of secular terms, which spoil the quality of the approximation to the solution if computed for many cycles of oscillation. However, secular terms can be avoided in the analysis of energy-conserving oscillatory systems in (1 + 1) and some (1 + 2)-dimensional systems, as the period of oscillations, T, is then either known in closed form (e.g., central-field motion in a plane [36] , the Duffing [37] [38] [39] [40] , quintic [41] [42] [43] [44] and relativistic harmonic [45] [46] [47] oscillators), or can be computed to any desired level of accuracy. This idea has been exploited in the analysis of the case of a polynomial nonlinearity when both the maximal oscillation amplitude, a(E), and the highest power in the polynomial, (2N + 1), both tend to infinity [48] .
This study focuses on the behavior of nonlinear oscillatory systems in the limit of extremely high oscillation energy, E (mathematically, E → ∞), when the nonlinear driving force grows with E much faster for x(t) close to the turning point, a(E), than at any position, x(t), that is not too close to a(E). In many such systems, this behavior of the driving force determines important aspects of the solution in the limit E → ∞. In some systems it even determines the full profile of the solution in that limit. Hence, such dominance of the nonlinear part of the driving force will be called “nonlinear violence”. The main message is that when nonlinear violence dominates, the solution develops boundary-layer characteristics, which emerge independently of whether the equation of motion of the oscillator can or cannot be transformed into the traditional form of a boundary-layer system—a boundary-value problem described by an ODE, in which the highest derivative is multiplied by a small parameter, e.g.:
(2)
A boundary-layer analysis has been applied in hundreds of papers to the study of oscillation-relaxation systems. For the study of such problems, see, e.g., [25] [49] [50] . What characterizes the solutions of boundary-layer problems is that an approximation to the solution can be found by solving Equation (2) over two parts of the range of x. One first solves Equation (2) over the range of values of x, for which the first term in the equation,
, is negligible owing to the smallness of ε, and
being bounded. The remainder of Equation (2) then generates an approximate solution that varies slowly. That part of the solution is called the “outer solution”; it covers most of the interval of time throughout the period of oscillations. One then solves for Equation (2) in the narrow range of t, called the “boundary layer” where
is large owing to the rapid change in x(t). The width of this range diminishes to zero as ε → 0 and
increases indefinitely as ε → 0. That part of the solution is called the “inner solution”. Both “inner” and “outer” solutions are determined up to some free parameters. The approximate solution for the system over the whole range of x is found by a “matching” procedure: requiring that the values of the “inner” and “outer” solutions and of their first derivatives coincide at some matching point, tm. This determines the free parameters, yields the value of the matching point, tm and the narrow range in x in the vicinity f tm, over which x(t) varies rapidly. As ε → 0, the variation of x(t) throughout the boundary layer becomes faster and the width of the layer tends to zero.
One purpose of this paper is to show that solutions of nonlinear oscillatory systems that are dominated by nonlinear violence exhibit characteristics that are similar to those of boundary-layer problems. When nonlinear violence occurs in the vicinity of a turning point, the solution varies rapidly around that point over a short time interval, the length of which tends to zero as E → ∞. In boundary-layer terminology, this would be called the “inner solution”. Beyond this interval, the solution rapidly evolves into a moderate profile, often, well approximated by linear time dependence, in other systems—by a constant. In boundary-layer terminology, this would be called the “outer solution”. Exploiting the analysis tools developed in boundary-layer theory facilitates understanding of the characteristics of the solutions, in particular, in the limit of very high oscillation energy.
Section 2 deals with (1 + 1)-dimensional conservative systems in the limit of total energy, E → ∞. Systems not transformable to boundary-layer formulation are analyzed in Section 2.1. The first to be discussed, are systems, the solutions of which are symmetric under space inversion and tend to saw-function patterns. The general case of potential energy, V(x), that is symmetric under reflection in x, when V(x) grows sufficiently rapidly near the turning points, x = ±a(E) is studied in Section 2.1.1. An example of a Gaussian potential is presented in Section 2.1.2. The case of a polynomial driving force with constant coefficients with a positive coefficient multiplying the leading power is reviewed in Section 2.1.3. The case of the extreme relativistic limit of the relativistic harmonic oscillator is discussed in Section 2.2. Unlike all previously discussed systems, in the extreme relativistic limit, the dynamical equation of the relativistic harmonic oscillator is transformable to a boundary-layer problem. In all the system discussed, the solutions are governed by nonlinear violence. The growth characteristics of V(a(E)) lead to a profile of x(t), which, to a good approximation, is comprised of rising and decreasing segments with linear time dependence throughout most of each half of the period, so that
is very close to a constant value rapidly changing sign over the vanishingly small interval around the turning points. Around the turning points, the solution is highly curved. Consequently, As E → ∞, x(t) tends to a periodic saw-tooth profile, the velocity,
,—to a periodic step-function, the acceleration,
,—to a periodic sequence of δ-functions and the phase-space plot,
vs. x(t),—to a rectangle.
A different effect of nonlinear violence emerges in some systems with polynomial nonlinearities, in which the coefficients of one or more of the lower powers in the polynomial grow with the maximal amplitude, a(E), so rapidly that they compete with the effect of the leading power of the polynomial when the solution, x(t), is close to a(E). For such systems, x(t) tends to a step function when a(E) grows indefinitely. The Duffing equation, in which the coefficient of the linear term grows with energy, is discussed in Section 2.3 because of its special physical significance. It was obtained by an ingenious transformation of the relativistic harmonic oscillator [46] [47] [51] [52] [53] .
The case, in which V(x) is not symmetric under reflection in x, and nonlinear violence emerges only near one turning point is discussed in Section 3.
The motion of a particle in a central field, V(r) (r is the distance of the particle from the origin) in (1 + 2) dimensions, when E → ∞ and the particle has non-zero angular momentum, opens the door to interesting phenomena in higher space dimensions. It is discussed in Section 4. The classical case of motion under a gravitational field, V(r) µ −1/r, follows a single “oscillation” with an infinite period of motion [36] . At high E, if the system has a non-zero angular momentum, l, the dominance of the latter in the effective potential leads to nonlinear violence around the single turning point, the smallest radius, rmin. If, instead, V(r) is chosen to grow rapidly at large r, the motion is periodic, with turning points at the smallest radius, rmin, and the largest radius, rmax. For E >> 1, the solution exhibits nonlinear violence at both turning points. At rmin, this behavior is caused by the dominance of the angular-momentum term in the effective potential (as E → ∞, rmin → 0). At rmax, this behavior is caused by the dominance of the potential, V(x), (as E → ∞, rmax → ∞). The boundary-layer characteristics of r(t) around rmin, leads to a rapid jump in the angle of rotation around the origin, φ, when crossing rmin, within a short time interval, which vanishes when E → ∞. This jump is followed by very slow variation of φ with time away from rmin. Summarizing comments are presented in Section 5.
In most systems discussed, when a turning point that is affected by nonlinear violence is crossed (crossing time t0), as E → ∞, the velocity jumps by a finite value, Δv. As a result, the acceleration tends to Δv ∙ δ(t − t0). For each system discussed, it is shown that, within the accuracy of calculational procedure, this jump is reproduced by the “inner solution” of the boundary layer-like analysis.
2. Nonlinear Violence in (1 + 1) Dimensional Oscillatory Systems
2.1. System Not Transformable to Boundary Layer Formulation
Consider an energy-conserving oscillatory system:
(3)
where V(x) → ∞ as x → ∞ and the maximal amplitude, a(E) tends to infinity as the total energy, E → ∞. To ensure oscillatory behavior, V(x), must be positive at both turning points. Nonlinear violence emerges in the vicinity of a turning point, x = a(E), if V(x) grows very rapidly in the vicinity of that point.
2.1.1. Reflection Invariant Potential—General Case
To see how a boundary layer-like behavior emerges when V(x) is symmetric under reflection in x, scaled coordinates will be used:
(4)
The motivation for the choice of the angle variable, θ, is that the motion is periodic with period T. As T can be found either in closed form, or computed to numerically to any desired accuracy (see discussion in Section 1), in the analysis presented here, T is treated as a known constant. Furthermore, as the oscillations vary in the range −a ≤ +a, to obtain a standard limiting profile form that applies to many systems, the coordinate, x(t) is replaced by the bounded scaled coordinate, η, −1 ≤ η ≤ +1, with period (2π) in θ. Equation (3) becomes a boundary value problem for η:
(5)
The expression for the period, T is:
(6)
Often the integral, I, defined in Equation (6) is O(1). For example, if the potential function obeys the requirement:
(7)
where C is a constant, independent of a, then
(8)
In many physical systems, C = 0, and then I → 1. When Equation (7) holds, T vanishes in proportion to (V(a))−1/2 = E−1/2.
Exploiting Equation (6), Equation (5) becomes:
(9)
When nonlinear violence dominates, Equation (9) leads to the emergence of behavior of the solution, which looks like the “inner” and “outer” solutions of boundary-layer problems. This happens if
(10)
, (11)
Note the dramatic difference between Equations (10) and (11). Equation (10) ensures that, as long as |η| < 1 (specifically, (1 − |η|) >> O(1/a)), namely, not too close to the turning points, the second term in Equation (9) becomes negligible, so that
is very small. Thus, Equation (11) ensures that, except near the turning points, Equation (9) tends to
. (12)
Equation (11) yields a linear profile for the “outer solution”:
. (13)
Consider the interval, 0 ≤ θ ≤ π, with initial conditions: η(0) = 1 and η'(0) = 0. Owing to symmetry and periodicity of the solution, one has ηouter(π/2) = 0. Also, although ηouter may miss the actual solution very close θ = 0 and π, up to errors of at most O(1/a), it must reach close to +1 as θ → 0 and -1 as θ → π: Hence, Equation (13) must be:
(14)
On the other hand, Equation (11) ensures that, very close to the turning point,
is very large, so that the maximum η(0) = +1 is very sharp. η(θ) varies rapidly over a very short range in θ. Hence, θ = 0 is the tip of a boundary layer. As η(0) = +1 is a point of a symmetric maximum for η(θ), approximate η(θ) for small θ by
. (15)
Choosing the nomenclature “inner” for this part of the solution is based on the observation presented below, that it varies rapidly over a narrow range in θ in the vicinity of θ = 0. This is just the characteristic of the inner solution in a boundary-layer problem. Replace η(θ) in Equation (9) by ηinner(θ) and expand the result around θ = 0 through θ4 to find:
. (16)
Requiring that the inner and outer solutions coincide and so do their derivatives at a matching point, θm, one obtains for θm, and for A (see Equation (13)):
. (17)
Equation (11) ensures that the range in θ, over which the “inner solution” evolves, is small and vanishes for a → ∞. Consequently, the matching point, θm of Equation (17), shrinks to zero and, concurrently, that A tends to (2/π) as a → ∞. With Equations (10) and (11) satisfied, as a → ∞, η(θ) tends to a periodic saw-tooth profile: Linear segments connect the turning points (η(2nπ) = +1, η((2n + 1)π) = −1), the velocity, η'(θ), tends to a periodic step-function profile, alternating between ±(2/π) and the phase-space plot, η'(θ) vs. η(θ), tends to a rectangle.
δ-function profile for acceleration The periodic step-function limit of the profile of η'(θ) as a → ∞ implies that, η''(θ), tends to a periodic sequence of δ-functions:
. (18)
The δ-functions arise from the cumulative change in velocity, η', across a turning point. Within a very a good approximation, the contribution of the “inner solution” generates the jump in η' as the turning point is crossed, as a → ∞. Using ηinner(θ) of Equation (14) through O(θ2) for η(θ) as well as Equations (15) and (16), yields the total change in the slope:
. (19)
Note that the result of Equation (19) does not depend on the asymptotic behavior of the period, T or the integral I (see Equation (6) and the following discussion).
2.1.2. Gaussian Potential
As an example to the statements made in Section 2.1.1, consider a Gaussian potential,
(20)
Equations (10) and (11) are satisfied for this potential as a → ∞.
Using the scaled coordinates of Equation (4), the period is given by:
(21)
For the error function, see [54] . As a grows, I tends to 1 rapidly. Already for a = 10, one finds I = 1.007. In addition, Equation (17) yields that the matching point, θm , is O(1/a2).
The equation for η(θ) is:
. (22)
Figures 1-3 present numerical solutions of Equation (22), showing, respectively, how η(θ) tends to a periodic saw-tooth profile, η'(θ)—to a periodic step-function and the phase plot, η'(θ) vs. η(θ), to a rectangle. Let us see how this happens.
Figure 1. Approach of η(θ), solution of Equation (22), to periodic saw-tooth profile as maximal amplitude, a, is increased.
Figure 2. Approach of η'(θ), solution of Equation (22), to periodic step-function profile as maximal amplitude, a, is increased.
Figure 3. Approach of η(θ) − η'(θ) phase-space plot of solution of Equation (22), to rectangle as maximal amplitude, a, is increased.
At the turning point |η| = 1 and |η'| = 0, energy conservation for Equation (22) yields:
. (23)
Equation (23) implies that max|η'(θ)| (attained at θ = π/2 and 3π/2, where η(θ) = 0), is:
. (24)
The limit of Equation (24) also holds for η'(θ) of Equation (23) as long as (1 − η(θ)2) = O(1):
. (25)
Thus, except near the turning points, where η(θ) ≈ ±1, the profile of the “outer solution” tends to a straight line, which, owing to the symmetry and periodicity of the solution, must pass through zero at θ = π/2:
. (26)
Only near the turning points, for (1 − η(θ)2) = O(1/a2), is this limit not attained, and the “inner solution” evolves. For η(0) = 1 and η'(0) = 0 at the turning point at θ = 0, the solution is maximal and symmetric in deviations around θ = 0. Hence, expand η(θ):
. (27)
Substituting Equation (27) in Equation (22), and expanding the equation through the O(θ2) yields:
. (28)
. (29)
Using Equations (21) and (27)-(29) yields the “inner solution”, whereas the “outer solution” is determined up to an overall coefficient:
. (30)
Requiring that the two solutions and their derivatives coincide at a matching angle, θm, in the vicinity of θ = 0, one obtains:
, (31)
. (32)
θm, the range, over which the “inner solution” evolves, shrinks to zero and A tends to (2/π) at an O(1/a2) rate.
In summary, as a → ∞, the solution for η(θ) tends to a periodic saw-tooth profile, which varies linearly from +1 at θ = 0, to −1 at θ = π and back to +1 at θ = 2π. η'(θ) tends to a periodic step function, with values of −(2/π) over 0 ≤ θ ≤ π and +(2/π) over π ≤ θ ≤ 2π.
δ-function profile for acceleration As in Section 2.1.1, η''(θ) tends to a periodic sequence of δ-functions with alternating signs (Equation (18)). Again, the jump in η'(θ) as the turning point is crossed is well approximated by the cumulative change of η'(θ) over the “inner solution”. Substituting ηinner(θ) of (24)-(29) for η(θ) through O(θ2) one obtains:
. (33)
This is the total change of the slope of the “outer solution” as the turning point is crossed.
2.1.3. Polynomial Nonlinearity: Constant Coefficients
This Section focuses on the case of a polynomial V(x):
. (34)
When all coefficients, a2n+1, in Equation (34), are constant, the leading power, x2N+1, dominates around the turning points in the limit of very high energy, E, as it grows like a2N+1. Still, the nonlinearity in Equation (34) does not obey conditions (10), (11) and (16) for fixed (2N + 1). The asymptotic saw-tooth profile of the solution emerges only in the double limit E → ∞ and N → ∞. As this case has been discussed in detail in [34] , the results are briefly reviewed here. Defining a small parameter,
, (35)
the period of oscillation, T, is found to be O(μN):
, (36)
where (2p + 1) is the power of the next-to the leading term.
Under the scaling of Equation (4)), Equation (34) becomes:
. (37)
In the double limit of a → ∞ (μ → 0) and N → ∞, for all |η(θ)| < 1, Equation (37) is reduced to
, (38)
leading to the emergence of the “outer solution”.
The boundary-layer-like analysis yields:
. (39)
Nonlinear violence emerges in the double limit of a → ∞ and N → ∞, leading to the same periodic saw-tooth profile for η(θ) as in previously discussed systems.
δ-function profile for acceleration
As in Sections 2.1.1 and 2.1.2, η''(θ) tends to a periodic sequence of δ-functions (Equation (18)). Again, the jump in η'(θ) generated by the δ-function is well approximated by cumulative change of η'(θ) over “ηinner” of Equation (39):
(40)
2.2. Relativistic Harmonic Oscillator: Transformable into Boundary-Layer Problem
2.2.1. Introductory Comments
In the extreme relativistic limit of the relativistic harmonic oscillator,
. (41)
x(t) tends to a periodic saw-tooth profile. Over the years, this has been demonstrated through numerical solutions of the Equation (41) [45] [46] [47] . A closed-form expression for the solution of Equation (41) [46] [51] [52] [53] , confirms the numerical results. For c, the speed of light in vacuum, there is little interest in the extreme relativistic limit of Equation (41): Near the turning points, the acceleration is so high, that General Relativistic effects become important. Hence, to perform an experiment, the oscillator must be placed in outer space far away from any large mass. Interest in the relativistic harmonic oscillator has recently emerged in low-temperature physics, in the oscillations of ions in a Bose-Einstein condensate [55] , where the speed of light is smaller than in vacuum by orders of magnitude.
To present numerical solutions of Equation (41) it pays to use dimensionless variables:
. (42)
Equation (41) is transformed into
(43)
Define
. (44)
An example of phase-space plots for
vs.
is shown in Figure 4 for β = 0.3, 0.7, 0.9999. As β → 1 (E → ∞), the plots tend to a rectangle. Such plots have been presented in many papers (see, e.g., [31] [32] [33] [37] [38] [39] .)
That the solution of Equation (41) tends in the extreme relativistic limit to the same limiting saw-tooth profile as in previous Sections may not be surprising, as
tends to ±c, the speed of light, throughout most of the period of oscillations. The characteristic that is common to the present system and the previously discussed systems, is that the solution of Equation (41) also exhibits a boundary-layer-like structure: most of the time, the solution is well approximated by the linear time dependence of an “outer solution”. Around each turning point, the solution of x(t) has a strongly curved profile over a very small fraction of the total period, T, a fraction that tends to zero as E → ∞. The curved part of the solution is very well approximated by an “inner solution”. Dominance of the nonlinearity near the turning points ensures that the asymptotic form of the solution of Equation (41) is identical to that found in Sections 2.1.1 - 2.1.3. This happens despite the fact that there are qualitative differences between the relativistic harmonic oscillator and the previously discussed systems. In the previous systems, one has:
, (45)
where E is the total energy. In the case of the relativistic harmonic oscillator, one has:
. (46)
Moreover, in the systems discussed in previous subsections, nonlinear violence emerged without resorting to the transformation of the dynamical equations into traditional boundary-layer systems of the form of Equation (2). In the case of Equation (41), nonlinear violence emerges in the extreme relativistic limit, through such a transformation.
2.2.2. Boundary Layer-Like Analysis
Energy conservation for Equation (41) yields:
. (47)
The transformation of Equation (4) yields a periodic boundary value problem:
Figure 4. Phase-space plots of solution of Equation (43). From inward to outward: β = 0.3, 0.7, 0.9999.
. (48)
η(θ) is (2π)-periodic in θ and bounded: |η(θ)| ≤ 1. The period, T, is given by [31] :
, (49)
where E(k) and K(k) are the complete elliptic integrals [54] :
. (50)
Applying energy conservation to Equation (48) and exploiting properties of the elliptic integrals [54] the maximal values of η and η' are found to be:
. (51)
To measure the proximity to the extreme relativistic limit, define a small parameter, μ:
. (52)
Expansion in power of μ, (for elliptic integrals, see [54] ) yields to lowest order in μ:
. (53)
The system exhibits the characteristics of a boundary layer problem (see Equation (2)): the highest derivative in the Equation (53) is multiplied by a small parameter.
To find the inner solution corresponding to Equation (53) near the turning point, θ = 0, (η(0) = 1, η'(0) = 0), observe that η' is very close to zero there. Replacing η' by zero in Equation (53), the O(μ2) approximation to Equation (53) becomes:
. (54)
For the boundary conditions at the turning point, Equation (54) is solved by
. (55)
To obtain the “outer solution” we focus on the vicinity of θ = (π/2). As the chosen turning point is at θ = 0, two observations hold near θ = (π/2): η(θ) ≈ 0 and η' ≈ (
). Substituting η(θ) = 0 in Equation (53), yields the approximate equation;
. (56)
Exploiting the leading-order term in
by expanding Equation (51) in powers of μ, yields:
. (57)
Using Equations (55) - (57), an O(μ2) approximation for the matching point, θ = θm, at which the slope of ηinner becomes equal to the constant slope of the ηouter, which prevails throughout most of the range 0 ≤ θ ≤ (π/2) one finds:
. (58)
Requiring that ηinner and ηouter coincide at θ = θm yields:
. (59)
δ-function profile for acceleration As in previous sections, η''(θ) tends to a periodic sequence of δ-functions with alternating signs (Equation (18)). Again, this is well approximated by the cumulative change of η'(θ) over “the inner solution” as the turning point is crossed. Using Equations (55) and (58), one obtains the total change in the slope of the “outer solution” as the turning point is crossed:
. (60)
In the extreme relativistic limit (μ → 0), the preceding analysis ensures that the linear “outer solution” describes the solution very well and the time interval around the turning points, where the solution is described by xinner tends to zero. As a result, the asymptotic profile of the solution is also a saw-tooth one, as in Section 2.1.
2.3. Polynomial Nonlinearity: Coefficients That Grow with Energy
A new vista opens up when coefficients of the non-leading monomials in Equation (34) grow with energy, i.e., if for some p < N,
, (61)
if, in addition, a2p+1 > 0 and a2N+1 < 0. More than one term in the polynomial participates in nonlinear violence, and owing to the opposite signs, they compete with one another. For a → ∞, an additional leading monomial shows up in Equation (38):
. (62)
Unlike the case of constant coefficients in a polynomial potential (Section 2.1.3), the additional requirement of N → ∞ is not necessary here. The characteristics of the solutions for different values of p and N are similar, provided α = 1. Hence, one example of physical significance is discussed, the Duffing equation with an energy-dependent liner term:
. (63)
Equation (63) was derived by an elegant transformation of the equation of motion of the relativistic harmonic oscillator, Equation (41), using proper time (time as measured on the moving mass rather than time measured by an external observer) and scaled, dimensionless coordinates [46] [47] [51] [52] [53] . The linear and cubic terms are comparable near the turning points, x = ±a(E). Employing Equation (4), the scaled equation is:
. (64)
The period T is given by [26] [27] [28] [29] :
, (65)
where properties of E(k) (see Equation (50)) are found in [54] . Using Equation (65), Equation (64) can be expressed as a traditional boundary-layer problem:
(66)
The highest derivative is multiplied by a small parameter, which varies slowly as a → ∞.
Equation (66) has a closed-form solution [40] [41] [42] [43] [44] [47] , given by:
. (67)
In Equation (67), sn is a Jacobi elliptic function [54] .
Applying energy conservation to Equation (66), or computing the derivative of η'(θ = π/2) through Equation (67), one obtains the maximal value of the derivative to be:
. (68)
The “inner solution” describes the rapid change around θ = π/2. It is obtained by solving for η(θ) in a power series in (π/2 − θ): Through O((π/2 − θ)2), it is given by:
. (69)
The “outer solution”, describing the moderate change around θ = 0 is found by solving Equation (66) around η(θ) = 1:
(70)
Substituting Equation (70) in Equation (66), the linearized equation for
becomes:
. (71)
Equation (71) is solved by
. (72)
The matching point, θm, is found by requiring
. (73)
One finds:
. (74)
The rate of approach as a → ∞ is slow. The matching point, θm, deviates from (π/2) by O(1/log2a). Figure 5 shows the explicit solution for η (Equation (67)), ηinner Equations (69) and (74), ηouter Equations (70) and (72) and the location of θm, Equation (74).
3. V(x) Not Symmetric under Reflection in x, (1 + 1) Dimensions
3.1. General Comments
If V(x) is not symmetric under reflection in x, the amplitudes at the two turning points are not equal up to a sign. Consequently, the analysis has to be performed separately close to any turning point, around which nonlinear violence prevails. The overall scaling of the solution (x = a ∙ η) in Equation (4) has to be replaced near each turning point by (x = xtur(E)×η), where xtur(E) is the amplitude at that turning point. Finally, while the transformation (t → θ) of Equation (4) yields
Figure 5. Solution of Equation (64). Full line: closed-form solution (Equation (67)); dotted line: outer solution (Equations (70) and (74)); dashed line: inner solution (Equations (69) and (72)).
(2π)-periodicity in θ, there is no symmetry between (0 ≤ θ ≤ π) and (π ≤ θ ≤ 2π). For example, if only one turning point exhibits nonlinear violence, T is a sum of two terms: a vanishingly small interval of time (spent around that turning point affected by nonlinear violence) and the rest of the time, spent over most of the oscillation.
Except for technical details, the analysis, specifically, the procedure for matching the “inner” and “outer” solutions, is the same as presented in Sections 2.1 - 2.4. Hence, instead of analyzing a general system, the analysis is carried out in Sections 3.2 - 3.3 for a specific example, where the solution is affected by nonlinear violence only around one turning point.
3.2. The System
Consider a harmonic oscillator with an exponential nonlinearity:
. (75)
The potential energy is asymmetric:
. (76)
For large E, at large negative x, V(x) of Equation (76) tends to that of a harmonic oscillator. Hence, the profile of oscillations is expected to be close to a trigonometric function. Nonlinear violence occurs only around the turning point, xmax(E) > 0, as there the potential energy grows rapidly. (xmax(E) → ∞ as E → ∞.) If the initial conditions are such the x(0) is that turning point, so that
= 0, turning points occur at t0 = 0, T, 2T, …
As E → ∞, for x(t) < 0, the exponential term in Equation (75) becomes rapidly negligible. Hence, the time spent over the x(t) < 0 branches of oscillations approaches π (half of a full sinusoidal oscillation). The time spent around the turning point, which is governed by nonlinear violence, is vanishingly small. Hence, T tends to π as E → ∞. For x(0) = 10, numerical integration yields T = 3.25.
Numerical Examples
The boundary-layer structure that emerges around xmax(E) is demonstrated Figure 6 and Figure 7 for an initial amplitude, x(0) = xmax(E) = 10. The details around a turning point are shown in Figure 8 and Figure 9. Near the turning point, x(t) is highly curved, soon evolving into a profile that is well approximated by linear time dependence. Over the same region,
varies rapidly between two constants, ±(
). This boundary-layer behavior is soon replaced by a trigonometric profile that dominates most of the period of oscillations.
Figure 6. x(t), solution of Equation (75). x(0) = 10,
.
Figure 7.
, solution of Equation (75). x(0) = 10,
.
Figure 8. x(t), solution of Equation (75) near positive turning point. x(0) = 10,
.
Figure 9.
, solution of Equation (75) near positive turning point. x(0) = 10,
.
3.3. Nonlinear Violence Analysis
To find an approximation for the solution of Equation (75) near the turning point, a = xmax(E) > 0, for large E, write x(t) as:
. (77)
Equation (75) becomes
. (78)
Energy conservation for Equation (75) yields
. (79)
To find α, substitute Equation (77) in Equation (75) and approximate the “inner solution” by:
. (80)
To leading order, one obtains
. (81)
As the solution evolves into a trigonometric profile once away from the turning point, we expect the “outer solution” to provide a good approximation (linear in time) for the full solution for a limited range in time. For this reason, scaling of the time variable with α of Equation (81), and a of Equation (79) will be employed also in the generation of the “outer solution”. Equation (75) becomes:
. (82)
For E → ∞ and (1 − |η(τ)|) = O(1)) > 0, Equation (82) tends to:
. (83)
This yields the “outer solution”:
. (84)
Matching of the “inner” and “outer” solutions by making their values and their derivatives equal at t = tm yields the matching point:
. (85)
For a turning point at t = t0, the “inner solution” is found to be
, (86)
and the outer solution:
. (87)
A comparison of the approximations provided by the “inner solution” and the “outer solution” with the numerical solution of Equation (75) around the turning point amplitude, x(t0) = xmax = 10 and t0 = 3.25, is shown in Figure 10. The “inner solution” well approximates the curved part of the full solution near the turning point and the “outer solution” does so (for a while) away from the turning point.
Velocity change across turning point Using Equations (83) and (84), the cumulative change of η'(t) over the diminishingly small interval, (t0 − tm, t0 + θm,), (θm → 0 for E → ∞) is found to coincide with the change in the slope of the “outer solution”, xout, as one crosses the turning point from left to right:
Figure 10. Comparison of quadratic and quartic approximations for “inner solution” (Equation (86), dashed and dotted lines, respectively) and “outer solution” (Equation (87), dot-dashed) with numerical solution of Equation (75) near t0 = 3.25 turning point. x(0) = 10,
.
, (88)
again, asymptotically yielding the overall change in x'(t) as the turning point is crossed.
4. (1 + 2) Dimensions: Motion in Central Force with Non-Zero Angular Momentum
At high E, the motion of a particle in a central field in two space dimensions offers new interesting boundary-layer phenomena if the angular momentum is not equal to zero. The evolution of the radius, r, of a unit mass particle is determined by [36] :
. (89)
Here V(r) is the potential and l is the conserved angular momentum:
(90)
φ is the angle of rotation of the moving mass around the origin.
4.1. Gravitational Potential
4.1.1. Analysis of Solution
For E < 0, the solution classical motion under a gravitational potential,
, (91)
is bounded and periodic. For E > 0, the solution visits the point of the minimal radius, rmin(E), once and wonders off to infinity [36] . Boundary-layer characteristics emerge in the vicinity of the smallest radius, rmin, when
if l ≠ 0.
4.1.2. Turning Point at rmin
Energy conservation for Equation (89) with V(r) of Equation (91),
, (92)
yields:
(93)
and
(94)
Hence, the “outer solution” is given by:
(95)
For the “inner solution”, as rmin is a local minimum around t = 0, write:
. (96)
Substituting Equation (96) in Equation (89) yields:
. (97)
Using Equations (95) - (97), matching the “inner” and “outer” solutions and their time derivatives at the matching time, tm, yields:
(98)
Thus, the matching point, tm, vanishes like E-1.
Velocity change across turning point, rmin Using Equations (90) and (93) - (98), the total change in the velocity
, around rmin over the time interval −tm ≤ t ≤ + tm, computed in terms of rinner, reproduces the full change in the slope of the “outer solution” from negative to positive linear slopes:
. (99)
Rapid rotation at turning point, rmin
The last interesting aspect is the variation of the angle of rotation, φ, in time as E → ∞. Close to the turning point, Equations (90), (94), (96) and (97) yield:
. (100)
Thus,
becomes infinite in the vicinity of the turning point as E → ∞. However, soon the “outer solution” takes over, for which:
. (101)
In summary, close to rmin, the angle, φ, varies rapidly between two constants over a vanishingly small time interval around t = 0. Away from that interval it varies very slowly with time up to the large turning point, rmax. In the limit E → ∞, φ exhibits a finite jump around the rmin. This point will be expanded in Section (4.2).
4.2. Positive Potential That Grows as r → ∞
4.2.1. Analysis of Solution
Consider in Equation (89) V(r), which grows rapidly at large r. To be specific choose:
. (102)
The motion is oscillatory as the total energy is:
. (103)
Nonlinear violence, leading to a boundary-layer phenomenon, emerges around the turning point at the maximal radius, rmax. For
:
. (104)
Nonlinear violence emerges also in the vicinity of rmin if l ≠ 0 and the contribution of the angular momentum term dominates:
. (105)
Equation (105) holds for V(r) of Equation (102). The analysis is identical to that presented in Section 4.1, hence, not repeated here.
Numerical examples are shown for E = 106 and l = 1. Figure 11 shows the numerical solution for r(t) of Equation (89) with the potential of Equation (102). The profile is almost linear in t, except for rapid bending close to the turning points. Figure 12 shows the numerical solution for
. The profile is almost constant, except for rapid change of sign close to the turning points. The rapid transition of r(t) from the highly curved “inner solution” to the linear profile of the “outer solution” near rmin and rmax is demonstrated in Figure 13 and Figure 14, respectively.
4.2.2. Turning Point at rmax
Energy conservation for Equation (89) with V(r) of Equation (102),
, (106)
yielding Equation (104) for rmax and
Figure 11. Numerical solution of Equation (89) for r(t) with exponential potential energy (Equation (102)), E = 106. Dashed lines represent rmin and rmax.
Figure 12. Numerical solution of Equation (89) for
with exponential potential energy (Equation (102)), E = 106 and l = 1. Dashed lines represent the constant values of
away from turning points.
Figure 13. Numerical solution of Equation (89) for r(t) with exponential potential energy (Equation (102)) in vicinity of rmin, E = 106 and l = 1.
Figure 14. Numerical solution of Equation (89) for r(t) with exponential potential energy (Equation (102)) in vicinity of rmax, E = 106 and l = 1.
. (107)
Hence, the “outer solution” is given by:
. (108)
For the “inner solution”, as rmax is a local minimum around t = 0, write:
. (109)
Substituting Equation (109) in Equation (89) with the potential of Equation (102) yields:
. (110)
Using Equations (108)-(110), matching the “inner” and “outer” solutions and their time derivatives at the matching time, tm, yields:
. (111)
Thus, the matching point, tm, vanishes as E-1/2.
Velocity change across turning point, rmax Using Equations (108) and (111), the total change in the velocity
around rmax as one crosses the turning point, is well reproduced by the contribution of rinner:
. (112)
Again, this is the total change in
.
4.2.3. Variation of Rotation Angle near Turning Points
Near rmin
The angular velocity,
, is shown in Figure 15 for the numerical solution of Equation (89) with the potential of (102). It is highly peeked near rmin and negligibly small elsewhere. Details are shown in Figure 16 and Figure 17 where
is plotted in the vicinity of each turning point. Using Equations (93) and (96)-(98), the change of
when crossing rmin over the time interval (−tm, +tm) is found to be:
Figure 15. Computation of
using numerical solution of Equation (89) with exponential potential energy (Equation (102)), E = 106 and l = 1.
Figure 16. Computation of
using numerical solution of Equation (89) with exponential potential energy (Equation (102)) in vicinity of rmin, E = 106 and l = 1.
Figure 17. Computation of
using numerical solution of Equation (89) with exponential potential energy (Equation (102)) in vicinity of rmax, E = 106 and l = 1.
(113)
If instead of rinner we use the numerical solution of Equation (89) for r(t), with l = 1, the potential of Equation (102), and E = 106 in Equation (113) the computed jump obtained is 1.5715. The exact result is that, as E → ∞, Δφ jumps across rmin by (π/2) = 1.5708 [36] . Both results are very good approximations deviating from (π/2) owing to the approximations employed in expressing rinner, router and tm.
Near rmax
Assuming that rmax occurs at t = 0, one finds in the vicinity of rmax:
. (114)
The integral of Equation (113) for rinner of Equation (114) reveals that the change in φ(t) across rmax in the time interval (−tm, + tm), with tm of Equation (114) vanishes as E → ∞:
, (115)
In summary, in agreement with what the numerical solution shows, φ(t) varies rapidly (a jump of (π/2) in the limit E → ∞) when crossing rmin and by a vanishingly small change when crossing rmax. Outside the short time intervals spent around rmin (Equation (93)) and rmax (Equation (104)), the total variation of
along each branch of the oscillation cycle is (3/4π), occurring at a very low rate, O(1/(E ∙ t2)).
Comment: Choice of the exponential potential of Equation (102), the contribution of which to Veff of Equation (89) becomes negligible at small r, was motivated by the role played by the angular momentum in generating rapid rotation around rmin. A detailed analysis (not presented here) shows that, if the exponential potential of Equation (89) is replaced by V(r) such that its contribution to Veff of Equation (89) dominates also at small r, then, depending on the details of V(r), rapid rotation around rmin may or may not emerge.
5. Concluding Comments
Nonlinear oscillations with high oscillation energy are of relevance in many areas of science and technology: communications science, electrical circuits, electronic devices, biological systems and more. A recent example has emerged in low-temperature physics, in the oscillations of ions in a Bose-Einstein condensate [55] , where the speed of light is smaller than in vacuum by orders of magnitude and a saw-tooth profile for the oscillations has been observed.
There is a large family of oscillatory systems whose solutions are affected by nonlinear violence, leading to boundary layer-like characteristics when the oscillation energy is extremely high. These characteristics emerge independently of whether the dynamical equation of the oscillator can or cannot be cast in the traditional form of boundary-layer problems.
Planar motion of a particle in a central field displays a new unique phenomenon, which opens the door to interesting behavior of oscillators in higher space dimensions. Boundary-layer characteristics emerge around the turning point of the largest radius owing to the growth of the potential energy. If the angular momentum, l, does not vanish, nonlinear violence may emerge around the turning point of the smallest radius, rmin. The latter is accompanied by a boundary-layer characteristic behavior of the angle that determines the position of the rotating particle in the plane. At large overall energy, the change in the angle when crossing rmin is very rapid and tends to a discrete jump of (π/2) as E → ∞.
In all the cases discussed here, the change of velocity and of rotation angle when a turning point is crossed are well approximated by computing them with the relevant “inner solution” over the time interval: (−tm ≤ t ≤ +tm).
In the case of (1 + 1)-dimensional energy conserving nonlinear oscillatory systems, the effect of nonlinear violence depends on the detailed structure of the equations of motion. If the potential energy of the system is not symmetric under space reflection, nonlinear violence may emerge around one or both turning points. If the potential energy is symmetric under reflection in x, these characteristics emerge at both turning points. The solution for x(t) then tends to a simple periodic profile: examples of saw-tooth and step-function profile were discussed.
The analysis discussed in this paper leads to the conclusion that, as E → ∞, the asymptotic profiles of solutions of a large family of oscillatory systems tend to similar standard forms, such as saw-tooth or step-function profiles. Hence, if such profiles are observed experimentally in oscillating physical systems, it may be difficult to decipher the equations of motion. A possible future research direction that will improve our understanding of high-energy nonlinear oscillations is the development of an expansion analysis around the asymptotic (saw-tooth or step-function) profiles, for the deviations of the actual solution from the asymptotic profile. This may be very important in deciphering the equations of motion governing the dynamical system, for which only experimental observations exist and theory has not been developed yet.
Author Responsibility
1) The author made all contributions to the conception or design of the work. No data acquisition was required. No new software was created;
2) The author drafted the work or revised it critically for important intellectual content;
3) The author approved the version to be published;
4) The author agrees to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.