88,904 research outputs found
Oscillation-free method for semilinear diffusion equations under noisy initial conditions
Noise in initial conditions from measurement errors can create unwanted
oscillations which propagate in numerical solutions. We present a technique of
prohibiting such oscillation errors when solving initial-boundary-value
problems of semilinear diffusion equations. Symmetric Strang splitting is
applied to the equation for solving the linear diffusion and nonlinear
remainder separately. An oscillation-free scheme is developed for overcoming
any oscillatory behavior when numerically solving the linear diffusion portion.
To demonstrate the ills of stable oscillations, we compare our method using a
weighted implicit Euler scheme to the Crank-Nicolson method. The
oscillation-free feature and stability of our method are analyzed through a
local linearization. The accuracy of our oscillation-free method is proved and
its usefulness is further verified through solving a Fisher-type equation where
oscillation-free solutions are successfully produced in spite of random errors
in the initial conditions.Comment: 19 pages, 9 figure
A discrete least squares collocation method for two-dimensional nonlinear time-dependent partial differential equations
In this paper, we develop regularized discrete least squares collocation and
finite volume methods for solving two-dimensional nonlinear time-dependent
partial differential equations on irregular domains. The solution is
approximated using tensor product cubic spline basis functions defined on a
background rectangular (interpolation) mesh, which leads to high spatial
accuracy and straightforward implementation, and establishes a solid base for
extending the computational framework to three-dimensional problems. A
semi-implicit time-stepping method is employed to transform the nonlinear
partial differential equation into a linear boundary value problem. A key
finding of our study is that the newly proposed mesh-free finite volume method
based on circular control volumes reduces to the collocation method as the
radius limits to zero. Both methods produce a large constrained least-squares
problem that must be solved at each time step in the advancement of the
solution. We have found that regularization yields a relatively
well-conditioned system that can be solved accurately using QR factorization.
An extensive numerical investigation is performed to illustrate the
effectiveness of the present methods, including the application of the new
method to a coupled system of time-fractional partial differential equations
having different fractional indices in different (irregularly shaped) regions
of the solution domain
A space-time pseudospectral discretization method for solving diffusion optimal control problems with two-sided fractional derivatives
We propose a direct numerical method for the solution of an optimal control
problem governed by a two-side space-fractional diffusion equation. The
presented method contains two main steps. In the first step, the space variable
is discretized by using the Jacobi-Gauss pseudospectral discretization and, in
this way, the original problem is transformed into a classical integer-order
optimal control problem. The main challenge, which we faced in this step, is to
derive the left and right fractional differentiation matrices. In this respect,
novel techniques for derivation of these matrices are presented. In the second
step, the Legendre-Gauss-Radau pseudospectral method is employed. With these
two steps, the original problem is converted into a convex quadratic
optimization problem, which can be solved efficiently by available methods. Our
approach can be easily implemented and extended to cover fractional optimal
control problems with state constraints. Five test examples are provided to
demonstrate the efficiency and validity of the presented method. The results
show that our method reaches the solutions with good accuracy and a low CPU
time.Comment: This is a preprint of a paper whose final and definite form is with
'Journal of Vibration and Control', available from
[http://journals.sagepub.com/home/jvc]. Submitted 02-June-2018; Revised
03-Sept-2018; Accepted 12-Oct-201
An advanced meshless method for time fractional diffusion equation
Recently, because of the new developments in sustainable engineering and renewable energy, which are usually governed by a series of fractional partial differential equations (FPDEs), the numerical modelling and simulation for fractional calculus are attracting more and more attention from researchers. The current dominant numerical method for modeling FPDE is Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings including difficulty in simulation of problems with the complex problem domain and in using irregularly distributed nodes. Because of its distinguished advantages, the meshless method has good potential in simulation of FPDEs. This paper aims to develop an implicit meshless collocation technique for FPDE. The discrete system of FPDEs is obtained by using the meshless shape functions and the meshless collocation formulation. The stability and convergence of this meshless approach are investigated theoretically and numerically. The numerical examples with regular and irregular nodal distributions are used to validate and investigate accuracy and efficiency of the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling and simulation of fractional partial differential equations
- …