Search results

The purpose of this paper is to solve the second-order nonlinear boundary value problem with nonlinear boundary conditions by an iterative numerical method.

The authors introduce eigenfunctions as test functions, such that a weak-form integral equation is derived. By expanding the numerical solution in terms of the weighted eigenfunctions and using the orthogonality of eigenfunctions with respect to a weight function, and together with the non-separated/mixed boundary conditions, one can obtain the closed-form expansion coefficients with the aid of Drazin inversion formula.

When the authors develop the iterative algorithm, removing the time-varying terms as well as the nonlinear terms to the right-hand sides, to solve the nonlinear boundary value problem, it is convergent very fast and also provides very accurate numerical solutions.

Basically, the authors’ strategy for the iterative numerical algorithm is putting the time-varying terms as well as the nonlinear terms on the right-hand sides.

Starting from an initial guess with zero value, the authors used the closed-form formula to quickly generate the new solution, until the convergence is satisfied.

Through the tests by six numerical experiments, the authors have demonstrated that the proposed iterative algorithm is applicable to the highly complex nonlinear boundary value problems with nonlinear boundary conditions. Because the coefficient matrix is set up outside the iterative loop, and due to the property of closed-form expansion coefficients, the presented iterative algorithm is very time saving and converges very fast.

The purpose of this paper is to propose the generalized integral transform technique (GITT) to the solution of convection-diffusion problems with nonlinear boundary conditions by employing the corresponding nonlinear eigenvalue problem in the construction of the expansion basis.

The original nonlinear boundary condition coefficients in the problem formulation are all incorporated into the adopted eigenvalue problem, which may be itself integral transformed through a representative linear auxiliary problem, yielding a nonlinear algebraic eigenvalue problem for the associated eigenvalues and eigenvectors, to be solved along with the transformed ordinary differential system. The nonlinear eigenvalues computation may also be accomplished by rewriting the corresponding transcendental equation as an ordinary differential system for the eigenvalues, which is then simultaneously solved with the transformed potentials.

An application on one-dimensional transient diffusion with nonlinear boundary condition coefficients is selected for illustrating some important computational aspects and the convergence behavior of the proposed eigenfunction expansions. For comparison purposes, an alternative solution with a linear eigenvalue problem basis is also presented and implemented.

This novel approach can be further extended to various classes of nonlinear convection-diffusion problems, either already solved by the GITT with a linear coefficients basis, or new challenging applications with more involved nonlinearities.

This paper aims to develop an efficient numerical method for nonlinear transient heat conduction problems with local radiation boundary conditions and nonlinear heat sources.

Based on the physical characteristic of the transient heat conduction and the distribution characteristic of the Green’s function, a quasi-superposition principle is presented for the transient heat conduction problems with local nonlinearities. Then, an efficient method is developed, which indicates that the solution of the original nonlinear problem can be derived by solving some nonlinear problems with small structures and a linear problem with the original structure. These problems are independent of each other and can be solved simultaneously by the parallel computing technique.

Within a small time step, the nonlinear thermal loads can only induce significant temperature responses of the regions near the positions of the nonlinear thermal loads, whereas the temperature responses of the remaining regions are very close to zero. According to the above physical characteristic, the original nonlinear problem can be transformed into some nonlinear problems with small structures and a linear problem with the original structure.

An efficient and accurate numerical method is presented for transient heat conduction problems with local nonlinearities, and some numerical examples demonstrate the high efficiency and accuracy of the proposed method.

The purpose of this paper is to study the effects of nonlinear partial slip on the walls for steady flow and heat transfer of an incompressible, thermodynamically compatible third grade fluid in a channel. The principal question the authors address in this paper is in regard to the applicability of the no‐slip condition at a solid‐liquid boundary. The authors present the effects of slip, magnetohydrodynamics (MHD) and heat transfer for the plane Couette, plane Poiseuille and plane Couette‐Poiseuille flows in a homogeneous and thermodynamically compatible third grade fluid. The problem of a non‐Newtonian plane Couette flow, fully developed plane Poiseuille flow and Couette‐Poiseuille flow are investigated.

The present investigation is an attempt to study the effects of nonlinear partial slip on the walls for steady flow and heat transfer of an incompressible, thermodynamically compatible third grade fluid in a channel. A very effective and higher order numerical scheme is used to solve the resulting system of nonlinear differential equations with nonlinear boundary conditions. Numerical solutions are obtained by solving nonlinear ordinary differential equations using Chebyshev spectral method.

Due to the nonlinear and highly complicated nature of the governing equations and boundary conditions, finding an analytical or numerical solution is not easy. The authors obtained numerical solutions of the coupled nonlinear ordinary differential equations with nonlinear boundary conditions using higher order Chebyshev spectral collocation method. Spectral methods are proven to offer a superior intrinsic accuracy for derivative calculations.

To the best of the authors' knowledge, no such analysis is available in the literature which can describe the heat transfer, MHD and slip effects simultaneously on the flows of the non‐Newtonian fluids.

The purpose of this paper is to present a new method based on the homotopy analysis method (HAM) with the aim of fast searching and calculating multiple solutions of nonlinear boundary value problems (NBVPs).

A major problem with the previously modified HAM, namely, predictor homotopy analysis method, which is used to predict multiplicity of solutions of NBVPs, is a time-consuming computation of high-order HAM-approximate solutions due to a symbolic variable namely “prescribed parameter”. The proposed new technique which is based on traditional shooting method, and the HAM cuts the dependency on the prescribed parameter.

To demonstrate the computational efficiency, the mentioned method is implemented on three important nonlinear exactly solvable differential equations, namely, the nonlinear MHD Jeffery–Hamel flow problem, the nonlinear boundary value problem arising in heat transfer and the strongly nonlinear Bratu problem.

The more high-order approximate solutions are computable, multiple solutions are easily searched and discovered and the more accurate solutions can be obtained depending on how nonhomogeneous boundary conditions are transcribed to the homogeneous boundary conditions.

The purpose of this paper is to determine both analytically and numerically the solution to a new one-dimensional equation for the propagation of small-amplitude waves in shallow waters that accounts for linear and nonlinear drift, diffusive attenuation, viscosity and dispersion, its dependence on the initial conditions, and its linear stability.

An implicit, finite difference method valid for both parabolic and second-order hyperbolic equations has been used to solve the equation in a truncated domain for five different initial conditions, a nil initial first-order time derivative and relaxation times linearly proportional to the viscosity coefficient.

A fast transition that depends on the coefficient of the linear drift, the diffusive attenuation and the power of the nonlinear drift are found for initial conditions corresponding to the exact solution of the generalized regularized long-wave equation. For initial Gaussian, rectangular and triangular conditions, the wave’s amplitude and speed increase as both the amplitude and the width of these conditions increase and decrease, respectively; wide initial conditions evolve into a narrow leading traveling wave of the pulse type and a train of slower oscillatory secondary ones. For the same initial mass and amplitude, rectangular initial conditions result in larger amplitude and velocity waves of the pulse type than Gaussian and triangular ones. The wave’s kinetic, potential and stretching energies undergo large changes in an initial layer whose thickness is on the order of the diffusive attenuation coefficient.

A new, one-dimensional equation for the propagation of small-amplitude waves in shallow waters is proposed and studied analytically and numerically. The equation may also be used to study the displacement of porous media subject to seismic effects, the dispersion of sound in tunnels, the attenuation of sound because of viscosity and/or heat and mass diffusion, the dynamics of second-order, viscoelastic fluids, etc., by appropriate choices of the parameters that appear in it.

In the last two decades with the rapid development of nonlinear science, there has appeared ever‐increasing interest of scientists and engineers in the analytical techniques for nonlinear problems. This paper considers linear and nonlinear systems that are not only regarded as general boundary value problems, but also are used as mathematical models in viscoelastic and inelastic flows. The purpose of this paper is to present the application of the homotopy‐perturbation method (HPM) and variational iteration method (VIM) to solve some boundary value problems in structural engineering and fluid mechanics.

Two new but powerful analytical methods, namely, He's VIM and HPM, are introduced to solve some boundary value problems in structural engineering and fluid mechanics.

Analytical solutions often fit under classical perturbation methods. However, as with other analytical techniques, certain limitations restrict the wide application of perturbation methods, most important of which is the dependence of these methods on the existence of a small parameter in the equation. Disappointingly, the majority of nonlinear problems have no small parameter at all. Furthermore, the approximate solutions solved by the perturbation methods are valid, in most cases, only for the small values of the parameters. In the present study, two powerful analytical methods HPM and VIM have been employed to solve the linear and nonlinear elastic beam deformation problems. The results reveal that these new methods are very effective and simple and do not require a large computer memory and can also be used for solving linear and nonlinear boundary value problems.

The results revealed that the VIM and HPM are remarkably effective for solving boundary value problems. These methods are very promoting methods which can be wildly utilized for solving mathematical and engineering problems.

Contact problems are among the most difficult ones in mechanics. Due to its practical importance, the problem has been receiving extensive research work over the years. The finite element method has been widely used to solve contact problems with various grades of complexity. Great progress has been made on both theoretical studies and engineering applications. This paper reviews some of the main developments in contact theories and finite element solution techniques for static contact problems. Classical and variational formulations of the problem are first given and then finite element solution techniques are reviewed. Available constraint methods, friction laws and contact searching algorithms are also briefly described. At the end of the paper, a bibliography is included, listing about seven hundred papers which are related to static contact problems and have been published in various journals and conference proceedings from 1976.

The numerical methods are of great importance for approximating the solutions of a system of nonlinear singular ordinary differential equations. In this paper, the authors present the biorthogonal flatlet multiwavelet collocation method (BFMCM) as a numerical scheme for a class of system of Lane–Emden equations with initial or boundary or four-point boundary conditions.

The approach is involved in combining the biorthogonal flatlet multiwavelet (BFM) with the collocation method. The authors investigate the properties and procedure of the BFMCM for first time on this class of equations. By using the BFM and the collocation points, the method is constructed and it transforms the nonlinear differential equations problem into a system of nonlinear algebraic equations. The unknown coefficients of the assuming solution are determined by solving the obtained system. Additionally, convergence analysis and numerical stability of the suggested method are provided.

According to the attained results, the proposed BFMCM has more accurate results in comparison with results of other methods. The maximum absolute errors are calculated by using the BFMCM for comparison purposes provided.

The key desirable properties of BFMCM are its efficiency, simple applicability and minimizes errors. Therefore, the proposed method can be used to solve nonlinear problems or problems with singular points.