scholarly journals Numerically Efficient Methods for Variational Fractional Wave Equations: An Explicit Four-Step Scheme

Mathematics ◽  
2019 ◽  
Vol 7 (11) ◽  
pp. 1095
Author(s):  
Jorge E. Macías-Díaz

In this work, we investigate numerically a one-dimensional wave equation in generalized form. The system considers the presence of constant damping and functional anomalous diffusion of the Riesz type. Reaction terms are also considered, in such way that the mathematical model can be presented in variational form when damping is not present. As opposed to previous efforts available in the literature, the reaction terms are not only functions of the solution. Instead, we consider the presence of smooth functions that depend on fractional derivatives of the solution function. Using a finite-difference approach, we propose a numerical scheme to approximate the solutions of the fractional wave equation. Along with this integrator, we propose discrete forms of the local and the total energy operators. In a first stage, we show rigorously that the energy properties of the continuous system are mimicked by our discrete methodology. In particular, we prove that the discrete system is dissipative (respectively, conservative) when damping is present (respectively, absent), in agreement with the continuous model. The theoretical numerical analysis of this system is more complicated in light of the presence of the functional form of the anomalous diffusion. To solve this problem, some novel technical lemmas are proved and used to establish the stability and the quadratic convergence of the scheme. Finally, we provide some computer simulations to show the capability of the scheme to conserve/dissipate the energy. Various fractional problems with functional forms of the anomalous diffusion of the solution are considered to that effect.

Author(s):  
N.H. Sweilam ◽  
T.A. Assiri

In this paper, the space fractional wave equation (SFWE) is numerically studied, where the fractional derivative is defined in the sense of Caputo. An explicit finite difference approximation (EFDA) for SFWE is presented. The stability and the error analysis of the EFDA are discussed. To demonstrate the effectiveness of the approximated method, some test examples are presented.   


Author(s):  
Peter Straka ◽  
Mark Meerschaert ◽  
Robert McGough ◽  
Yuzhen Zhou

AbstractFractional wave equations with attenuation have been proposed by Caputo [5], Szabo [28], Chen and Holm [7], and Kelly et al. [11]. These equations capture the power-law attenuation with frequency observed in many experimental settings when sound waves travel through inhomogeneous media. In particular, these models are useful for medical ultrasound. This paper develops stochastic solutions and weak solutions to the power law wave equation of Kelly et al. [11].


2021 ◽  
Vol 5 (1) ◽  
pp. 314-336
Author(s):  
Tristram de Piro ◽  

We clarify some arguments concerning Jefimenko’s equations, as a way of constructing solutions to Maxwell’s equations, for charge and current satisfying the continuity equation. We then isolate a condition on non-radiation in all inertial frames, which is intuitively reasonable for the stability of an atomic system, and prove that the condition is equivalent to the charge and current satisfying certain relations, including the wave equations. Finally, we prove that with these relations, the energy in the electromagnetic field is quantised and displays the properties of the Balmer series.


Author(s):  
P. Metallidis ◽  
S. Natsiavas

Abstract The present study generalises previous research work on the dynamics of discrete oscillators with piecewise linear characteristics and investigates the response of a continuous model system with clearance and motion-limiting constraints. More specifically, in the first part of this work, an analysis is presented for determining exact periodic response of a periodically excited deformable rod, whose motion is constrained by a flexible obstacle. This methodology is based on the exact solution form obtained within response intervals where the system parameters remain constant and its behavior is governed by a linear equation of motion. The unknowns of the problem are subsequently determined by imposing an appropriate set of periodicity and matching conditions. The analytical part is complemented by a suitable method for determining the stability properties of the located periodic motions. In the second part of the study, the analysis is applied to several cases in order to investigate the effect of the system parameters on its dynamics. Special emphasis is placed on comparing these results with results obtained for similar but rigid rods. Finally, direct integration of the equation of motion in selected areas reveals the existence of motions, which are more complicated than the periodic motions determined analytically.


Mathematics ◽  
2020 ◽  
Vol 8 (6) ◽  
pp. 874
Author(s):  
Francesco Iafrate ◽  
Enzo Orsingher

In this paper we study the time-fractional wave equation of order 1 < ν < 2 and give a probabilistic interpretation of its solution. In the case 0 < ν < 1 , d = 1 , the solution can be interpreted as a time-changed Brownian motion, while for 1 < ν < 2 it coincides with the density of a symmetric stable process of order 2 / ν . We give here an interpretation of the fractional wave equation for d > 1 in terms of laws of stable d−dimensional processes. We give a hint at the case of a fractional wave equation for ν > 2 and also at space-time fractional wave equations.


Geophysics ◽  
2020 ◽  
Vol 85 (1) ◽  
pp. T1-T13 ◽  
Author(s):  
Ning Wang ◽  
Tieyuan Zhu ◽  
Hui Zhou ◽  
Hanming Chen ◽  
Xuebin Zhao ◽  
...  

The spatial derivatives in decoupled fractional Laplacian (DFL) viscoacoustic and viscoelastic wave equations are the mixed-domain Laplacian operators. Using the approximation of the mixed-domain operators, the spatial derivatives can be calculated by using the Fourier pseudospectral (PS) method with barely spatial numerical dispersions, whereas the time derivative is often computed with the finite-difference (FD) method in second-order accuracy (referred to as the FD-PS scheme). The time-stepping errors caused by the FD discretization inevitably introduce the accumulative temporal dispersion during the wavefield extrapolation, especially for a long-time simulation. To eliminate the time-stepping errors, here, we adopted the [Formula: see text]-space concept in the numerical discretization of the DFL viscoacoustic wave equation. Different from existing [Formula: see text]-space methods, our [Formula: see text]-space method for DFL viscoacoustic wave equation contains two correction terms, which were designed to compensate for the time-stepping errors in the dispersion-dominated operator and loss-dominated operator, respectively. Using theoretical analyses and numerical experiments, we determine that our [Formula: see text]-space approach is superior to the traditional FD-PS scheme mainly in three aspects. First, our approach can effectively compensate for the time-stepping errors. Second, the stability condition is more relaxed, which makes the selection of sampling intervals more flexible. Finally, the [Formula: see text]-space approach allows us to conduct high-accuracy wavefield extrapolation with larger time steps. These features make our scheme suitable for seismic modeling and imaging problems.


Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. C85-C94 ◽  
Author(s):  
Houzhu (James) Zhang ◽  
Hongwei Liu ◽  
Yang Zhao

Seismic anisotropy is an intrinsic elastic property. Appropriate accounting of anisotropy is critical for correct and accurate positioning seismic events in reverse time migration. Although the full elastic wave equation may serve as the ultimate solution for modeling and imaging, pseudoelastic and pseudoacoustic wave equations are more preferable due to their computation efficiency and simplicity in practice. The anisotropic parameters and their relations are not arbitrary because they are constrained by the energy principle. Based on the investigation of the stability condition of the pseudoelastic wave equations, we have developed a set of explicit formulations for determining the S-wave velocity from given Thomsen’s parameters [Formula: see text] and [Formula: see text] for vertical transverse isotropy and tilted transverse isotropy media. The estimated S-wave velocity ensures that the wave equations are stable and well-posed in the cases of [Formula: see text] and [Formula: see text]. In the case of [Formula: see text], a common situation in carbonate, a positive value of S-wave velocity is needed to avoid the wavefield instability. Comparing the stability constraints of the pseudoelastic- with the full-elastic wave equation, we conclude that the feasible range of [Formula: see text] and [Formula: see text] was slightly larger for the pseudoelastic assumption. The success of achieving high-accuracy images and high-quality angle gathers using the proposed constraints is demonstrated in a synthetic example and a field example from Saudi Arabia.


2020 ◽  
Vol 14 (4) ◽  
pp. 523-534
Author(s):  
Faiqul Fikri ◽  
Eddy Djauhari ◽  
Endang Rusyaman

Non-linear differential equations with fractional derivative order are mathematical models that are widely used in modeling physical phenomena, one of the applications of these models is non-linear fractional wave equations. Many methods for solving non-linear fractional partial differential equations, one of which is the New Version of Optimal Homotopy Asymptotic Method which is developed by Liaqat Ali in 2016. The author will use this method to solve non-linear fractional wave equations predetermined, so that the convergence of function of the approximation solution non-linear fractional wave equation can be observed and it can be observed that the function of approximation solution of non-linear fractional wave equation solution using the New Version of Optimal Homotopy Asymptotic Method is simple and has a value error using Mean Absolute Percentage Error which is categorized very well


Author(s):  
Jia Wei He ◽  
Yong Zhou

In this paper, we concern with a backward problem for a nonlinear time fractional wave equation in a bounded domain. By applying the properties of Mittag-Leffler functions and the method of eigenvalue expansion, we establish some results about the existence and uniqueness of the mild solutions of the proposed problem based on the compact technique. Due to the ill-posedness of backward problem in the sense of Hadamard, a general filter regularization method is utilized to approximate the solution and further we prove the convergence rate for the regularized solutions.


2020 ◽  
Vol 18 (1) ◽  
pp. 67-86
Author(s):  
Xiaoyong Xu ◽  
Fengying Zhou

Abstract In this paper, a discrete orthogonal spline collocation method combining with a second-order Crank-Nicolson weighted and shifted Grünwald integral (WSGI) operator is proposed for solving time-fractional wave equations based on its equivalent partial integro-differential equations. The stability and convergence of the schemes have been strictly proved. Several numerical examples in one variable and in two space variables are given to demonstrate the theoretical analysis.


Sign in / Sign up

Export Citation Format

Share Document