Integral equation has been one of the essential tools for various areas of applied mathematics. In this paper, we review different numerical methods for solving both linear and nonlinear Fredholm integral equations of second kind. The goal is to categorize the selected methods and assess their accuracy and efficiency. We discuss challenges faced by researchers in this field, and we emphasize the importance of interdisciplinary effort for advancing the study on numerical methods for solving integral equations.

The classification of exact solutions, including solitons and elliptic solutions, to the generalized $K(m,n)$ equation by the complete discrimination system for polynomial method has been obtained. From here, we find some interesting results for nonlinear partial differential equations with generalized evolution.

A numerical method is presented to obtain the approximate solutions of the fractional partial differential equations (FPDEs). The basic idea of this method is to achieve the approximate solutions in a generalized expansion form of two-dimensional fractional-order Legendre functions (2D-FLFs). The operational matrices of integration and derivative for 2D-FLFs are first derived. Then, by these matrices, a system of algebraic equations is obtained from FPDEs. Hence, by solving this system, the unknown 2D-FLFs coefficients can be computed. Three examples are discussed to demonstrate the validity and applicability of the proposed method.

It is shown that a strong solution of the Degasperis-Procesi equation possesses persistence property in the sense that the solution with algebraically decaying initial data and its spatial derivative must retain this property. Moreover, we give estimates of measure for the momentum support.

The initial-boundary value problem for a class of nonlinear wave equations system in bounded domain is studied. The existence of global solutions for this problem is proved by constructing a stable set and obtain the asymptotic stability of global solutions through the use of a difference inequality.

This paper discusses spectral and spectral element methods with Legendre-Gauss-Lobatto nodal basis for general 2nd-order elliptic eigenvalue problems. The special work of this paper is as follows. (1) We prove a priori and a posteriori error estimates for spectral and spectral element methods. (2) We compare between spectral methods, spectral element methods, finite element methods and their derived $p$-version, $h$-version, and $hp$-version methods from accuracy, degree of freedom, and stability and verify that spectral methods and spectral element methods are highly efficient computational methods.

The nonlinear hydroelastic waves propagating beneath an infinite ice sheet floating on an inviscid fluid of finite depth are investigated analytically. The approximate series solutions for the velocity potential and the wave surface elevation are derived, respectively, by an analytic approximation technique named homotopy analysis method (HAM) and are presented for the second-order components. Also, homotopy squared residual technique is employed to guarantee the convergence of the series solutions. The present formulas, different from the perturbation solutions, are highly accurate and uniformly valid without assuming that these nonlinear partial differential equations (PDEs) have small parameters necessarily. It is noted that the effects of water depth, the ice sheet thickness, and Young’s modulus are analytically expressed in detail. We find that, in different water depths, the hydroelastic waves traveling beneath the thickest ice sheet always contain the largest wave energy. While with an increasing thickness of the sheet, the wave elevation tends to be smoothened at the crest and be sharpened at the trough. The larger Young’s modulus of the sheet also causes analogous effects. The results obtained show that the thickness and Young’s modulus of the floating ice sheet all greatly affect the wave energy and wave profile in different water depths.

A numerical method for solving nonlinear Fredholm integrodifferential equations is proposed. The method is based on hybrid functions approximate. The properties of hybrid of block pulse functions and orthonormal Bernstein polynomials are presented and utilized to reduce the problem to the solution of nonlinear algebraic equations. Numerical examples are introduced to illustrate the effectiveness and simplicity of the present method.

Two strategies for estimating open boundary conditions (OBCs) with adjoint method are compared by carrying out semi-idealized numerical experiments. In the first strategy, the OBC is assumed to be partly space varying and generated by linearly interpolating the values at selected feature points. The advantage is that the values at feature points are taken as control variables so that the variations of the curves can be reproduced by the minimum number of points. In the second strategy, the OBC is assumed to be fully space varying and the values at every open boundary points are taken as control variables. A series of semi-idealized experiments are carried out to compare the effectiveness of two inversion strategies. The results demonstrate that the inversion effect is in inverse proportion to the number of feature points which characterize the spatial complexity of open boundary forcing. The effect of ill-posedness of inverse problem will be amplified if the observations contain noises. The parameter estimation problems with more control variables will be much more sensitive to data noises, and the negative effects of noises can be restricted by reducing the number of control variables. This work provides a concrete evidence that ill-posedness of inverse problem can generate wrong parameter inversion results and produce an unreal “good data fitting.”

This paper proposes a domain decomposition method for the coupled stationary Navier-Stokes and Darcy equations with the Beavers-Joseph-Saffman interface condition in order to improve the efficiency of the finite element method. The physical interface conditions are directly utilized to construct the boundary conditions on the interface and then decouple the Navier-Stokes and Darcy equations. Newton iteration will be used to deal with the nonlinear systems. Numerical results are presented to illustrate the features of the proposed method.

From rotational potential vorticity-conserved equation with topography effect and dissipation effect, with the help of the multiple-scale method, a new integro-differential equation is constructed to describe the Rossby solitary waves in deep rotational fluids. By analyzing the equation, some conservation laws associated with Rossby solitary waves are derived. Finally, by seeking the numerical solutions of the equation with the pseudospectral method, by virtue of waterfall plots, the effect of detuning parameter and dissipation on Rossby solitary waves generated by topography are discussed, and the equation is compared with KdV equation and BO equation. The results show that the detuning parameter $\alpha$ plays an important role for the evolution features of solitary waves generated by topography, especially in the resonant case; a large amplitude nonstationary disturbance is generated in the forcing region. This condition may explain the blocking phenomenon which exists in the atmosphere and ocean and generated by topographic forcing.

The aim of this paper is to apply OHAM to solve numerically the problem of harmonic wave propagation in a nonlinear thermoelasticity under influence of rotation, thermal relaxation times, and magnetic field. The problem is solved in one-dimensional elastic half-space model subjected initially to a prescribed harmonic displacement and the temperature of the medium. The HAM contains a certain auxiliary parameter which provides us with a simple way to adjust and control the convergence region and rate of convergence of the series solution. This optimal approach has a general meaning and can be used to get fast convergent series solutions of the different type of nonlinear fractional differential equation. The displacement and temperature are calculated for the models with the variations of the magnetic field, relaxation times, and rotation. The results obtained are displayed graphically to show the influences of the new parameters.

We investigate the effect of boundary slip on the transient pulsatile fluid flow through a vessel with body acceleration. The Fahraeus-Lindqvist effect, expressing the fluid behavior near the wall by the Newtonian fluid while in the core by a non-Newtonian fluid, is also taken into account. To describe the non-Newtonian behavior, we use the modified second-grade fluid model in which the viscosity and the normal stresses are represented in terms of the shear rate. The complete set of equations are then established and formulated in a dimensionless form. For a special case of the material parameter, we derive an analytical solution for the problem, while for the general case, we solve the problem numerically. Our subsequent analytical and numerical results show that the slip parameter has a very significant influence on the velocity profile and also on the convergence rate of the numerical solutions.

The homotopy decomposition method, a relatively new analytical method, is used to solve the 2D and 3D Poisson equations and biharmonic equations. The method is chosen because it does not require the linearization or assumptions of weak nonlinearity, the solutions are generated in the form of general solution, and it is more realistic compared to the method of simplifying the physical problems. The method does not require any corrected function or any Lagrange multiplier and it avoids repeated terms in the series solutions compared to the existing decomposition method including the variational iteration method, the Adomian decomposition method, and Homotopy perturbation method. The approximated solutions obtained converge to the exact solution as $N$ tends to infinity.

We investigate the spatial dynamics of a predator-prey system with Allee effect. By using bifurcation analysis, the exact Turing domain is found in the parameters space. Furthermore, we obtain the amplitude equations and determine the stability of different patterns. In Turing space, it is found that predator-prey systems with Allee effect have rich dynamics. Our results indicate that predator mortality plays an important role in the pattern formation of populations. More specifically, as predator mortality rate increases, coexistence of spotted and stripe patterns, stripe patterns, spotted patterns, and spiral wave emerge successively. The results enrich the finding in the spatial predator-prey systems well.

We introduce the rudiments of fractional calculus and the consequent applications of the Sumudu transform on fractional derivatives. Once this connection is firmly established in the general setting, we turn to the application of the Sumudu transform method (STM) to some interesting nonhomogeneous fractional ordinary differential equations (FODEs). Finally, we use the solutions to form two-dimensional (2D) graphs, by using the symbolic algebra package Mathematica Program 7.

The fractional complex transformation is used to transform nonlinear partial differential equations to nonlinear ordinary differential equations. The improved ($G^{\prime }/G$)-expansion method is suggested to solve the space and time fractional foam drainage and KdV equations. The obtained results show that the presented method is effective and appropriate for solving nonlinear fractional differential equations.

From the viewpoint of operator theory, we deal with the temperature dependence of the solution to the BCS gap equation for superconductivity. When the potential is a positive constant, the BCS gap equation reduces to the simple gap equation. We first show that there is a unique nonnegative solution to the simple gap equation, that it is continuous and strictly decreasing, and that it is of class $C^2$ with respect to the temperature. We next deal with the case where the potential is not a constant but a function. When the potential is not a constant, we give another proof of the existence and uniqueness of the solution to the BCS gap equation, and show how the solution varies with the temperature. We finally show that the solution to the BCS gap equation is indeed continuous with respect to both the temperature and the energy under a certain condition when the potential is not a constant.

A method is proposed to determine the numerical solution of system of linear Volterra integrodifferential equations (IDEs) by using Bezier curves. The Bezier curves are chosen as piecewise polynomials of degree n, and Bezier curves are determined on $[t_0, t_f ]$ by $n+1$ control points. The efficiency and applicability of the presented method are illustrated by some numerical examples.

We present the analytical and multishaped solitary wave solutions for extended reduced Ostrovsky equation (EX-ROE). The exact solitary (traveling) wave solutions are expressed by three types of functions which are hyperbolic function solution, trigonometric function solution, and rational solution. These results generalized the previous results. Multishape solitary wave solutions such as loop-shaped, cusp-shaped, and hump-shaped can be obtained as well when the special values of the parameters are taken. The $(G^{\prime }/G)$-expansion method presents a wide applicability for handling nonlinear partial differential equations.

We study a generalized double sinh-Gordon equation, which has applications in various fields, such as fluid dynamics, integrable quantum field theory, and kink dynamics. We employ the Exp-function method to obtain new exact solutions for this generalized double sinh-Gordon equation. This method is important as it gives us new solutions of the generalized double sinh-Gordon equation.

The aim of this study is to present the use of a semi analytical method called the optimal homotopy asymptotic method (OHAM) for solving the linear Fredholm integral equations of the first kind. Three examples are discussed to show the ability of the method to solve the linear Fredholm integral equations of the first kind. The results indicated that the method is very effective and simple.

The two-dimensional Volterra integral equations are solved using more recent semianalytic method, the reduced differential transform method (the so-called RDTM), and compared with the differential transform method (DTM). The concepts of DTM and RDTM are briefly explained, and their application to the two-dimensional Volterra integral equations is studied. The results obtained by DTM and RDTM together are compared with exact solution. As an important result, it is depicted that the RDTM results are more accurate in comparison with those obtained by DTM applied to the same Volterra integral equations. The numerical results reveal that the RDTM is very effective, convenient, and quite accurate compared to the other kind of nonlinear integral equations. It is predicted that the RDTM can be found widely applicable in engineering sciences.

A pressure-stabilized Lagrange-Galerkin method is implemented in a parallel domain decomposition system in this work, and the new stabilization strategy is proved to be effective for large Reynolds number and Rayleigh number simulations. The symmetry of the stiffness matrix enables the interface problems of the linear system to be solved by the preconditioned conjugate method, and an incomplete balanced domain preconditioner is applied to the flow-thermal coupled problems. The methodology shows good parallel efficiency and high numerical scalability, and the new solver is validated by comparing with exact solutions and available benchmark results. It occupies less memory than classical product-type solvers; furthermore, it is capable of solving problems of over 30 million degrees of freedom within one day on a PC cluster of 80 cores.

We introduced a relatively new operator called the triple Laplace transform. We presented some properties and theorems about the relatively new operator. We examine the triple Laplace transform of some function of three variables. We make use of the operator to solve some kind of third-order differential equation called “Mboctara equations.”

We examine possible approximate solutions of both integer and noninteger systems of nonlinear differential equations describing tuberculosis disease population dynamics. The approximate solutions are obtained via the relatively new analytical technique, the homotopy decomposition method (HDM). The technique is described and illustrated with numerical example. The numerical simulations show that the approximate solutions are continuous functions of the noninteger-order derivative. The technique used for solving these problems is friendly, very easy, and less time consuming.

We introduce a mixed finite element method for an elliptic equation modelling Darcy flow in porous media. We use a staggered mesh where the two components of the velocity and the pressure are defined on three different sets of grid nodes. In the present mixed finite element, the approximate velocity is continuous and the conservation law still holds locally. The LBB consistent condition is established, while the $L\mathrm{2}$ error estimates are obtained for both the velocity and the pressure. Numerical examples are presented to confirm the theoretical analysis.