Impulsive boundary value problems with nonlinear boundary conditions

In this paper, the upper and lower solution method and Schauder’s fixed point theorem are employed in the study of boundary value problems for a class of second-order impulsive ordinary differential equations with nonlinear boundary conditions. We prove the existence of solutions to the problem under the assumption that there exist lower and upper solutions associated with the problem.     

A meshless method based on boundary integral equations and radial basis functions for biharmonic-type problems

This paper presents a meshless method, which replaces the inhomogeneous biharmonic equation by two Poisson equations in terms of an intermediate function. The solution of the Poisson equation with the intermediate function as the right-hand term may be written as a sum of a particular solution and a homogeneous solution of a Laplace equation. The intermediate function is approximated by a series of radial basis functions. Then the particular solution is obtained via employing Kansa’s method, while the homogeneous solution is approximated by using the boundary radial point interpolation method by means of boundary integral equations. Besides, the proposed meshless method, in conjunction with the analog equation method, is further developed for solving generalized biharmonic-type problems. Some numerical tests illustrate the efficiency of the method proposed.     

Solution of Poisson’s equation by analytical boundary element integration

The solution of Poisson’s equation is essential for many branches of science and engineering such as fluid-mechanics, geosciences, and electrostatics. Solution of two-dimensional Poisson’s equations is carried out by BEM based on Galerkin Vector Method in which the integrals that appear in the boundary element method are expressed by analytical integration. In this paper, the Galerkin vector method is developed for more general case than presented in the previous works. The integrals are computed for constant and linear elements in BEM. By employing analytical integration in BEM computation, the numerical schemes and coordinate transformations can be avoided. The presented method can also be used for the multiple domain case. The results of the analytical integration are employed in BEM code and the obtained analytical expression will be applied to several examples where the exact solution exists. The produced results are in good agreement with the exact solution.     

First-order three-point boundary value problems at resonance

We consider three-point boundary value problems for a system of first-order equations in perturbed systems of ordinary differential equations at resonance. We obtain new results for the above boundary value problems with nonlinear boundary conditions. The existence of solutions is established by applying a version of Brouwer’s Fixed Point Theorem which is due to Miranda.     

Existence of solutions for anti-periodic boundary value problems

By using Schauder’s fixed point theorem and the lower and upper solutions method, we discuss the anti-periodic boundary value problem for a class of second-order differential equations. Some sufficient conditions for the existence of solutions are obtained.     

Solving two-point boundary value problems using combined homotopy perturbation method and Green’s function method

In this paper, an algorithm is proposed for the solution of second-order boundary value problems with two-point boundary conditions. The Green’s function method is applied first to transform the ordinary differential equation into an equivalent integral one, which has already satisfied the boundary conditions. And then, the homotopy perturbation method is used to the resulting equation to construct the numerical solution for such problems. Numerical examples demonstrate the efficiency and reliability of the algorithm developed, it is quite accurate and readily implemented for both linear and nonlinear differential equations with homogeneous and nonhomogeneous boundary conditions. Furthermore, the lower order approximation is of higher accuracy for most cases. Some other extended applications of this algorithm are also exhibited.     

Efficient long-time computations of time-domain boundary integrals for 2D and dissipative wave equation

Linear hyperbolic partial differential equations in a homogeneous medium, e.g., the wave equation describing the propagation and scattering of acoustic waves, can be reformulated as time-domain boundary integral equations. We propose an efficient implementation of a numerical discretization of such equations when the strong Huygens’ principle does not hold.For the numerical discretization, we make use of convolution quadrature in time and standard Galerkin boundary element method in space. The quadrature in time results in a discrete convolution of weights W_j with the boundary density evaluated at equally spaced time points. If the strong Huygens’ principle holds, W_j converge to 0 exponentially quickly for large enough j. If the strong Huygens’ principle does not hold, e.g., in even space dimensions or when some damping is present, the weights are never zero, thereby presenting a difficulty for efficient numerical computation.In this paper we prove that the kernels of the convolution weights approximate in a certain sense the time domain fundamental solution and that the same holds if both are differentiated in space. The tails of the fundamental solution being very smooth, this implies that the tails of the weights are smooth and can efficiently be interpolated. Further, we hint on the possibility to apply the fast and oblivious convolution quadrature algorithm of Schädle et al. to further reduce memory requirements for long-time computation. We discuss the efficient implementation of the whole numerical scheme and present numerical experiments.     

Periodic boundary value problems of second-order impulsive differential equations

This paper discusses periodic boundary value problems of second-order impulsive differential equations. By using Schaeffer’s fixed-point theorem and monotone iterative technique, some existence results are obtained.     

Approximate controllability of a class of semilinear degenerate systems with boundary control

This paper concerns a class of control systems governed by semilinear degenerate equations with boundary control in one-dimensional space. The control is proposed on the ‘degenerate’ part of the boundary. The control systems are shown to be approximately controllable by Kakutani's fixed point theorem.     

A new modification of the Adomian decomposition method for solving boundary value problems for higher order nonlinear differential equations

In this paper we propose a new modified recursion scheme for the resolution of multi-order and multi-point boundary value problems for nonlinear ordinary and partial differential equations by the Adomian decomposition method (ADM). Our new approach, including Duan’s convergence parameter, provides a significant computational advantage by allowing for the acceleration of convergence and expansion of the interval of convergence during calculations of the solution components for nonlinear boundary value problems, in particular for such cases when one of the boundary points lies outside the interval of convergence of the usual decomposition series. We utilize the boundary conditions to derive an integral equation before establishing the recursion scheme for the solution components. Thus we can derive a modified recursion scheme without any undetermined coefficients when computing successive solution components, whereas several prior recursion schemes have done so. This modification also avoids solving a sequence of nonlinear algebraic equations for the undetermined coefficients fraught with multiple roots, which is required to complete calculation of the solution by several prior modified recursion schemes using the ADM.     

A preconditioning strategy for boundary element Galerkin methods

The Dirichlet and Neumann problems for the Laplacian are reformulated in the usual way as boundary integral equations of the first kind with symmetric kernels. These integral equations are solved using Galerkin's method with piecewise-constant and piecewise-linear boundary elements, respectively. In both cases, the stiffness matrix is symmetric and positive-definite, and has a condition number of order N, the number of degrees of freedom. By contrast, the condition number of the product of the two stiffness matrices is bounded independently of N. Hence, we can use the Neumann stiffness matrix to precondition the Dirichlet stiffness matrix, and vice versa. © 1997 John Wiley & Sons, Inc.     

Reconstruction of closely spaced small inhomogeneities via boundary measurements for the full time-dependent Maxwell’s equations

We consider for the full time-dependent Maxwell’s equations the inverse problem of identifying locations and certain properties of small electromagnetic inhomogeneities in a homogeneous background medium from dynamic boundary measurements on the boundary for a finite time interval.     

Dual reciprocity hybrid radial boundary node method for the analysis of Kirchhoff plates

A meshless method of dual reciprocity hybrid radial boundary node method (DHRBNM) for the analysis of arbitrary Kirchhoff plates is presented, which combines the advantageous properties of meshless method, radial point interpolation method (RPIM) and BEM. The solution in present method comprises two parts, i.e., the complementary solution and the particular solution. The complementary solution is solved by hybrid radial boundary node method (HRBNM), in which a three-field interpolation scheme is employed, and the boundary variables are approximated by RPIM, which is applied instead of moving least square (MLS) and obtains the Kronecker’s delta property where the traditional HBNM does not satisfy. The internal variables are interpolated by two groups of symmetric fundamental solutions. Based on those, a hybrid displacement variational principle for Kirchhoff plates is developed, and a meshless method of HRBNM for solving biharmonic problems is obtained, by which the complementary solution can be solved.     

Special Green’s function boundary element approach for steady-state axisymmetric heat conduction across low and high conducting planar interfaces

The problem of determining the steady-state axisymmetric temperature distribution in a bimaterial with a planar interface is considered here. The interface is either low or high conducting. Special Green’s functions satisfying the thermal conditions on the interface are derived and employed to obtain boundary integral equations whose path of integration does not include the interface. Boundary element procedures that do not require the interface to be discretized into elements are proposed for solving the problem under consideration.     

Series solutions of nano boundary layer flows by means of the homotopy analysis method

We present here a ‘similar’ solution for the nano boundary layer with nonlinear Navier boundary condition. Three types of flows are considered: (i) the flow past a wedge; (ii) the flow in a convergent channel; (iii) the flow driven by an exponentially-varying outer flows. The resulting differential equations are solved by the homotopy analysis method. Different from the perturbation methods, the present method is independent of small physical parameters so that it is applicable for not only weak but also strong nonlinear flow phenomena. Numerical results are compared with the available exact results to demonstrate the validity of the present solution. The effects of the slip length ?, the index parameters n and m on the velocity profile and the tangential stress are investigated and discussed.     

Discontinuous Legendre wavelet element method for elliptic partial differential equations

By incorporating the Legendre multiwavelet into the discontinuous Galerkin (DG) method, this paper presents a novel approach for solving Poisson’s equation with Dirichlet boundary, which is known as the discontinuous Legendre multiwavelet element (DLWE) method, derive an adaptive algorithm for the method, and estimate the approximating error of its numerical fluxes. One striking advantage of our method is that the differential operator, boundary conditions and numerical fluxes involved in the elementwise computation can be done with lower time cost. Numerical experiments demonstrate the validity of this method. Furthermore, this paper generalizes the DLWE method to the general elliptic equations defined on a bounded domain and describes the possibilities of constructing optimal adaptive algorithm. The proposed method and its generalizations are also applicable to some other kinds of partial differential equations.     

Exact solutions of a class of second-order nonlocal boundary value problems and applications

In this paper, solutions of a class of second-order differential equations with some multi-point boundary conditions are studied. We give exact expressions of the solutions for the linear m-point boundary problems by the Green’s functions. As applications, we study uniqueness and iteration of the positive solutions for a nonlinear singular second-order m-point boundary value problem.     

Positive solutions of a nonlinear impulsive equation with periodic boundary conditions

In this paper, we investigate nonlinear second order differential equations subject to linear impulse conditions and periodic boundary conditions. Sign properties of an associated Green’s function are exploited to get existence results for positive solutions of the nonlinear boundary value problem with impulse. Upper and lower bounds for positive solutions are also given. The results obtained yield periodic positive solutions of the corresponding periodic impulsive nonlinear differential equation on the whole real axis.     

Approximation of boundary values in the boundary element method

Various techniques may be applied to the approximation of the unknown boundary functions involved in the boundary element method (BEM). Several techniques have been examined numerically to find the most efficient. Techniques considered were: Lagrangian polynomials of the zeroth, first and second orders; spline functions; and the novel weighted minimization technique used successfully in the finite difference method (FDM) for arbitrarily irregular meshes. All these approaches have been used in the BEM for the numerical analysis of plates with various boundary conditions.Both coarse and fine grids on the boundary have been assumed. Maximal errors of the deflections of each plate and the bending moments have been found and the effective computer CPU times determined.Analysis of the results showed that, for the same computer time, the greatest accuracy was obtained by the weighted FDM approach. In the case of the Lagrange approximation, higher order polynomials have proved more efficient. The spline technique yielded more accurate results, but with a higher CPU time.Two discretization approaches have been investigated: the least-squares technique and the collocation method. Despite the fact that the simultaneous algebraic equations obtained were not symmetric, the collocation approach has been confirmed as clearly superior to the least-squares technique, because of the amount of computation time used.