This article addresses a new numerical method to find a numerical solution of the linear delay differential equation of fractional order , the fractional derivatives described in the Caputo sense. The new approach is to approximating second and third derivatives. A backward finite difference method is used. Besides, the composite Trapezoidal rule is used in the Caputo definition to match the integral term. The accuracy and convergence of the prescribed technique are explained. The results are shown through numerical examples.
This paper examines the finding of spacewise dependent heat source function in pseudoparabolic equation with initial and homogeneous Dirichlet boundary conditions, as well as the final time value / integral specification as additional conditions that ensure the uniqueness solvability of the inverse problem. However, the problem remains ill-posed because tiny perturbations in input data cause huge errors in outputs. Thus, we employ Tikhonov’s regularization method to restore this instability. In order to choose the best regularization parameter, we employ L-curve method. On the other hand, the direct (forward) problem is solved by a finite difference scheme while the inverse one is reformulated as an optimization problem. The
... Show MoreThis paper deals with finding the approximation solution of a nonlinear parabolic boundary value problem (NLPBVP) by using the Galekin finite element method (GFEM) in space and Crank Nicolson (CN) scheme in time, the problem then reduce to solve a Galerkin nonlinear algebraic system(GNLAS). The predictor and the corrector technique (PCT) is applied here to solve the GNLAS, by transforms it to a Galerkin linear algebraic system (GLAS). This GLAS is solved once using the Cholesky method (CHM) as it appear in the matlab package and once again using the Cholesky reduction order technique (CHROT) which we employ it here to save a massive time. The results, for CHROT are given by tables and figures and show
... Show MoreIn this article, an inverse problem of finding timewise-dependent thermal conductivity has been investigated numerically. Numerical solution of forward (direct) problem has been solved by finite-difference method (FDM). Whilst, the inverse (indirect) problem solved iteratively using Lsqnonlin routine from MATLAB. Initial guess for unknown coefficient expressed by explicit relation based on nonlocal overdetermination conditions and intial input data .The obtained numrical results are presented and discussed in several figures and tables. These results are accurate and stable even in the presense of noisy data.
The goal of this research is to develop a numerical model that can be used to simulate the sedimentation process under two scenarios: first, the flocculation unit is on duty, and second, the flocculation unit is out of commission. The general equation of flow and sediment transport were solved using the finite difference method, then coded using Matlab software. The result of this study was: the difference in removal efficiency between the coded model and operational model for each particle size dataset was very close, with a difference value of +3.01%, indicating that the model can be used to predict the removal efficiency of a rectangular sedimentation basin. The study also revealed