Two-Step Hybrid Block Method for Solving Nonlinear Jerk Equations ()
1. Introduction
Jerk is the rate of acceleration change in physics; that is, the time derivative of acceleration, and as such the second velocity derivative, or the third time position derivative. The jerk is important in several mechanics and acoustics applications. The Jerk vector is here resolved into tangential-normal and radial-transverse components for planar motion, and the normal component is expressed as an affine differential invariant recognized as the aberrancy. Several geometric properties of the Jerk vector are established for plane motion using known aberrancy properties of curves [1] .
Nonlinear third-order differential equations, known as nonlinear Jerk equations, involving the third temporal displacement derivative, are of great interest in analyzing some structures which exhibit rotating and translating movements, such as robots or machine tools, where excessive Jerk leads to accelerated wear of transmissions and bearing elements, noisy operations and large contouring errors in discontinuities (such as corners) in the machining path [2] .
Many authors have studied the numerical solutions of the Jerk equation, harmonic balance approach to periodic solutions is used in [3] , in [4] they have written the high-order ordinary differential equation in terms of its differential invariants. New algorithm for the numerical solutions of nonlinear third-order differential equations was used jacobi-gauss collocation method in [5] , He’s variational iteration method was used in [6] for nonlinear Jerk equations. Modified harmonic balance method was used for nonlinear Jerk equations in [7] . In this paper, we consider a Jerk equation of the form
(1)
With initial conditions
. The most general function of Jerk, which is invariant in time-reversal and space-reversal and only nonlinear cubic, can be written as
(2)
where the parameters are
and
is constants. The current work is motivated by optimizing local truncation errors in order to find a hybrid point in a two-step block method to have the most accurate solution for Jerk Equation (1). We organize this paper as follows: The next section illustrates the method derivation, Section 3 presents the analysis of the method involving order four, Section 4 presents the numerical examples showing the productivity of the new technique when it is contrasted with the different strategies proposed in the scientific writing.
2. Methodology
To derive two steps hybrid block method with one off-step point we have a polynomial of degree 6 as follows:
(3)
with third derivative given by
(4)
Substituting Equation (4) into Equation (3) gives
(5)
By interpolating Equation (3) at
and collocating Equation (5) at
we get a system of equations written in the matrix form
,
Solving the above system (by using Gaussian elimination method for solving the system of linear equations [8] ) gives us the coefficients of the polynomial
.
By making the substitution
, the polynomial in Equation (3) may be written in the form:
(6)
where
(7)
Now, by evaluating the solution approximation at the point
:
(8)
Assess the approximation at the point of first derivative
:
(9)
Evaluate the second derivative approximation at the point
:
(10)
We choose to optimize the local truncation errors in the Equation (8), Equation (9) and Equation (10) to determine appropriate values for r Equation (8), Equation (9) and Equation (10). This choice at the end of the block
,
and
, which result respectively in
(11)
Determine the r values equating to zero the
coefficients in the local truncation error formulas in Equation (11), and we obtain the system
(12)
By solving Equation (12) and substituting r’s in Equation (11) we can choose the value r that gives the least truncation errors
Substituting r in the local truncation errors gives:
To obtain a two-step hybrid block solution method for solving Equation (1), we evaluate
at the points
, we get the following block of six equations
(13)
3. Characteristics of the Method
This section is presented the basic properties of the main method and analyzed it to establish their validity.
3.1. Order of the Method
We can rewrite the hybrid block method in the form
(14)
where
are matrices of coefficients of dimensions
,
,
,
,
.
If
is a sufficiently differentiable function, the linear difference operator
associated with the implicit two-step block hybrid method is considered in Equation (8), Equation (9), Equation (10), Equation (13), that is given
(15)
where the
and
are respectively the vector columns of the matrices
and V.
Expanding
and
in Taylor series about
we get
(16)
with
and
Note that the proposed method has order
at least [9] .
3.2. Zero Stability
We can write the method as a vector As
in Equation (14), we can write the method as a vector form.
where
The first characteristic polynomial
, the roots of polynomials are
. Hence the block method is zero -stable [9] .
3.3. Consistency
The block method has order
, in case
, this is a sufficient condition to be consistent with the associated block method [10] .
3.4. Convergence
We can establish the convergence of the two-step with three points hybrid block method if and only if it is consistent and zero stable [11] .
3.5. Region of Absolute Stability
As we mentioned earlier, zero-stability is a concept of the numerical method behavior for
. To decide whether a numerical method will produce good results with a given value of
, we need a concept of stability that is different from zero-stability. In most numerical methods intended to solve problems of third order, the stability properties are usually analyzed by considering the linear equation given by the Dalquist test [11] .
(17)
This problem has bounded solutions for
that tend to zero for
. We will define the region where the numerical method reproduces the manner of the exact solutions. Let us explain the procedure for obtaining such a region. Our method has nine equations in which there are four different terms of first derivatives:
, and second derivative:
and one intermediate values
. Let us depict the procedure to gain such a region [12] , We eliminated these terms from the equations system by using mathematica, and get a recurrence equation in the terms
. This recurrence equation reads
(18)
where
, and
(19)
Through its characteristic equation, we study the extent boundedness of their solutions to define the region of stability. The root causes of the characteristic equation
(20)
are
The roots of the characteristic equation must be less than 1, for the method to be stable. The stability region for the method has shown in Figure 1.
4. Numerical Examples
Example 1: Consider the Jerk equation in the following form:
(21)
With initial conditions:
, and exact solution:
where
.
Tables 1-3 show the absolute Errors at
and
respectively. Example 1 was solved in [13] , the comparison between our method and result in [13] is proposed in Table 4 and Table 5. It is seen in Figures 2-4 the method gives a good approximation in the interval
and
respectively.
Example 2: Consider the Jerk equation in the following form:
(22)
With initial conditions:
, and exact solution:
where
.
Figure 2. The solution of Example 1 at
,
.
Figure 3. The solution of Example 1 at
,
.
Figure 4. The solution of Example 1 at
,
.
Table 1. Absolute error of numerical solutions at
.
Table 2. Absolute error of numerical solutions at
.
Table 3. Absolute error of numerical solutions at
.
Table 4. Comparison between the numerical solution in our method and the method in [13] ,
.
Table 5. Comparsion between the numerical solution in our method and the method in [13] ,
.
Tables 6-8 show the absolute Errors at
and
respectively. Example 2 was solved in [13] , the comparison between our method and result in [13] is proposed in Table 9 and Table 10. It is seen in Figures 5-7 the method gives a good approximation in the interval
and
respectively.
Figure 5. The solution of Example 2 at
,
.
Figure 6. The solution of Example 2 at
,
.
Figure 7. The solution of Example 2 at
,
.
Table 6. Absolute error of numerical solutions at
.
Table 7. Absolute error of numerical solutions at
.
Table 8. Absolute error of numerical solutions at
.
Table 9. Comparsion between the numerical solution in our method and the method in [13] ,
.
Table 10. Comparison between the numerical solution in our method and the method in [13] ,
.
5. Conclusion
A two-step with one hybrid point was proposed and proceeded as a self-starting method for solving nonlinear Jerk equations. We considered one hybrid point and specified for approximation after optimizing local truncation errors related to the main formula. Therefore, our method’s good convergent and stability properties make it attractive for the numerical solution of nonlinear problems. The method presented is zero stable, consistent and convergence of four-algebraic order. The numerical results and figures show their efficiency and precision compared to other methods in the literature.