Comparative Studies between Picard’s and Taylor’s Methods of Numerical Solutions of First Ordinary Order Differential Equations Arising from Real-Life Problems

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.

Efficient Numerical Solution of the Modified Mild-Slope Equation

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.

Numerical solutions to heat transfer of nanofluid flow over stretching sheet subjected to variations of nanoparticle volume fraction and wall temperature

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.

Numerical Solution of Constrained Mechanical System Motions Equations and Inverse Problems of Dynamics

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.

ON THE NUMERICAL SOLUTION OF QUASILINEAR WAVE EQUATION WITH STRONG DISSIPATIVE TERM

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.

NUMERICAL SOLUTIONS OF DISCONTINUOUS BOUNDARY VALUE PROBLEMS FOR GENERAL ELLIPTIC COMPLEX EQUATIONS OF FIRST ORDER

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.

Existence and Numerical Solution for a Coupled System of Multi-term Fractional Differential Equations

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.

Numerical Solution of the Weight Function for Electromagnetic Flowmeter

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.

Numerical Solution for Fractional-Order Differential Systems with Time-Domain and Frequency-Domain Methods

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.

ON NUMERICAL SOLUTIONS OF PERIODICALLY PERTURBED CONSERVATIVE SYSTEMS

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.

Approximate Formulation and Numerical Solution for Hypersingular Boundary Integral Equations in Plane Elasticity

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.

Numerical solution of a mathematical model for water wavesin large coastal areas

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.

Governing equations and numerical solutions of tension leg platform with finite amplitude motion

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.

AN APPLICATION OF THE EXPONENTIAL CUBIC SPLINES TO NUMERICAL SOLUTION OF A SELF-ADJOINT PERTURBATION PROBLEM

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.

NUMERICAL SOLUTION OF QUASILINEAR SINGULARLY PERTURBED ORDINARY DIFFERENTIAL EQUATION WITHOUT TURNING POINTS

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 .

NUMERICAL SOLUTION OF THE SINGULARLY PERTURBED PROBLEM FOR THE HYPERBOLIC EQUATION WITH INITIAL JUMP

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.

THE STABILITY AND CONVERGENCE OF THE FINITE ANALYTIC METHOD FOR THE NUMERICAL SOLUTION OF CONVECTIVE DIFFUSION EQUATION

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.

Numerical Solution for Fractional Partial Differential Equation with Bernstein Polynomials

A framework to obtain numerical solution of the fractional partial differential equation using Bernstein polynomials is presented. The main characteristic behind this approach is that a fractional order operational matrix of Bernstein polynomials is derived. With the operational matrix, the equation is transformed into the products of several dependent matrixes which can also be regarded as the system of linear equations after dispersing the variable. By solving the linear equations, the numerical solutions are acquired. Only a small number of Bernstein polynomials are needed to obtain a satisfactory result. Numerical examples are provided to show that the method is computationally efficient.

The θ－Methods in Numerical Solution of Systems of Differential Equations with Two Delay Terms

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.