A high-precision and space-time fully decoupled numerical method is developed for a class of nonlinear initial boundary value problems. It is established based on a proposed Coiflet-based approximation scheme with an ...A high-precision and space-time fully decoupled numerical method is developed for a class of nonlinear initial boundary value problems. It is established based on a proposed Coiflet-based approximation scheme with an adjustable high order for the functions over a bounded interval, which allows the expansion coefficients to be explicitly expressed by the function values at a series of single points. When the solution method is used, the nonlinear initial boundary value problems are first spatially discretized into a series of nonlinear initial value problems by combining the proposed wavelet approximation and the conventional Galerkin method, and a novel high-order step-by-step time integrating approach is then developed for the resulting nonlinear initial value problems with the same function approximation scheme based on the wavelet theory. The solution method is shown to have the N th-order accuracy, as long as the Coiflet with [0, 3 N-1]compact support is adopted, where N can be any positive even number. Typical examples in mechanics are considered to justify the accuracy and efficiency of the method.展开更多
This paper proposes a new collocation method for initial value problems of second order ODEs based on the Laguerre-Gauss interpolation. It provides the global numerical solutions and possesses the spectral accuracy. N...This paper proposes a new collocation method for initial value problems of second order ODEs based on the Laguerre-Gauss interpolation. It provides the global numerical solutions and possesses the spectral accuracy. Numerical results demonstrate its high efficiency.展开更多
In this paper, the standard homotopy analysis method was applied to initial value problems of the second order with some types of discontinuities, for both linear and nonlinear cases. To show the high accuracy of the ...In this paper, the standard homotopy analysis method was applied to initial value problems of the second order with some types of discontinuities, for both linear and nonlinear cases. To show the high accuracy of the solution results compared with the exact solution, a comparison of the numerical results was made applying the standard homotopy analysis method with the iteration of the integral equation and the numerical solution with the Simpson rule. Also, the maximum absolute error, , the maximum relative error, the maximum residual error and the estimated order of convergence were given. The research is meaningful and I recommend it to be published in the journal.展开更多
This paper mainly presents Euler method and fourth-order Runge Kutta Method (RK4) for solving initial value problems (IVP) for ordinary differential equations (ODE). The two proposed methods are quite efficient and pr...This paper mainly presents Euler method and fourth-order Runge Kutta Method (RK4) for solving initial value problems (IVP) for ordinary differential equations (ODE). The two proposed methods are quite efficient and practically well suited for solving these problems. In order to verify the ac-curacy, we compare numerical solutions with the exact solutions. The numerical solutions are in good agreement with the exact solutions. Numerical comparisons between Euler method and Runge Kutta method have been presented. Also we compare the performance and the computational effort of such methods. In order to achieve higher accuracy in the solution, the step size needs to be very small. Finally we investigate and compute the errors of the two proposed methods for different step sizes to examine superiority. Several numerical examples are given to demonstrate the reliability and efficiency.展开更多
In this paper, some modifications of Adomian decomposition method are presented for solving initial value problems in ordinary differential equations. Also, the restarted and two-step methods are applied to the proble...In this paper, some modifications of Adomian decomposition method are presented for solving initial value problems in ordinary differential equations. Also, the restarted and two-step methods are applied to the problem. The effectiveness of the each modified is verified by several examples.展开更多
This paper is concerned with Godunov＇s scheme for the initial-boundary value problem of scalar conservation laws. A kind of Godunov＇s scheme, which satisfies the boundary entropy condition, was given. By use of the ...This paper is concerned with Godunov＇s scheme for the initial-boundary value problem of scalar conservation laws. A kind of Godunov＇s scheme, which satisfies the boundary entropy condition, was given. By use of the scheme, numerical simulation for the weak entropy solution to the initial-boundary value problem of scalar conservation laws is conducted.展开更多
A time-spectral method for solution of initial value partial differential equations is outlined. Multivariate Chebyshev series are used to represent all temporal, spatial and physical parameter domains in this general...A time-spectral method for solution of initial value partial differential equations is outlined. Multivariate Chebyshev series are used to represent all temporal, spatial and physical parameter domains in this generalized weighted residual method (GWRM). The approximate solutions obtained are thus analytical, finite order multivariate polynomials. The method avoids time step limitations. To determine the spectral coefficients, a system of algebraic equations is solved iteratively. A root solver, with excellent global convergence properties, has been developed. Accuracy and efficiency are controlled by the number of included Chebyshev modes and by use of temporal and spatial subdomains. As examples of advanced application, stability problems within ideal and resistive magnetohydrodynamics (MHD) are solved. To introduce the method, solutions to a stiff ordinary differential equation are demonstrated and discussed. Subsequently, the GWRM is applied to the Burger and forced wave equations. Comparisons with the explicit Lax-Wendroff and implicit Crank-Nicolson finite difference methods show that the method is accurate and efficient. Thus the method shows potential for advanced initial value problems in fluid mechanics and MHD.展开更多
In this article,we consider to solve the inverse initial value problem for an inhomogeneous space-time fractional diffusion equation.This problem is ill-posed and the quasi-boundary value method is proposed to deal wi...In this article,we consider to solve the inverse initial value problem for an inhomogeneous space-time fractional diffusion equation.This problem is ill-posed and the quasi-boundary value method is proposed to deal with this inverse problem and obtain the series expression of the regularized solution for the inverse initial value problem.We prove the error estimates between the regularization solution and the exact solution by using an a priori regularization parameter and an a posteriori regularization parameter choice rule.Some numerical results in one-dimensional case and two-dimensional case show that our method is efficient and stable.展开更多
This research paper deals with the boundary and initial value problems for the Bratu-type model by using the New Improved Variational Homotopy Perturbation Method. The New Method does not require discritization, linea...This research paper deals with the boundary and initial value problems for the Bratu-type model by using the New Improved Variational Homotopy Perturbation Method. The New Method does not require discritization, linearization or any restrictive assumption of any form in providing analytical or approximate solutions to linear and nonlinear equation without the integral related with nonlinear term. Theses virtues make it to be reliable and its efficiency is demonstrated with numerical examples.展开更多
This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditiona...This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.展开更多
A uniform high-order method is. presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order pro...A uniform high-order method is. presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order problems ( 1.4), i.e., the solution of (1.1) is a linear combination of the solutions of (1.4). Then we derive a uniformly O (hm+1) accurate scheme for the first-order problems (1.4), where m is an arbitrary nonnegative integer, so we can get a uniformly O (hm+1) accurate solution of the original problem (1.1) by relation (1.3). Some illustrative numerical results are also given.展开更多
We employ the Duan-Rach-Wazwaz modified Adomian decomposition method for solving initial value problems for the systems of nonlinear ordinary differential equations numerically. In order to confirm practicality, robus...We employ the Duan-Rach-Wazwaz modified Adomian decomposition method for solving initial value problems for the systems of nonlinear ordinary differential equations numerically. In order to confirm practicality, robustness and reliability of the method, we compare the results from the modified Adomian decomposition method with those from the MATHEMATICA solutions and also from the fourth-order Runge Kutta method solutions in some cases. Furthermore, we apply Padé approximants technique to improve the solutions of the modified decomposition method whenever the exact solutions exist.展开更多
In this paper,we apply Fokas unified method to study the initial boundary value(IBV)problems for nonlinear integrable equation with 3×3 Lax pair on the finite interval[0,L].The solution can be expressed by the so...In this paper,we apply Fokas unified method to study the initial boundary value(IBV)problems for nonlinear integrable equation with 3×3 Lax pair on the finite interval[0,L].The solution can be expressed by the solution of a 3×3 Riemann-Hilbert(RH)problem.The relevant jump matrices are written in terms of matrix-value spectral functions s(k),S(k),S_(l)(k),which are determined by initial data at t=0,boundary values at x=0 and boundary values at x=L,respectively.What's more,since the eigenvalues of 3×3 coefficient matrix of k spectral parameter in Lax pair are three different values,search for the path of analytic functions in RH problem becomes a very interesting thing.展开更多
In this paper, a new Fourier-differential transform method (FDTM) based on differential transformation method (DTM) is proposed. The method can effectively and quickly solve linear and nonlinear partial differential e...In this paper, a new Fourier-differential transform method (FDTM) based on differential transformation method (DTM) is proposed. The method can effectively and quickly solve linear and nonlinear partial differential equations with initial boundary value (IBVP). According to boundary condition, the initial condition is expanded into a Fourier series. After that, the IBVP is transformed to an iterative relation in K-domain. The series solution or exact solution can be obtained. The rationality and practicability of the algorithm FDTM are verified by comparisons of the results obtained by FDTM and the existing analytical solutions.展开更多
The multisplitting algorithm for solving large systems of ordinary differential equations on parallel computers was introduced by Jeltsch and Pohl in [1]. On fixed time intervals conver gence results could be derived ...The multisplitting algorithm for solving large systems of ordinary differential equations on parallel computers was introduced by Jeltsch and Pohl in [1]. On fixed time intervals conver gence results could be derived if the subsystems are solving exactly.Firstly,in theis paper,we deal with an extension of the waveform relaxation algorithm by us ing multisplittin AOR method based on an overlapping block decomposition. We restricted our selves to equidistant timepoints and dealed with the case that an implicit integration method was used to solve the subsystems numerically in parallel. Then we have proved convergence of multi splitting AOR waveform relaxation algorithm on a fixed window containing a finite number of timepoints.展开更多
In this paper, the random Euler and random Runge-Kutta of the second order methods are used in solving random differential initial value problems of first order. The conditions of the mean square convergence of the nu...In this paper, the random Euler and random Runge-Kutta of the second order methods are used in solving random differential initial value problems of first order. The conditions of the mean square convergence of the numerical solutions are studied. The statistical properties of the numerical solutions are computed through numerical case studies.展开更多
In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor depos...In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor deposition) processes, means the flow of species to a gas-phase, which are influenced by an electric field. Such a field we can model by wave equations. The main contributions are to improve the standard discretization schemes of each part of the coupling equation. So we discuss an improvement with implicit Runge- Kutta methods instead of the Yee’s algorithm. Further we balance the solver method between the Maxwell and Transport equation.展开更多
We present the optimal homotopy asymptotic method (OHAM) to find the numerical solution of the second order initial value problems of Bratu-type. We solve some examples to illustrate the validity and efficiency of the...We present the optimal homotopy asymptotic method (OHAM) to find the numerical solution of the second order initial value problems of Bratu-type. We solve some examples to illustrate the validity and efficiency of the method.展开更多
Fourth order differential equations are considered to develop the class of methods for the numerical solution of boundary value problems. In this paper, we have discussed the regions of absolute stability of fourth or...Fourth order differential equations are considered to develop the class of methods for the numerical solution of boundary value problems. In this paper, we have discussed the regions of absolute stability of fourth order boundary value problems. Methods proposed and derived in this paper are applied to solve a fourth-order boundary value problem. Numerical results are given to illustrate the efficiency of our methods and compared with exact solution.展开更多
In this paper, we extend the reliable modification of the Adomian Decom-position Method coupled to the Lesnic’s approach to solve boundary value problems and initial boundary value problems with mixed boundary condit...In this paper, we extend the reliable modification of the Adomian Decom-position Method coupled to the Lesnic’s approach to solve boundary value problems and initial boundary value problems with mixed boundary conditions for linear and nonlinear partial differential equations. The method is applied to different forms of heat and wave equations as illustrative examples to exhibit the effectiveness of the method. The method provides the solution in a rapidly convergent series with components that can be computed iteratively. The numerical results for the illustrative examples obtained show remarkable agreement with the exact solutions. We also provide some graphical representations for clear-cut comparisons between the solutions using Maple software.展开更多
基金Project supported by the National Natural Science Foundation of China(No.11472119)the Fundamental Research Funds for the Central Universities(No.lzujbky-2017-ot11)the 111 Project(No.B14044)
文摘A high-precision and space-time fully decoupled numerical method is developed for a class of nonlinear initial boundary value problems. It is established based on a proposed Coiflet-based approximation scheme with an adjustable high order for the functions over a bounded interval, which allows the expansion coefficients to be explicitly expressed by the function values at a series of single points. When the solution method is used, the nonlinear initial boundary value problems are first spatially discretized into a series of nonlinear initial value problems by combining the proposed wavelet approximation and the conventional Galerkin method, and a novel high-order step-by-step time integrating approach is then developed for the resulting nonlinear initial value problems with the same function approximation scheme based on the wavelet theory. The solution method is shown to have the N th-order accuracy, as long as the Coiflet with [0, 3 N-1]compact support is adopted, where N can be any positive even number. Typical examples in mechanics are considered to justify the accuracy and efficiency of the method.
基金supported by the National Natural Science Foundation of China(No.11171227)the Ph.D.Programs Foundation of Ministry of Education of China(No.20080270001)+2 种基金the Shanghai Leading Academic Discipline Project(No.S30405)the Fund for E-Institute of Shanghai Universities(No.E03004)the Foundation for Distinguished Young Talents in Higher Education of Guangdong of China(No.LYM09138)
文摘This paper proposes a new collocation method for initial value problems of second order ODEs based on the Laguerre-Gauss interpolation. It provides the global numerical solutions and possesses the spectral accuracy. Numerical results demonstrate its high efficiency.
文摘In this paper, the standard homotopy analysis method was applied to initial value problems of the second order with some types of discontinuities, for both linear and nonlinear cases. To show the high accuracy of the solution results compared with the exact solution, a comparison of the numerical results was made applying the standard homotopy analysis method with the iteration of the integral equation and the numerical solution with the Simpson rule. Also, the maximum absolute error, , the maximum relative error, the maximum residual error and the estimated order of convergence were given. The research is meaningful and I recommend it to be published in the journal.
文摘This paper mainly presents Euler method and fourth-order Runge Kutta Method (RK4) for solving initial value problems (IVP) for ordinary differential equations (ODE). The two proposed methods are quite efficient and practically well suited for solving these problems. In order to verify the ac-curacy, we compare numerical solutions with the exact solutions. The numerical solutions are in good agreement with the exact solutions. Numerical comparisons between Euler method and Runge Kutta method have been presented. Also we compare the performance and the computational effort of such methods. In order to achieve higher accuracy in the solution, the step size needs to be very small. Finally we investigate and compute the errors of the two proposed methods for different step sizes to examine superiority. Several numerical examples are given to demonstrate the reliability and efficiency.
文摘In this paper, some modifications of Adomian decomposition method are presented for solving initial value problems in ordinary differential equations. Also, the restarted and two-step methods are applied to the problem. The effectiveness of the each modified is verified by several examples.
基金Project supported by the National Natural Science Foundation of China （Grant No. 10671120）
文摘This paper is concerned with Godunov＇s scheme for the initial-boundary value problem of scalar conservation laws. A kind of Godunov＇s scheme, which satisfies the boundary entropy condition, was given. By use of the scheme, numerical simulation for the weak entropy solution to the initial-boundary value problem of scalar conservation laws is conducted.
文摘A time-spectral method for solution of initial value partial differential equations is outlined. Multivariate Chebyshev series are used to represent all temporal, spatial and physical parameter domains in this generalized weighted residual method (GWRM). The approximate solutions obtained are thus analytical, finite order multivariate polynomials. The method avoids time step limitations. To determine the spectral coefficients, a system of algebraic equations is solved iteratively. A root solver, with excellent global convergence properties, has been developed. Accuracy and efficiency are controlled by the number of included Chebyshev modes and by use of temporal and spatial subdomains. As examples of advanced application, stability problems within ideal and resistive magnetohydrodynamics (MHD) are solved. To introduce the method, solutions to a stiff ordinary differential equation are demonstrated and discussed. Subsequently, the GWRM is applied to the Burger and forced wave equations. Comparisons with the explicit Lax-Wendroff and implicit Crank-Nicolson finite difference methods show that the method is accurate and efficient. Thus the method shows potential for advanced initial value problems in fluid mechanics and MHD.
基金The project is supported by the National Natural Science Foundation of China(11561045,11961044)the Doctor Fund of Lan Zhou University of Technology.
文摘In this article,we consider to solve the inverse initial value problem for an inhomogeneous space-time fractional diffusion equation.This problem is ill-posed and the quasi-boundary value method is proposed to deal with this inverse problem and obtain the series expression of the regularized solution for the inverse initial value problem.We prove the error estimates between the regularization solution and the exact solution by using an a priori regularization parameter and an a posteriori regularization parameter choice rule.Some numerical results in one-dimensional case and two-dimensional case show that our method is efficient and stable.
文摘This research paper deals with the boundary and initial value problems for the Bratu-type model by using the New Improved Variational Homotopy Perturbation Method. The New Method does not require discritization, linearization or any restrictive assumption of any form in providing analytical or approximate solutions to linear and nonlinear equation without the integral related with nonlinear term. Theses virtues make it to be reliable and its efficiency is demonstrated with numerical examples.
基金The project supported by the National Key Basic Research and Development Foundation of the Ministry of Science and Technology of China （G2000048702, 2003CB716707）the National Science Fund for Distinguished Young Scholars （10025208）+1 种基金 the National Natural Science Foundation of China （Key Program） （10532040） the Research Fund for 0versea Chinese （10228028）.
文摘This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.
文摘A uniform high-order method is. presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order problems ( 1.4), i.e., the solution of (1.1) is a linear combination of the solutions of (1.4). Then we derive a uniformly O (hm+1) accurate scheme for the first-order problems (1.4), where m is an arbitrary nonnegative integer, so we can get a uniformly O (hm+1) accurate solution of the original problem (1.1) by relation (1.3). Some illustrative numerical results are also given.
文摘We employ the Duan-Rach-Wazwaz modified Adomian decomposition method for solving initial value problems for the systems of nonlinear ordinary differential equations numerically. In order to confirm practicality, robustness and reliability of the method, we compare the results from the modified Adomian decomposition method with those from the MATHEMATICA solutions and also from the fourth-order Runge Kutta method solutions in some cases. Furthermore, we apply Padé approximants technique to improve the solutions of the modified decomposition method whenever the exact solutions exist.
基金supported by the National Natural Science Foundation of China(11901167,11971313 and 51879045)Key scientific research projects of higher education institutions in Henan,China(18B110008).
文摘In this paper,we apply Fokas unified method to study the initial boundary value(IBV)problems for nonlinear integrable equation with 3×3 Lax pair on the finite interval[0,L].The solution can be expressed by the solution of a 3×3 Riemann-Hilbert(RH)problem.The relevant jump matrices are written in terms of matrix-value spectral functions s(k),S(k),S_(l)(k),which are determined by initial data at t=0,boundary values at x=0 and boundary values at x=L,respectively.What's more,since the eigenvalues of 3×3 coefficient matrix of k spectral parameter in Lax pair are three different values,search for the path of analytic functions in RH problem becomes a very interesting thing.
文摘In this paper, a new Fourier-differential transform method (FDTM) based on differential transformation method (DTM) is proposed. The method can effectively and quickly solve linear and nonlinear partial differential equations with initial boundary value (IBVP). According to boundary condition, the initial condition is expanded into a Fourier series. After that, the IBVP is transformed to an iterative relation in K-domain. The series solution or exact solution can be obtained. The rationality and practicability of the algorithm FDTM are verified by comparisons of the results obtained by FDTM and the existing analytical solutions.
文摘The multisplitting algorithm for solving large systems of ordinary differential equations on parallel computers was introduced by Jeltsch and Pohl in [1]. On fixed time intervals conver gence results could be derived if the subsystems are solving exactly.Firstly,in theis paper,we deal with an extension of the waveform relaxation algorithm by us ing multisplittin AOR method based on an overlapping block decomposition. We restricted our selves to equidistant timepoints and dealed with the case that an implicit integration method was used to solve the subsystems numerically in parallel. Then we have proved convergence of multi splitting AOR waveform relaxation algorithm on a fixed window containing a finite number of timepoints.
文摘In this paper, the random Euler and random Runge-Kutta of the second order methods are used in solving random differential initial value problems of first order. The conditions of the mean square convergence of the numerical solutions are studied. The statistical properties of the numerical solutions are computed through numerical case studies.
文摘In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor deposition) processes, means the flow of species to a gas-phase, which are influenced by an electric field. Such a field we can model by wave equations. The main contributions are to improve the standard discretization schemes of each part of the coupling equation. So we discuss an improvement with implicit Runge- Kutta methods instead of the Yee’s algorithm. Further we balance the solver method between the Maxwell and Transport equation.
文摘We present the optimal homotopy asymptotic method (OHAM) to find the numerical solution of the second order initial value problems of Bratu-type. We solve some examples to illustrate the validity and efficiency of the method.
文摘Fourth order differential equations are considered to develop the class of methods for the numerical solution of boundary value problems. In this paper, we have discussed the regions of absolute stability of fourth order boundary value problems. Methods proposed and derived in this paper are applied to solve a fourth-order boundary value problem. Numerical results are given to illustrate the efficiency of our methods and compared with exact solution.
文摘In this paper, we extend the reliable modification of the Adomian Decom-position Method coupled to the Lesnic’s approach to solve boundary value problems and initial boundary value problems with mixed boundary conditions for linear and nonlinear partial differential equations. The method is applied to different forms of heat and wave equations as illustrative examples to exhibit the effectiveness of the method. The method provides the solution in a rapidly convergent series with components that can be computed iteratively. The numerical results for the illustrative examples obtained show remarkable agreement with the exact solutions. We also provide some graphical representations for clear-cut comparisons between the solutions using Maple software.
