To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’...To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.展开更多
An efficient numerical model for wave refraction, diffraction and reflection is presented in this paper. In the model the modified time-dependent mild-slope equation is transformed into an evolution equation and an im...An efficient numerical model for wave refraction, diffraction and reflection is presented in this paper. In the model the modified time-dependent mild-slope equation is transformed into an evolution equation and an improved ADI method involving a relaxation factor is adopted to solve it. The method has the advantage of improving the numerical stability and convergence rate by properly determining the relaxation factor. The range of the relaxation factor making the differential scheme unconditionally stable is determined by stability analysis. Several verifications are performed to examine the accuracy of the present model. The numerical results coincide with the analytic solutions or experimental data very well and the computer time is reduced.展开更多
The numerical analysis of heat transfer of laminar nanofluid flow over a fiat stretching sheet is presented. Two sets of boundary conditions (BCs) axe analyzed, i.e., a constant (Case 1) and a linear streamwise va...The numerical analysis of heat transfer of laminar nanofluid flow over a fiat stretching sheet is presented. Two sets of boundary conditions (BCs) axe analyzed, i.e., a constant (Case 1) and a linear streamwise variation of nanopaxticle volume fraction and wall temperature (Case 2). The governing equations and BCs axe reduced to a set of nonlinear ordinary differential equations (ODEs) and the corresponding BCs, respectively. The dependencies of solutions on Prandtl number Pr, Lewis number Le, Brownian motion number Nb, and thermophoresis number Nt are studied in detail. The results show that the reduced Nusselt number and the reduced Sherwood number increase for the BCs of Case 2 compared with Case 1. The increases of Nb, Nt, and Le numbers cause a decrease of the reduced Nusselt number, while the reduced Sherwood number increases with the increase of Nb and Le numbers. For low Prandtl numbers, an increase of Nt number can cause to decrease in the reduced Sherwood number, while it increases for high Prandtl numbers.展开更多
In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarant...In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.展开更多
The numerical solution for a type of quasilinear wave equation is studied.The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved.The...The numerical solution for a type of quasilinear wave equation is studied.The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved.The error of the difference solution is estimated.The theoretical results are controlled on a numerical example.展开更多
This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solutio...This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.展开更多
In this paper, authors discuss the numerical methods of general discontinuous boundary value problems for elliptic complex equations of first order, They first give the well posedness of general discontinuous boundary...In this paper, authors discuss the numerical methods of general discontinuous boundary value problems for elliptic complex equations of first order, They first give the well posedness of general discontinuous boundary value problems, reduce the discontinuous boundary value problems to a variation problem, and then find the numerical solutions of above problem by the finite element method. Finally authors give some error-estimates of the foregoing numerical solutions.展开更多
An initial value problem was considered for a coupled differential system with multi-term Caputo type fractional derivatives. By means of nonlinear alternative of Leray-Schauder and Banach contraction principle,the ex...An initial value problem was considered for a coupled differential system with multi-term Caputo type fractional derivatives. By means of nonlinear alternative of Leray-Schauder and Banach contraction principle,the existence and uniqueness of solutions for the system were derived. Using a fractional predictorcorrector method, a numerical method was presented for the specified system. An example was given to illustrate the obtained results.展开更多
A numerical solution of the weight function for a two-electrode electromagnetic flowmeter was proposed. The solution was obtained by using the finite element method based on the basic equation of a traditional two-ele...A numerical solution of the weight function for a two-electrode electromagnetic flowmeter was proposed. The solution was obtained by using the finite element method based on the basic equation of a traditional two-electrode electromagnetic flowmeter. The two-dimensional distribution of the weight function of the electromagnetic flowmeter obtained was verified by the analytical solution. Three-dimensional distribution of the weight function was also presented in the paper. It can be employed to analyze the sensitivity and linearity of the electromagnetic flowmeter with non-uniform magnetic field, and even to assist the design of the excitation coil pair.展开更多
For a general nonlinear fractional-order differential equation, the numerical solution is a good way to approximate the trajectory of such systems. In this paper, a novel algorithm for numerical solution of fractional...For a general nonlinear fractional-order differential equation, the numerical solution is a good way to approximate the trajectory of such systems. In this paper, a novel algorithm for numerical solution of fractional-order differential equations based on the definition of Grunwald-Letnikov is presented. The results of numerical solution by using the novel method and the frequency-domain method are compared, and the limitations of frequency-domain method are discussed.展开更多
A nonlinear perturbed conservative system is discussed. BY means of Hadamard's theorem. the existence and uniqueness of the solution of the continuous problem arc proved. When the equation is discreted on the unif...A nonlinear perturbed conservative system is discussed. BY means of Hadamard's theorem. the existence and uniqueness of the solution of the continuous problem arc proved. When the equation is discreted on the uniform meshes, it is proved that the corresponding discrete problem has a unique solution, Finally, the accuracy of the numerical solution is considered and a simple algorithm is provided for solving the nonlinear difference equation.展开更多
In this paper, the variable-coefficient diffusion-advection (DA) equation, which arises in modeling various physical phenomena, is studied by the Lie symmetry approach. The similarity reductions are derived by deter...In this paper, the variable-coefficient diffusion-advection (DA) equation, which arises in modeling various physical phenomena, is studied by the Lie symmetry approach. The similarity reductions are derived by determining the complete sets of point symmetries of this equation, and then exact and numerical solutions are reported for the reduced second-order nonlinear ordinary differential equations. Further, an extended (Gl/G)-expansion method is applied to the DA equation to construct some new non-traveling wave solutions.展开更多
Based on the fact that the singular boundary integrals in the sense of Cauchy principal value can be represented approximately by the mean values of two companion nearly singular boundary integrals, a vary general app...Based on the fact that the singular boundary integrals in the sense of Cauchy principal value can be represented approximately by the mean values of two companion nearly singular boundary integrals, a vary general approach was developed in the paper. In the approach, the approximate formulation before discretization was constructed to cope with the difficulties encountered in the corner treatment in the formulations of hypersingular boundary integral equations. This makes it possible to solve the hypersingular boundary integral equation numerically in a non regularized form and in a local manner by using conforming C 0 quadratic boundary elements and standard Gaussian quadratures similar to those employed in the conventional displacement BIE formulations. The approximate formulation is very convenient to use because the corner information is comprised naturally in the representations of those approximate integrals. Numerical examples in plane elasticity show that with the present approach, the compatible or better results can be achieved in comparison with those of the conventional BIE formulations.展开更多
Based on the evolution equation for water waves, a mathematical model for wave propaga tion in large mild - slope areas is derived. The model is solved by the finite difference method with the staggered grid system. ...Based on the evolution equation for water waves, a mathematical model for wave propaga tion in large mild - slope areas is derived. The model is solved by the finite difference method with the staggered grid system. The computational results are in good agreement with experimental data and show that the model can obtain better results with relatively coarser grids. The model can be used to simulate water wave propagation in large coastal areas and can be efficiently solved without much programming effort.展开更多
It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP,...It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP, The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theo- retical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displace- ment, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force, Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.展开更多
We present a class of the second order optimal splines difference schemes derived from ex- ponential cubic splines for self-adjoint singularly perturbed 2-point boundary value problem. We prove an optimal error estima...We present a class of the second order optimal splines difference schemes derived from ex- ponential cubic splines for self-adjoint singularly perturbed 2-point boundary value problem. We prove an optimal error estimate and give illustrative numerical example.展开更多
In this paper we consider a quasilinear second order ordinary diferential equation with a small parameter Firstly an approximate problem is constructed. Then an iterative procedure is developed. Finally we give an alg...In this paper we consider a quasilinear second order ordinary diferential equation with a small parameter Firstly an approximate problem is constructed. Then an iterative procedure is developed. Finally we give an algorithm whose accuracy is good for arbitrary e>0 .展开更多
In this paper we consider the initial-boundary value problem for a second order hyperbolic equation with initial jump. The bounds on the derivatives of the exact solution are given. Then a difference scheme is constru...In this paper we consider the initial-boundary value problem for a second order hyperbolic equation with initial jump. The bounds on the derivatives of the exact solution are given. Then a difference scheme is constructed on a non-uniform grid. Finally, uniform convergence of the difference solution is proved in the sense of the discrete energy norm.展开更多
In this paper we make a close study of the finite analytic method by means of the maximum principles in differential equations and give the proof of the stability and convergence of the finite analytic method.
文摘To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.
文摘An efficient numerical model for wave refraction, diffraction and reflection is presented in this paper. In the model the modified time-dependent mild-slope equation is transformed into an evolution equation and an improved ADI method involving a relaxation factor is adopted to solve it. The method has the advantage of improving the numerical stability and convergence rate by properly determining the relaxation factor. The range of the relaxation factor making the differential scheme unconditionally stable is determined by stability analysis. Several verifications are performed to examine the accuracy of the present model. The numerical results coincide with the analytic solutions or experimental data very well and the computer time is reduced.
文摘The numerical analysis of heat transfer of laminar nanofluid flow over a fiat stretching sheet is presented. Two sets of boundary conditions (BCs) axe analyzed, i.e., a constant (Case 1) and a linear streamwise variation of nanopaxticle volume fraction and wall temperature (Case 2). The governing equations and BCs axe reduced to a set of nonlinear ordinary differential equations (ODEs) and the corresponding BCs, respectively. The dependencies of solutions on Prandtl number Pr, Lewis number Le, Brownian motion number Nb, and thermophoresis number Nt are studied in detail. The results show that the reduced Nusselt number and the reduced Sherwood number increase for the BCs of Case 2 compared with Case 1. The increases of Nb, Nt, and Le numbers cause a decrease of the reduced Nusselt number, while the reduced Sherwood number increases with the increase of Nb and Le numbers. For low Prandtl numbers, an increase of Nt number can cause to decrease in the reduced Sherwood number, while it increases for high Prandtl numbers.
基金Supported by Russian Fund of Fund amental Investigations(Pr.990101064)and Russian Minister of Educatin
文摘In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.
文摘The numerical solution for a type of quasilinear wave equation is studied.The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved.The error of the difference solution is estimated.The theoretical results are controlled on a numerical example.
文摘This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.
文摘In this paper, authors discuss the numerical methods of general discontinuous boundary value problems for elliptic complex equations of first order, They first give the well posedness of general discontinuous boundary value problems, reduce the discontinuous boundary value problems to a variation problem, and then find the numerical solutions of above problem by the finite element method. Finally authors give some error-estimates of the foregoing numerical solutions.
基金National Natural Science Foundation of China(No.11371087)
文摘An initial value problem was considered for a coupled differential system with multi-term Caputo type fractional derivatives. By means of nonlinear alternative of Leray-Schauder and Banach contraction principle,the existence and uniqueness of solutions for the system were derived. Using a fractional predictorcorrector method, a numerical method was presented for the specified system. An example was given to illustrate the obtained results.
基金Supported by Chinese National Programs for High Technology Research and Development(2008AA042207)
文摘A numerical solution of the weight function for a two-electrode electromagnetic flowmeter was proposed. The solution was obtained by using the finite element method based on the basic equation of a traditional two-electrode electromagnetic flowmeter. The two-dimensional distribution of the weight function of the electromagnetic flowmeter obtained was verified by the analytical solution. Three-dimensional distribution of the weight function was also presented in the paper. It can be employed to analyze the sensitivity and linearity of the electromagnetic flowmeter with non-uniform magnetic field, and even to assist the design of the excitation coil pair.
基金the Natural Science Foundation of CQ CSTC under Grant No. 2007BB2161.
文摘For a general nonlinear fractional-order differential equation, the numerical solution is a good way to approximate the trajectory of such systems. In this paper, a novel algorithm for numerical solution of fractional-order differential equations based on the definition of Grunwald-Letnikov is presented. The results of numerical solution by using the novel method and the frequency-domain method are compared, and the limitations of frequency-domain method are discussed.
文摘A nonlinear perturbed conservative system is discussed. BY means of Hadamard's theorem. the existence and uniqueness of the solution of the continuous problem arc proved. When the equation is discreted on the uniform meshes, it is proved that the corresponding discrete problem has a unique solution, Finally, the accuracy of the numerical solution is considered and a simple algorithm is provided for solving the nonlinear difference equation.
文摘In this paper, the variable-coefficient diffusion-advection (DA) equation, which arises in modeling various physical phenomena, is studied by the Lie symmetry approach. The similarity reductions are derived by determining the complete sets of point symmetries of this equation, and then exact and numerical solutions are reported for the reduced second-order nonlinear ordinary differential equations. Further, an extended (Gl/G)-expansion method is applied to the DA equation to construct some new non-traveling wave solutions.
文摘Based on the fact that the singular boundary integrals in the sense of Cauchy principal value can be represented approximately by the mean values of two companion nearly singular boundary integrals, a vary general approach was developed in the paper. In the approach, the approximate formulation before discretization was constructed to cope with the difficulties encountered in the corner treatment in the formulations of hypersingular boundary integral equations. This makes it possible to solve the hypersingular boundary integral equation numerically in a non regularized form and in a local manner by using conforming C 0 quadratic boundary elements and standard Gaussian quadratures similar to those employed in the conventional displacement BIE formulations. The approximate formulation is very convenient to use because the corner information is comprised naturally in the representations of those approximate integrals. Numerical examples in plane elasticity show that with the present approach, the compatible or better results can be achieved in comparison with those of the conventional BIE formulations.
基金National Natural Science Foundation of China under contract No.59839330by ChinaPostdoctoral Science Foundation.
文摘Based on the evolution equation for water waves, a mathematical model for wave propaga tion in large mild - slope areas is derived. The model is solved by the finite difference method with the staggered grid system. The computational results are in good agreement with experimental data and show that the model can obtain better results with relatively coarser grids. The model can be used to simulate water wave propagation in large coastal areas and can be efficiently solved without much programming effort.
基金Project supported by "Creativeness Project of the Tenth Five-Year Plan" of Chinese Academy of Sciences (No.KJCX2-SW-L03)the National High-Tech Research and Development Program of China (863 Program) (No.2004AA617010)
文摘It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP, The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theo- retical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displace- ment, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force, Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.
文摘We present a class of the second order optimal splines difference schemes derived from ex- ponential cubic splines for self-adjoint singularly perturbed 2-point boundary value problem. We prove an optimal error estimate and give illustrative numerical example.
文摘In this paper we consider a quasilinear second order ordinary diferential equation with a small parameter Firstly an approximate problem is constructed. Then an iterative procedure is developed. Finally we give an algorithm whose accuracy is good for arbitrary e>0 .
文摘In this paper we consider the initial-boundary value problem for a second order hyperbolic equation with initial jump. The bounds on the derivatives of the exact solution are given. Then a difference scheme is constructed on a non-uniform grid. Finally, uniform convergence of the difference solution is proved in the sense of the discrete energy norm.
文摘In this paper we make a close study of the finite analytic method by means of the maximum principles in differential equations and give the proof of the stability and convergence of the finite analytic method.