# Parametric Identification for the Biased Ship Roll Motion Model Using Genocchi Polynomials

Authors: G. Swaminathan; G. Hariharan; S. A. Mohiuddine; Kandhasamy Tamilvanan; Masho Jima Kabeto
Journal: Journal of Mathematics
 (2022)

---

## Abstract

Roll motion is one of the key motions related to a vessel’s dynamic stability. It is essential for the dynamic stability of ships in the realistic sea. For this research study, we have investigated the parameters involved in damping of the ship. In general, mathematical modelling of the rolling response of a ship can be formulated by the linear, nonlinear, and fractional differential equations because the amplitude of oscillation is increased. An efficient Genocchi polynomial approximation method (GPAM) is successfully applied for the biased ship roll motion model. The basic idea of the collocation method together with the operational matrices of derivatives used for nonlinear differential equation and convert it into a system of algebraic equations. The convergence and error analysis of the proposed method are also discussed. A few numerical experiments are carried out for some specific and important types of problems including the biased roll motion equations. The results are compared to those produced using the Legendre wavelet method (LWM) and the homotopy perturbation method (HPM). It is observed that the proposed spectral algorithm is robust, accurate, and easy to apply.


---

## Body

## 1. Introduction

Several research papers have been studied on the ship roll motion models because they are highly nonlinear differential equations. It is an essential for the dynamic stability of ships in realistic sea. For this research study, GPAM has been applied to estimate the parameters involved in damping of the ship. In general, mathematical modelling of the rolling response of a ship can be formulated by the nonlinear differential equations because the amplitude of oscillation is increased. In ship roll motion models, the nonlinear term occurs due to the damping and restoring moments.In pioneering work, Spyrou [1] established the ship roll motion model for large regular waves. Several research papers have been published related to ship dynamics models [2, 3]. Nonlinear mathematical formulation of the ship roll model is discussed in [4]. In general, the linear differential model of the motion has been violated by the nonlinear features of the motion.The theoretical models of ship dynamics having nonlinear terms such as damping and restoring moments are discussed in [2, 5]. A few mathematical methods had been established to investigate the approximate solutions of nonlinear models.During the last decades, Appell polynomials are employed in analytic number theory as well as asymptotic approximation theory. They discovered differential equations of Bernoulli and Euler polynomials as a special case. For Appell polynomials, a collection of finite-order differential equations were established by [6]. The 2D Bernoulli and 2D Euler polynomials were introduced for solving differential, partial differential, and integrodifferential equations. Hermite polynomials based Appell polynomial were presented by [7] and some of the references related with Appell polynomials are [6, 8–15]. Several analytical and approximation schemes are proposed to solve a system of differential equations in nonlinear and fractional order [16, 17]. Hariharan et al. [18, 19] had developed the wavelet method for solving reaction-diffusion problems and integral equations.Nonlinear problem research has advanced at a rapid and intensive rate during the forecast period. For example, approaches based on operational matrices for integration of large numbers of polynomials and functions can be mentioned. Over the last four decades, numerical techniques based on operational matrices of integration (especially for orthogonal polynomials and functions) have gained a lot of attention for dealing with variational problems.The authors of [20, 21] have recognised the significance of using operational approaches in the study of special functions and their applications. The majority of the attention is focused on operational identities related to ordinary and multivariable versions of Hermite and Laguerre polynomials.A few researcher developed Frobenius–Euler and Genocchi polynomials and its application introduced to solve nonlinear differential equations. The novel identities of umbral calculus application majorly using the Genocchi numbers and polynomials which were derived by [8]. Recently, Swaminathan et al. [22] established the Genocchi polynomial approximation method for Bratu-type differential equations. Some of the researcher solved fractional-order differential equations [23] and nonlinear reaction-diffusion problems [24]. And also he derived some more special polynomials and function arising in application of fractional calculus [23, 25].This study provides a numerical analysis of ship roll motion in terms of concentrated parameters within the time domain using a Genocchi polynomial operational matrix method.This study is summarized as follows. Ship roll motion in biased type is discussed in Section2. Genocchi polynomial and its properties are given in Section 3. In Section 4, a few illustrative examples are presented to show the efficiency and accuracy of the proposed Genocchi polynomial method. Conclusion is given in Section 5.


## 2. Ship Roll Motion in Biased Type

The biased roll motion is described theoretically in [26]:(1)θ″+βθ′+θ1−θ1+aθ=Fsinωt,where a is a measure of the amount of bias, and its parameter lies between 0 and 1. Due to large angles of roll, the hull may affected by waves, winds, or cargo shifting. The perspective of biased roll motion is to approximate the ship. There is no bias in the Falzarano model [27]. If a=1, then the restoring moment is a cubic polynomial of symmetric type. If a=0, then the restoring moment becomes the second-order polynomial. Here, we consider the mathematical model in extreme case of bias. The ship can only capsize from one side. Then, it is called Helmholtz–Thompson equation, and Spyrou et al. [26] pointed out the value of a is closer to 1 in realistic situation. Therefore, 1−a is a small number. Equation (1) becomes(2)θ″+θ−θ3+βθ′=1−aθ2−θ3+Fsinωt.


## 3. Genocchi Polynomial and Its Properties

Genocchi polynomials have been used in various fields of mathematics. The classical Genocchi polynomials is denoted byGmx, and it is defined as [8, 28](3)2textet+1=∑m=0∞Gmxtnn!,t<π,where Gmx is the mth degree Genocchi polynomial and is given by(4)Gmx=∑p=0mmpGm−pxp.Gp is called the Genocchi number. Differentiate (4), with respect to x; then,(5)dGmxdx=mGm−1x,m≥1,Gm1+Gm0=0,m>1.The above equality (6) is an important property. Let Gx=G1x,G2x,G3x,…,GNx be the Gennochi vectors, then its derivative G′x, by the reference of equation (5), can be expressed by(6)G′xT=SGTx,where(7)G′xT=G1′xG2′xG3′xG4′x…GR−1′xGR′x,E=000…000200…000030…000004…000⋮⋮⋮…⋮⋮⋮000…R−100000…0R0,GxT=G1xG2xG3xG4x…GR−1xGRx,where S operational matrix of derivative of order R×R. Based on that procedure, the pth derivative of Gx can be expressed as(8)Gpx=GxETp.


### 3.1. Function Approximation

LetM=L20,1 [23] (section 2.4, page 3039) and G1x,G2x,…,GRx⊂M be the set of all Genocchi polynomials, and the SpanG1x,G2x,…,GRx=X. Consider g be an arbitrary element of M, since X is a finite dimensional vector space, and g∗ be an unique better approximation in R such that(9)g−g∗2≤g−x2∀x∈X.Sinceg∗∈N, then there exist an unique coefficients c1,c2,…,cR such that(10)g≈g∗=∑r=1RcrGrx=CGx,where(11)C=c1,c2,…,cRT,Gx=G1x,G2x,…,GRx.The coefficientscn can be derived using the following lemma.Lemma 1.
Letg∈M=L20,1 be an arbitrary function approximated by the truncated Genocchi series ∑r=1NcrGrx; then, the coefficients cr, for r=1,2,…,N, can be calculated from the following relation:(12)cr=12r!gr−11+gr−10,where gr−1x denotes the r−1th derivative of g.Proof.
See Lemma 4.1, page 2125 in [29].Theorem 1.
Considergx∈C∞0,1 and g∗x be the approximated value of gx by using Genocchi polynomials; then, the error bound can be expressed as(13)errorgx∞≤1R!GRFR,where GR and FR denote the maximum value of GRx and gr−1x for all x∈0,1, respectively.Proof.
The proof is obvious by using the above lemma. See [29] (Theorem 2, page 2125).


## 3.1. Function Approximation

LetM=L20,1 [23] (section 2.4, page 3039) and G1x,G2x,…,GRx⊂M be the set of all Genocchi polynomials, and the SpanG1x,G2x,…,GRx=X. Consider g be an arbitrary element of M, since X is a finite dimensional vector space, and g∗ be an unique better approximation in R such that(9)g−g∗2≤g−x2∀x∈X.Sinceg∗∈N, then there exist an unique coefficients c1,c2,…,cR such that(10)g≈g∗=∑r=1RcrGrx=CGx,where(11)C=c1,c2,…,cRT,Gx=G1x,G2x,…,GRx.The coefficientscn can be derived using the following lemma.Lemma 1.
Letg∈M=L20,1 be an arbitrary function approximated by the truncated Genocchi series ∑r=1NcrGrx; then, the coefficients cr, for r=1,2,…,N, can be calculated from the following relation:(12)cr=12r!gr−11+gr−10,where gr−1x denotes the r−1th derivative of g.Proof.
See Lemma 4.1, page 2125 in [29].Theorem 1.
Considergx∈C∞0,1 and g∗x be the approximated value of gx by using Genocchi polynomials; then, the error bound can be expressed as(13)errorgx∞≤1R!GRFR,where GR and FR denote the maximum value of GRx and gr−1x for all x∈0,1, respectively.Proof.
The proof is obvious by using the above lemma. See [29] (Theorem 2, page 2125).


## 4. Illustrative Examples

Example 1.
Consider the equations with initial conditions:(14)θ″+1.888267θ−0.081886cos0.527t=0,t∈0,1,θ0=1,θ0=0.
By using the procedure of Section3 for R=12. The general collocation points of Gennochi polynomial is ti=i/R−2,i=1,2,…,R−2. We find the approximated values of θ″t,θ′t, and θt was substituted in (15); we obtain(15)GtET2C+1.888267GtC−0.081886cos0.527t=0,with the initial conditions,(16)G0C=1,G0ETC=0.
Apply the collocation point (16) at ti=i/10; we get ten different algebraic equations. The obtained equations are solved together with (17) for the values of cj,j=1,2,…,12, and we obtain(17)c1=0.6340705,c2=−0.3659295,c3=−0.2439530,c4=c5=c6=c7=c8=c9=c10=c11=c12=0.
Hence,(18)θt=1−0.731859t2.
The Genocchi polynomial results are compared with HPM and LWM in Table1. The approximate solution has been illustrated in Figure 1.Table 1
Comparative Study of Genocchi polynomial, LWM, and HPM for roll angle with respect to the timet (in sec).
TimeEx:1Ex:2Ex:3tGenocchi polynomialLWMHPMGenocchi polynomialLWMHPMGenocchi polynomialLWMHPM0.01111111110.10.99268140.99268140.99308140.99929890.99929890.99933590.99931050.99931050.99933650.20.97072560.97072560.97232560.99719560.99719560.99734360.99724230.99724230.99734630.30.93413260.93413260.93773260.99369020.99369020.99402320.99379530.99379530.99402930.40.88290250.88290240.88930240.98878260.98878260.98937460.98896950.98896950.98938550.50.81703520.81703500.82703500.98247290.98247290.98339790.98276480.98276480.98341480.60.73653070.73653040.75093040.97476100.97476100.97609300.97518130.97518140.97611730.70.64138900.64138860.66098860.96564690.96564690.96745990.96621910.96621910.96749310.80.53161020.53160960.55720960.95513070.95513070.95749870.95587800.95587800.95754200.90.40719420.40719340.43959340.94321230.94321230.94620930.94415810.94415810.94626411.00.26814100.26814040.30814020.92989180.92989180.93359180.93105940.93105950.9336594Figure 1
Comparative study of Genocchi polynomial, LWM, and HPM for roll angle with respect to the timet (in sec) for Ex.1.Example 2.
Consider the equation with initial conditions:(19)θ″+1.888267θ−0.044582cos0.527tθ=0,t∈0,1,θ0=1,θ′0=0.
By Section3 procedure, this example have been done with N=12.
After simplification, we obtain(20)θt=1−0.0701082t2.
The Genocchi polynomial results are compared with HPM and LWM in Table1. The approximate solution has been illustrated in Figure 2.Figure 2
Comparative study of Genocchi polynomial, LWM, and HPM for roll angle with respect to the timet (in sec) for Ex.2.Example 3.
Consider the equation with initial conditions:(21)θ″+0.034480θ′+0.185761θ−0.044582cos0.527tθ=0,t∈0,1,θ0=1,θ′0=0.
By Section3 procedure, this example has been done with N=12.
After simplification, we obtain(22)θt=1−0.0689406t2.
The Genocchi polynomial results are compared with LWM and HPM in Table1. The approximate solution has been illustrated in Figure 3.Figure 3
Comparative study of Genocchi polynomial, LWM, and HPM for roll angle with respect to the timet values (in sec) for Ex.3.Example 4.
Consider the equation with initial conditions:(23)θ″+0.01θ′+θ−0.025θ2−0.975θ3−0.0195sin0.527t=0,t∈0,1,θ0=1,θ′0=0.
By Section3 procedure, this example has been carried out with N=12.
After simplification, we obtain(24)θt=1+0.00316317t2.
The Genocchi polynomial results are compared with LWM and HPM in Table1. The approximate solution has been illustrated in Figure 4.Figure 4
Comparative study of Genocchi polynomial, LWM, and HPM for roll angle with respect to the timet values (in sec) for Ex.4.Example 5.
Consider the equation with initial conditions:(25)θ″+0.01θ′+θ−0.05θ2−0.95θ3−0.0195sin0.527t=0,t∈0,1,θ0=1,θ′0=0.
By using Section3 procedure, this example have been carried out with N=12.
After simplification, we obtain(26)θt=1+0.00315087t2.
The Genocchi polynomial results are compared with LWM and HPM in Table2. The approximate solution has been illustrated in Figure 5.Table 2
Comparative Study of Genocchi polynomial, LWM, and HPM for roll angle with respect to the timet (in sec.).
TimeEx:4Ex:5Ex:6Ex:7tGenocchi polynomialLWMHPMGenocchi polynomialLWMHPMGenocchi polynomialLWMHPMGenocchi polynomialLWMHPM0.01111111111110.11.0000311.0000311.0000281.0000311.0000311.0000281.0000321.0000321.0000341.0000321.0000321.0000340.21.0001261.0001261.0001141.0001261.0001261.0001141.0001301.0001301.0001381.0001301.0001301.0001380.31.0002841.0002841.0002571.0002831.0002831.0002571.0002931.0002931.0003111.0002921.0002921.0003100.41.0005061.0005061.0004581.0005041.0005041.0004571.0005221.0005221.0005541.0005201.0005201.0005520.51.0007901.0007901.0007151.0007871.0007871.0007151.0008161.0008161.0008661.0008121.0008121.0008620.61.0011381.0011381.0010301.0011341.0011341.0010291.0011751.0011751.0012471.0011701.0011701.0012420.71.0015491.0015491.0014021.0015431.0015431.0014011.0015991.0015991.0016971.0015931.0015931.0016910.81.0020241.0020241.0018321.0020161.0020161.0018301.0020891.0020891.0022171.0020811.0020811.0022090.91.0025621.0025621.0023191.0025521.0025521.0023171.0026441.0026441.0028061.0026331.0026331.0027951.01.0031631.0031631.0028631.0031501.0031501.0028601.0032641.0032641.0034641.0032511.0032511.003451Figure 5
Comparative Study of Genocchi polynomial, LWM, and HPM for roll angle with respect to the timet (in sec) for Ex.5.Example 6.
Consider the equation with initial conditions:(27)θ″+0.05θ′+θ−0.025θ2−0.975θ3−0.0195sin0.527t=0,t∈0,1,and with initial conditions,(28)θ0=1,θ′0=0.
By using Section3 procedure, this example have been carried out with N=12.
After simplification, we obtain(29)θt=1+0.00326487t2.
The Genocchi polynomial results are compared with LWM and HPM in Table2. The approximate solution has been illustrated in Figure 6.Figure 6
Comparative Study of Genocchi polynomial, LWM, and HPM for roll angle with respect to the timet (in sec) for Ex.6.Example 7.
Consider the equation with initial conditions:(30)θ″+0.05θ′+θ−0.05θ2−0.95θ3−0.0195sin0.527t=0,t∈0,1,which subject to(31)θ0=1,θ′0=0.
By using Section3 procedure, this example have been carried out with N=12.
After simplification, we obtain(32)θt=1+0.00325179t2.
The Genocchi polynomial results are compared with LWM and HPM in Table2. The approximate solution has been illustrated in Figure 7.
In the tables, we have introduced the Genocchi polynomial-based spectral methods for the estimation of the parameters in an equivalent ship roll motion model using only the roll motion response (roll angle and roll velocity), which can be easily identified for a ship sailing at sea. There are several advantages of the proposed Genocchi polynomial approach method by comparison to other existing methods in the literature. The absolute error has been calculated for all the above examples, and it was listed in Table3.Figure 7
Comparative Study of Genocchi polynomial, LWM, and HPM for roll angle with respect to the timet (in sec) for Ex.7.Table 3
Absolute error of Genocchi polynomial for roll angle with respect to the timet (in sec).
TimeEx.1Ex.2Ex.3Ex.4Ex.5Ex.6Ex.70.000000000.100.00000370.0000260.0000030.0000030.0000020.0000020.200.00014800.0001040.0000080.0000120.0000080.0000080.30.00360.0003330.0002340.0000270.0000260.0000180.0000180.40.006390.0005920.0004160.0000480.0000470.0000320.0000320.50.009990.0009250.000650.0000750.0000720.000050.000050.60.014390.0013320.0009360.0000920.0001050.0000720.0000720.70.0196080.0018130.0012740.0001470.0001420.0000980.0000980.80.02559940.0023680.0016640.0001920.0001860.0001280.0001280.90.03239920.0029970.0021060.0002430.0002350.0001620.0001621.00.03999920.00370.00260.00030.0002900.00020.0003


## 5. Conclusion

In this study, the nonlinear equations in the roll motion of ships are solved by the Genocchi polynomial operational matrix approach with respect to the suitable collocation points. We successfully applied the proposed Genocchi polynomial operational matrix approach for estimating the ship roll angle. It has been observed that the proposed method is reliable to solve the biased ship roll motion equations approximately. Also, the obtained results have been validated with the results of HPM and LWM. It can also be concluded that the Genocchi polynomial operational matrix approach will be powerful for finding the approximation solutions for a large class of nonlinear differential equations. To avoid the computation complexity, we impose the collocation points to form algebraic equations. Satisfactory agreement with the results and this method is computationally attractive and effective while comparing with other methods. The advantage of this method is that highly accurate approximate solutions are achieved using a few numbers of terms of the approximate expansion. Comparative study of Genocchi polynomial, LWM, and HPM for roll angle and its figure give the accuracy by increasing the value of R. The proposed method can also be extended to solve other types of ODEs and PDEs arising in science and engineering.



---

*Source: 7918725-2022-05-14.xml*