Explore the vast range of titles available.

For the splitting iterative scheme to solve the system of linear equations, an equivalent form in terms of descent and residual vectors is formulated. We propose an extrapolation technique using the new formulation, such that a new splitting iterative scheme (NSIS) can be simply generated from the original one by inserting an acceleration parameter preceding the descent vector. The spectral radius of the NSIS is proven to be smaller than the original one, and so has a faster convergence speed. The orthogonality of consecutive residual vectors is coined into the second NSIS, from which a stepwise varying orthogonalization factor can be derived explicitly. Multiplying the descent vector by the factor, the second NSIS is proven to be absolutely convergent. The modification is based on the maximal reduction of residual vector norm. Two-parameter and three-parameter NSIS are investigated, wherein the optimal value of one parameter is obtained by using the maximization technique. The splitting iterative schemes are unified to have the same iterative form, but endowed with different governing equations for the descent vector. Some examples are examined to exhibit the performance of the proposed extrapolation techniques used in the NSIS.

For nonlinear third-order three-point boundary value problems (BVPs), we develop two algorithms to find solutions, which automatically satisfy the specified three-point boundary conditions. We construct a boundary shape function (BSF), which is designed to automatically satisfy the boundary conditions and can be employed to develop new algorithms by assigning two different roles of free function in the BSF. In the first algorithm, we let the free functions be complete functions and the BSFs be the new bases of the solution, which not only satisfy the boundary conditions automatically, but also can be used to find solution by a collocation technique. In the second algorithm, we let the BSF be the solution of the BVP and the free function be another new variable, such that we can transform the BVP to a corresponding initial value problem for the new variable, whose initial conditions are given arbitrarily and terminal values are determined by iterations; hence, we can quickly find very accurate solution of nonlinear third-order three-point BVP through a few iterations. Numerical examples confirm the performance of the new algorithms.

Besides the closed-form expansion coefficients of a weak-form numerical differentiator (WFND), we introduce a cubic boundary shape function with the aid of two parameters for reducing the boundary errors of fourth-order numerical derivatives to zero. So that the accuracy of numerical derivatives obtained by the new WFND can be improved significantly. The fourth-order numerical derivation can be modeled as a linear beam equation subjecting to specified boundary conditions and displacements to recover an unknown forcing term. By means of boundary shape functions, two numerical collocation methods automatically satisfying the boundary conditions are developed. For a simply supported linear Euler–Bernoulli beam with an elastic foundation, the unknown spatially–temporally dependent force is retrieved. The displacement at a final time and strain on the right-boundary of the beam are over-specified to recover the external force using the method of superposition of boundary shape functions (MSBSF). When the displacement is determined to satisfy the prescribed right-boundary strain, we can recover an unknown spatially–temporally dependent force by inserting the displacement into the linear beam equation. An embedded method (EM) is developed to transform the linear beam model into a vibrating linear beam equation, and then we can develop a robust technique to compute the fourth-order derivative of noisy data by using the EM and MSBSF. The four proposed methods for evaluating the fourth-order derivatives of noisy data are efficient and accurate.

This paper introduces a singular distance function rs in terms of a symmetric non-negative metric tensor S. If S satisfies a quadratic matrix equation involving a parameter β, then for the Laplace equation rsβ is a non-singular generalized radial basis function solution if 2 > β > 0, and a weaker singularity fundamental solution if −1 < β < 0. With a unit vector as a medium to express S, we can derive the metric tensor in closed form and prove that S is a singular projection operator. For the anisotropic Laplace equation, the corresponding closed-form representation of S is also derived. The concept of non-singular generalized radial basis function solution for the Laplace-type equations is novel and useful, which has not yet appeared in the literature. In addition, a logarithmic type method of fundamental solutions is developed for the anisotropic Laplace equation. Owing to non-singularity and weaker singularity of the bases of solutions, numerical experiments verify the accuracy and efficiency of the proposed methods.

The Navier equations are reformulated to be third-order partial differential equations. New anti-Cauchy-Riemann equations can express a general solution in 2D space for incompressible materials. Based on the third-order solutions in 3D space and the Boussinesq–Galerkin method, a third-order method of fundamental solutions (MFS) is developed. For the 3D Navier equation in linear elasticity, we present three new general solutions, which have appeared in the literature for the first time, to signify the theoretical contributions of the present paper. The first one is in terms of a biharmonic function and a harmonic function. The completeness of the proposed general solution is proven by using the solvability conditions of the equations obtained by equating the proposed general solution to the Boussinesq–Galerkin solution. The second general solution is expressed in terms of a harmonic vector, which is simpler than the Slobodianskii general solution, and the traditional MFS. The main achievement is that the general solution is complete, and the number of harmonic functions, three, is minimal. The third general solution is presented by a harmonic vector and a biharmonic vector, which are subjected to a constraint equation. We derive a specific solution by setting the two vectors in the third general solution as the vectorizations of a single harmonic potential. Hence, we have a simple approach to the Slobodianskii general solution. The applications of the new solutions are demonstrated. Owing to the minimality of the harmonic functions, the resulting bases generated from the new general solution are complete and linearly independent. Numerical instability can be avoided by using the new bases. To explore the efficiency and accuracy of the proposed MFS variant methods, some examples are tested.

The development of simple and yet accurate formulations of frequency–amplitude relationships for non-conservative nonlinear oscillators is an important issue. The present paper is concerned with integral-type frequency–amplitude formulas in the dimensionless time domain and time domain to accurately determine vibrational frequencies of nonlinear oscillators. The novel formulation is a balance of kinetic energy and the work during motion of the nonlinear oscillator within one period; its generalized formulation permits a weight function to appear in the integral formula. The exact values of frequencies can be obtained when exact solutions are inserted into the formulas. In general, the exact solution is not available; hence, low-order periodic functions as trial solutions are inserted into the formulas to obtain approximate values of true frequencies. For conservative nonlinear oscillators, a powerful technique is developed in terms of a weighted integral formula in the spatial domain, which is directly derived from the governing ordinary differential equation (ODE) multiplied by a weight function, and integrating the resulting equation after inserting a general trial ODE to acquire accurate frequency. The free parameter is involved in the frequency–amplitude formula, whose optimal value is achieved by minimizing the absolute error to fulfill the periodicity conditions. Several examples involving two typical non-conservative nonlinear oscillators are explored to display the effectiveness and accuracy of the proposed integral-type formulations.

For a non-conservative nonlinear oscillator (NCNO) having a periodic solution, the existence of the first integral is a certain symmetry of the nonlinear dynamical system, which signifies the balance of kinetic energy and potential energy. A first-order nonlinear ordinary differential equation (ODE) is used to derive the first integral, which, equipped with a right-end boundary condition, can determine an implicit potential function for computing the period by an exact integral formula. However, the integrand is singular, which renders a less accurate value of the period. A generalized integral conservation law endowed with a weight function is constructed, which is proved to be equivalent to the exact integral formula. Minimizing the error to satisfy the periodicity conditions, the optimal initial value of the weight function is determined. Two non-iterative methods are developed by integrating three first-order ODEs or two first-order ODEs to compute the period. Very accurate value of the period can be observed upon testing five examples. For the NCNO without having the first integral, the integral-type period formula is derived. Four examples belong to the Liénard equation, involving the van der Pol equation, are evaluated by the proposed iterative method to determine the oscillatory amplitude and period. For the case with one or more limit cycles, the amplitude and period can be estimated very accurately. For the NCNO of a broad type with or without having the first integral, the present paper features a solid theoretical foundation and contributes integral-type formulations for the determination of the oscillatory period. The development of new numerical algorithms and extensive validation across a diverse set of examples is given.

It is difficult to exactly and automatically satisfy nonseparable multipoint boundary conditions by numerical methods. With this in mind, we develop a novel algorithm to find solution for a second-order nonlinear boundary value problem (BVP), which automatically satisfies the multipoint boundary conditions prescribed. A novel concept of boundary shape function (BSF) is introduced, whose existence is proven, and it can satisfy the multipoint boundary conditions a priori. In the BSF, there exists a free function, from which we can develop an iterative algorithm by letting the BSF be the solution of the BVP and the free function be another variable. Hence, the multipoint nonlinear BVP is properly transformed to an initial value problem for the new variable, whose initial conditions are given arbitrarily. The BSF method (BSFM) can find very accurate solution through a few iterations.

A new concept of projective solution is introduced for the multi-dimensional Laplace equations. We project the field point onto a characteristic vector to obtain a projective variable, which can be used to reduce the Laplace equations to a second-order ordinary differential equation with only a leading term multiplied by the squared norm of the characteristic vector. The projective solutions involve characteristic vectors as parameters, which must be complex numbers to satisfy a null equation. Since the projective variable is a complex variable, we can construct the analytic function based on the conventional complex analytic function theory. Both the analytic function and the Cauchy–Riemann equations are generalized for the multi-dimensional Laplace equations. A powerful numerical technique to solve the 3D Laplace equation with high accuracy is available by further developing the Trefftz-type bases. Numerical experiments confirm the accuracy and efficiency of the projective solutions method (PSM).

A new frequency–amplitude formula by improving an ancient Chinese mathematics method results in a modification of He’s formula. The Chinese mathematics method that expresses via a fixed-point Newton form is proven to be equivalent to the original nonlinear frequency equation. We modify the fixed-point Newton method by adding a term in the denominator, and then a new frequency–amplitude formula including a parameter is derived. Upon using the new frequency formula with the parameter by minimizing the absolute error of the periodicity condition, one can significantly raise the accuracy of the frequency several orders. The innovative idea of a linearly perturbed frequency equation is a simple extension of the original frequency equation, which is supplemented by a linear term to acquire a highly precise frequency for the nonlinear oscillators. In terms of a differentiable weight function, an integral-type formula is coined to expeditiously estimate the frequency; it is a generalized conservation law for the damped nonlinear oscillator. To seek second-order periodic solutions of nonlinear oscillators, a linearized residual Galerkin method (LRGM) is developed whose process to find the second-order periodic solution and the vibrational frequency is quite simple. A hybrid method is achieved through a combination of the linearly perturbed frequency equation and the LRGM; very accurate frequency and second-order periodic solutions can be obtained. Examples reveal high efficacy and accuracy of the proposed methods; the mathematical reliability of these methods is clarified.