Continuous Galerkin Petrov time discretization scheme is tested on some Hamiltonian systems including simple harmonic oscillator, Kepler’s problem with different eccentricities and molecular dynamics problem. In particular, we implement the fourth order Continuous Galerkin Petrov time discretization scheme and analyze numerically, the efficiency and conservation of Hamiltonian. A numerical comp...

The equations governing the flow of an electrically conducting, incompressible viscous fluid over an infinite flat plate in the presence of a magnetic field are investigated using the homotopy perturbation method (HPM) with Padé approximants (PA) and 4 order Runge–Kutta method (4RKM). Approximate analytical and numerical solutions for the velocity field and heat transfer are obtained and compar...

A B-spline finite element method is used to solve the equal width equation numerically. This approach involves a collocation method using quintic B-splines at the knot points as element shape. Time integration of the resulting system of ordinary differential equations is effected using the fourth order Runge–Kutta method, instead of the finite difference method, the resulting system of ordinary...

Despite the popularity of high-order explicit Runge–Kutta (ERK) methods for integrating semi-discrete systems of equations, ERK methods suffer from severe stability-based time step restrictions for very stiff problems. We implement a discontinuous Galerkin finite element method (DGFEM) along with recently introduced high-order implicit–explicit Runge–Kutta (IMEX-RK) schemes to overcome geometry...

We now begin the definition and construction of Runge-Kutta methods. These one–step methods are essentially always stable, but designing Runge– Kutta methods which are consistent to high order can be difficult. This theory is presented in Sec. We have already seen several examples of Runge–Kutta methods: explicit and implicit Euler, the implicit midpoint rule, the explicit midpoint rule with Eu...

The purpose of this study is to approximate the stream function and temperature distribution of the MHD flow in a laminar liquid film from a horizontal stretching surface. In this paper DTM-Padé method was used which is a combination of differential transform method DTM and Padé approximant. The DTM solutions are only valid for small values of independent variables. Comparison between the solut...

In this paper, an effective algorithm which is a combination of the spectral collocation method and the parametric iteration method (PIM) is proposed for simulating the solution of the chaotic Genesio system (CGS). Comparison with the fourth-order Runge–Kutta method (RK4) confirms the very high accuracy of the presented algorithm. The obtained results reveal that the proposed algorithm is a pro...

The numerical behavior of two implicit time-marching methods is investigated in solving two-dimensional unsteady compressible flows. The two methods are the second-order multistep backward differencing formula and the fourth-order multistage explicit first stage, single-diagonal coefficient, diagonally implicit Runge-Kutta scheme. A Newton-Krylov method is used to solve the nonlinear problem ar...