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.

INVERSE PROBLEM FOR THE VISCOELASTIC MEDIUM WITH DISCONTINUOUS WAVE IMPEDANCE

In this paper, the inverse problem for the viscoelastic medium is investigated in the time domain, in which the wave impedance of the medium is discontinuous at the rear interface. The differentio-integral equations governing the behavior of the scattering and propagation operators are utilized to reconstruct the relaxation modulus of the viscoelastic medium. A new approach, in which only the one-side measurement reflection data for one round trip through the viscoelastic layer, is developed. The numerical examples are given at the end of the paper. Ir is shown that the curves of the reconstructed moduli coincide very well with the original relaxation moduli.

An Upwind Mixed Finite Volume Element-fractional Step Method and Convergence Analysis for Three-dimensional Compressible Contamination Treatment from Nuclear Waste

In this paper the authors discuss a numerical simulation problem of three-dimensional compressible contamination treatment from nuclear waste. The mathematical model, a nonlinear convection-diffusion system of four PDEs, determines four major physical unknowns: the pressure, the concentrations of brine and radionuclide, and the temperature. The pressure is solved by a conservative mixed finite volume element method, and the computational accuracy is improved for Darcy velocity. Other unknowns are computed by a composite scheme of upwind approximation and mixed finite volume element. Numerical dispersion and nonphysical oscillation are eliminated, and the convection-dominated diffusion problems are solved well with high order computational accuracy. The mixed finite volume element is conservative locally, and get the objective functions and their adjoint vector functions simultaneously. The conservation nature is an important character in numerical simulation of underground fluid. Fractional step difference is introduced to solve the concentrations of radionuclide factors, and the computational work is shortened significantly by decomposing a three-dimensional problem into three successive one-dimensional problems. By the theory and technique of a priori estimates of differential equations, we derive an optimal order estimates in L2norm. Finally, numerical examples show the effectiveness and practicability for some actual problems.

TWO-GRID CHARACTERISTIC FINITE VOLUME METHODS FOR NONLINEAR PARABOLIC PROBLEMS＊

In this work, two-grid characteristic finite volume schemes for the nonlinear parabolic problem are considered. In our algorithms, the diffusion term is discretized by the finite volume method, while the temporal differentiation and advection terms are treated by the characteristic scheme. Under some conditions about the coefficients and exact solution, optimal error estimates for the numerical solution are obtained. Furthermore, the two- grid characteristic finite volume methods involve solving a nonlinear equation on coarse mesh with mesh size H, a large linear problem for the Oseen two-grid characteristic finite volume method on a fine mesh with mesh size h = O（H2） or a large linear problem for the Newton two-grid characteristic finite volume method on a fine mesh with mesh size h = 0（I log hll/2H3）. These methods we studied provide the same convergence rate as that of the characteristic finite volume method, which involves solving one large nonlinear problem on a fine mesh with mesh size h. Some numerical results are presented to demonstrate the efficiency of the proposed methods.

An Upwind Mixed Volume Element-Fractional Step Method on a Changing Mesh for Compressible Contamination Treatment from Nuclear Waste

In this paper the authors discuss the numerical simulation problem of threedimensional compressible contamination treatment from nuclear waste.The mathematical model is defined by an initial-boundary nonlinear convection-diffusion system of four partial differential equations:a parabolic equation for the pressure,two convection-diffusion equations for the concentrations of brine and radionuclide and a heat conduction equation for the temperature.The pressure appears within the concentration equations and heat conduction equation,and the Darcy velocity controls the concentrations and the temperature.The pressure is solved by the conservative mixed volume element method,and the order of the accuracy is improved by the Darcy velocity.The concentration of brine and temperature are computed by the upwind mixed volume element method on a changing mesh,where the diffusion is discretized by a mixed volume element and the convection is treated by an upwind scheme.The composite method can solve the convection-dominated diffusion problems well because it eliminates numerical dispersion and nonphysical oscillation and has high order computational accuracy.The mixed volume element has the local conservation of mass and energy,and it can obtain the brine and temperature and their adjoint vector functions simultaneously.The conservation nature plays an important role in numerical simulation of underground fluid.The concentrations of radionuclide factors are solved by the method of upwind fractional step difference and the computational work is decreased by decomposing a three-dimensional problem into three successive one-dimensional problems and using the method of speedup.By the theory and technique of a priori estimates of differential equations,we derive an optimal order result in L^(2) norm.Numerical examples are given to show the effectiveness and practicability and the composite method is testified as a powerful tool to solve the well-known actual problem.