scholarly journals A finite difference method for an initial–boundary value problem with a Riemann–Liouville–Caputo spatial fractional derivative

2021 ◽  
Vol 381 ◽  
pp. 113020
Author(s):  
José Luis Gracia ◽  
Martin Stynes
2020 ◽  
Vol 20 (4) ◽  
pp. 815-825 ◽  
Author(s):  
Chaobao Huang ◽  
Xiaohui Liu ◽  
Xiangyun Meng ◽  
Martin Stynes

AbstractAn initial-boundary value problem, whose differential equation contains a sum of fractional time derivatives with orders between 0 and 1, is considered. Its spatial domain is {(0,1)^{d}} for some {d\in\{1,2,3\}}. This problem is a generalisation of the problem considered by Stynes, O’Riordan and Gracia in SIAM J. Numer. Anal. 55 (2017), pp. 1057–1079, where {d=1} and only one fractional time derivative was present. A priori bounds on the derivatives of the unknown solution are derived. A finite difference method, using the well-known L1 scheme for the discretisation of each temporal fractional derivative and classical finite differences for the spatial discretisation, is constructed on a mesh that is uniform in space and arbitrarily graded in time. Stability and consistency of the method and a sharp convergence result are proved; hence it is clear how to choose the temporal mesh grading in a optimal way. Numerical results supporting our theoretical results are provided.


Author(s):  
А.М. Слиденко ◽  
В.М. Слиденко

Приводится анализ механических колебаний элементов ударного устройства с помощью модели стержневого типа. Ударник и инструмент связаны упругими и диссипативными элементами, которые имитируют их взаимодействие. Аналогично моделируется взаимодействие инструмента с рабочей средой. Сформулирована начально-краевая задача для системы двух волновых уравнений с учетом переменных поперечных сечений стержней. Площади поперечных сечений определяются параметрическими формулами при сохранении объемов стержней. Параметрические формулы позволяют получать различного вида зависимости площади поперечного сечения стержня от его длины. Начальные условия отражают физическую картину взаимодействия инструмента с ударником и рабочей средой. Краевые условия описывают контактные взаимодействия ударника с инструментом и последнего с рабочей средой. В качестве модельной задачи рассматривается соударение ударника и инструмента через элемент большой жесткости. Начально-краевая задача исследуется разностным методом. Проводится сравнение решений задачи, полученных с помощью двухслойной и трехслойной разностных схем. Такие схемы реализованы в общей компьютерной программе в системе Mathcad. Показано, что при вычислениях распределения нормальных напряжений по длине стержня лучшими свойствами относительно устойчивости обладает двухслойная схема The article gives the analysis of mechanical vibrations of the impact device elements using the model of the rod type. The hammer and the tool are connected by elastic and dissipative elements that simulate their interaction. The interaction of the tool with the processing medium is simulated in a similar way. An initial boundary-value problem is formulated for a system of two wave equations taking into account the variable cross sections of the rods. Cross-sectional areas are determined by parametric formulas maintaining the volume of the rods. Parametric formulas allow one to obtain various dependence types of the cross-sectional area of the rod on its length. The initial and boundary conditions reflect the physical phenomenon of the tool interaction with the processing medium, and also describe the contact interactions of the hammer with the tool. The impacting of the hammer and the tool through an element of high rigidity is considered as a model problem. To control the limiting values, the solution of the model problem by the Fourier method is used. The initial-boundary-value problem is investigated by the difference method. A comparison of solutions obtained for the two-layer and three-layer difference schemes is given. Such schemes are realized in a common computer program in the Mathcad. It is shown that the two-layer scheme has the best properties in relation to stability while calculating the distribution of normal voltage along the length of the rod


2021 ◽  
Vol 65 (3) ◽  
pp. 166-174

An initial-boundary value problem for the subdiffusion equation with an elliptic operator A(D) in RN is studied in the article. Existence and uniqueness theorems for the problem under study are proved by the Fourier method. Considering the order of the Riemann-Liouville time-fractional derivative as an unknown parameter, an inverse problem of determining this parameter is investigated. Likewise, the initial-boundary value problem was considered in the case of replacing the operator A(D) with its power Aσ.Then, existence and uniqueness theorems were proved for the solution of the inverse problem of determining the order of the fractional derivative and the power σ.


Sign in / Sign up

Export Citation Format

Share Document