Diffusion-wave Equation Research Articles

A Local Radial Basis Function Method for Numerical Approximation of Multidimensional Multi-Term Time-Fractional Mixed Wave-Diffusion and Subdiffusion Equation Arising in Fluid Mechanics

This article develops a simple hybrid localized mesh-free method (LMM) for the numerical modeling of new mixed subdiffusion and wave-diffusion equation with multi-term time-fractional derivatives. Unlike conventional multi-term fractional wave-diffusion or subdiffusion equations, this equation features a unique time–space coupled derivative while simultaneously incorporating both wave-diffusion and subdiffusion terms. Our proposed method follows three basic steps: (i) The given equation is transformed into a time-independent form using the Laplace transform (LT); (ii) the LMM is then used to solve the transformed equation in the LT domain; (iii) finally, the time domain solution is obtained by inverting the LT. We use the improved Talbot method and the Stehfest method to invert the LT. The LMM is used to circumvent the shape parameter sensitivity and ill-conditioning of interpolation matrices that commonly arise in global mesh-free methods. Traditional time-stepping methods achieve accuracy only with very small time steps, significantly increasing the computational time. To overcome these shortcomings, the LT is used to provide a more powerful alternative by removing the need for fine temporal discretization. Additionally, the Ulam–Hyers stability of the considered model is analyzed. Four numerical examples are presented to illustrate the effectiveness and practical applicability of the method.

Quantitative photothermal investigation of nonradiative recombination parameters in GaAs/InAs(QD)/GaAs quantum dot structures using a three-layer laser beam deflection model

In this paper, we developed a theoretical model for the photothermal deflection technique in order to investigate the electronic parameters of three-layer semiconductor structures. This model is based on the resolution of thermal and photogenerated carrier diffusion-wave equations in different media. Theoretical results show that the amplitude and phase of the photothermal deflection signal is very sensitive to the nonradiative recombination parameters. The theoretical model is applied to one layer of InAs quantum dots (QDs) inserted in GaAs matrix InAs/GaAs QDs in order to investigate the QD density effects on nonradiative recombination parameters in InAs through fitting the theoretical photothermal beam deflection signal to the experimental data. It was found that the minority carrier lifetime and the electronic diffusivity decrease as functions of increasing InAs QD density. This result is also related to the decrease in the mobility from 21.58 to 4.17 (±12.9%) cm2/V s and the minority carrier diffusion length from 0.62 (±5.8%) to 0.14 (±10%) μm, respectively. Furthermore, both interface recombination velocities S2/3 of GaAs/InAs (QDs) and S1/2 of InAs (QDs)/GaAs increase from 477.7 (±6.2%) to 806.5 (±4%) cm/s and from 75 (±7.8%) to 148.1 (±5.5%) cm/s, respectively.

A Novel and Accurate Algorithm for Solving Fractional Diffusion-Wave Equations

The main objective of this work is to apply a novel and accurate algorithm for solving the second-order and fourth-order fractional diffusion-wave equations (FDWEs). First, the desired equation is reduced to the corresponding Volterra integral equation (VIE). Then, the collocation method is applied, for which the Chebyshev cardinal functions (CCFs) have been considered as the bases. In this paper, the CCFs based on a Lobatto grid are introduced and used for the first time to solve these kinds of equations. To this end, the derivative and fractional integral operators are represented in CCFs. The main features of the method are simplicity, compliance with boundary conditions, and good accuracy. An exact analysis to show the convergence of the scheme is presented, and illustrative examples confirm our investigation.

An efficient temporal approximation for weakly singular time-fractional semilinear diffusion-wave equation with variable coefficients

An efficient temporal approximation for weakly singular time-fractional semilinear diffusion-wave equation with variable coefficients

A novel explicit fast numerical scheme for the Cauchy problem for integro-differential equations with a difference kernel and its application

The present study focuses on designing a second-order novel explicit fast numerical scheme for the Cauchy problem incorporating memory associated with an evolutionary equation, where the integral term's kernel is a discrete difference operator. The Cauchy problem under consideration is related to a real finite-dimensional Hilbert space and includes a self-adjoint operator that is both positive and definite. We introduce a transformative technique for converting the Cauchy problem incorporating memory, into a local evolutionary system of equations by approximating the difference kernel using the sum of exponentials (SoE) approach. A second-order explicit scheme is then constructed to solve the local system. We thoroughly investigate the stability of this explicit scheme, and present the necessary conditions for the stability of the scheme. Moreover, we extended our investigation to encompass time-fractional diffusion-wave equations (TFDWEs) involving a fractional Caputo derivative with an order ranging between (1,2). Initially, we transform the main TFDWE model into a new model that incorporates the fractional Riemann-Liouville integral. Subsequently, we expand the applicability of our idea to develop an explicit fast numerical algorithm for approximating the model. The stability properties of this fast scheme for solving TFDWEs are assessed. Numerical simulations including a two-dimensional Cauchy problem as well as one-dimensional and two-dimensional TFDWE models are provided to validate the accuracy and experimental order of convergence of the schemes.

Numerical simulation of an effective transform mechanism with convergence analysis of the fractional diffusion-wave equations

In the current study, we solve two very important mathematical models, such as the time fractional-order space-fractional telegraph and diffusion-wave equations using a reliable technique called the Adomian decomposition natural method (ADNM), which combines Adomian decomposition and natural transform. The diffusion wave equation describes the flood wave propagation, which is used in solving overland and open channel flow problems. For this reason, it is critical to fully understand and effectively solve the diffusion wave equations. Because telegraph equations are crucial for modeling and developing voltage or frequency transmission, they are widely used in physics and engineering. Furthermore, the designing process is greatly impacted by the uncertainty in the system parameters. For nonlinear ordinary differential equations based on the theorem of Banach fixed point, we provide existence and uniqueness theorem proofs. The present approach has been successfully used to explore exact solutions for time fractional-order and space fractional-order applications. The results show how effective and valuable the ADNM. This paper presents a methodology that will be used in future work to address similar nonlinear problems related to fractional calculus.

H3N3-2$$_\sigma $$-based difference schemes for time multi-term fractional diffusion-wave equation

H3N3-2$$_\sigma $$-based difference schemes for time multi-term fractional diffusion-wave equation

Mainardi smoothing homotopy method for solving nonlinear optimal control problems

This paper proposes a new type of smoothing homotopy method for solving general nonlinear optimal control problems via the indirect method, leveraging the Mainardi kernel as the smoothing kernel. The Mainardi kernel is firstly derived from the fundamental solution of the time-fractional diffusion-wave equation, representing a generalized form of the Gaussian kernel. By altering the fractional derivative order, the kernel can seamlessly switch between non-Gaussian and Gaussian forms. Then, parts of the two-point boundary value problems are convolved with the smoothing kernel, and the resulting surrogates are incorporated into the necessary conditions, replacing the terminal state and costate variables. The homotopy process is used for the smoothing parameter: increasing this parameter enhances the smoothing effect, simplifying the homotopy problems, while a parameter value of zero implies no smoothing, reverting to the original problems. Additionally, explicit expressions for the Mainardi kernel at specific derivative orders are derived, avoiding the inefficiencies associated with fractional derivatives. Simulation examples demonstrate that the proposed method offers significant advantages in flexibility and convergence compared to the traditional Gaussian smoothing homotopy method.

A New Perspective for Scientific Modelling: Sparse Reconstruction Based Approach for Learning Time-Space Fractional Differential Equations

Abstract This paper studies a sparse reconstruction based approach to learn time-space fractional differential equations, i.e. to identify parameter values and particularly the order of the fractional derivatives. The approach uses a generalized Taylor series expansion to generate, in every iteration, a feature matrix which is used to learn the fractional orders of both, temporal and spatial derivatives by minimizing the LASSO operator using Differential Evolution algorithm. To verify the robustness of the method, numerical results for time-space fractional diffusion equation, wave equation, and Burgers? equation at different noise levels in the data are presented. Finally, the methodology is applied to a realistic example where underlying fractional differential equation associated to published experimental data obtained from an in-vitro cell culture assay is learned.

Analysis of a meshless generalized finite difference method for the time-fractional diffusion-wave equation

In this paper, a generalized finite difference method (GFDM) is proposed and analyzed for meshless numerical solution of the time-fractional diffusion-wave equation. Two (3−α)-order accurate temporal discretization schemes are presented by using the L1 formula and the original H2N2 or fast H2N2 formulas to discretize the time-fractional derivative of order α∈(1,2). The stability of the temporal discretization schemes is analyzed. Then, the time-fractional diffusion-wave initial-boundary value problem is transformed into a series of time-independent integer-order boundary value problems, and discrete linear algebraic systems are built by the application of the GFDM. Accuracy analysis of the GFDM with both original H2N2 and fast H2N2 formulas is presented in theory, and numerical experimental results are provided to verify the theoretical results and the effectiveness of the proposed meshless method.

Regularization of an inverse source problem for fractional diffusion-wave equations under a general noise assumption

Regularization of an inverse source problem for fractional diffusion-wave equations under a general noise assumption

Enhancing flood wave modelling of reservoir failure: a comparative study of structure-from-motion based 2D and 3D methodologies

Predicting flood wave propagation from reservoir failures is critical to practical flood hazard assessment and risk management. Flood waves are sensitive to topography, channel geometry, structures, and natural features along floodplain paths. Thus, the accuracy of flood wave modelling depends on how precisely those features are represented. This study introduces an enhancing approach to flood wave modelling by accurately representing three-dimensional objects in floodplains using the structure-from-motion (SfM). This method uses an unmanned aerial vehicle to capture topographic complexities and account for ground objects that impact flood propagation. Using the three-dimensional volume of fluid numerical approach significantly improves an enhanced representation of turbulent flow dynamics and computational efficiency, especially in handling large topography datasets. Reproductions from this enhanced three-dimensional approach were validated against recent reservoir failure observations and contrasted with traditional two-dimensional models. The results revealed that the suggested three-dimensional methodology achieved a significant 84.4% reproducibility when juxtaposed with actual inundation traces. It was 35.5%p more accurate than the two-dimensional diffusion wave equation (DWE) and 17.1%p more than the shallow water equation (SWE) methods in predicting flood waves. This suggests that the reproducibility of the DWE and SWE decreases compared to the three-dimensional approach when considering more complex floodplains. These results demonstrate that three-dimensional flood wave analysis with the SfM methodology is optimal for effectively minimising topographic and flood wave reproduction errors across extensive areas. This dual reduction in errors significantly enhances the reliability of flood hazard assessments and improves risk management by providing more precise and realistic predictions of flood waves.

Superconvergence analysis of spectral volume methods for one-dimensional diffusion and third-order wave equations

Superconvergence analysis of spectral volume methods for one-dimensional diffusion and third-order wave equations

Meshless analysis of fractional diffusion-wave equations by generalized finite difference method

In this paper, a meshless generalized finite difference method (GFDM) is proposed to solve the time fractional diffusion-wave (TFDW) equations. A second-order temporal discretization scheme is developed to tackle the Caputo fractional derivative, and then spatial discretization formulas are derived by the GFDM. Theoretical accuracy and convergence of the GFDM for TFDW equations are analyzed. Numerical results verify the theoretical results and the efficiency of the method.

Convergence analysis of a fully discrete scheme for diffusion-wave equation forced by tempered fractional Brownian motion

In this paper, we investigate a fully discrete approximation of stochastic diffusion-wave equation driven by additive tempered fractional Gaussian noise. This additive noise exhibits semi-long range dependence. The model involves two nonlocal terms in time, i.e., a Caputo fractional derivative and a tempered fractional Brownian motion. The space discretization is achieved via the spectral Galerkin method. The Caputo fractional derivative of order α∈(1,2) is approximated by using the Grünwald–Letnikov formula. Combining Mittag-Leffler function, Laplace transform and z-transform, we establish the error estimates of the full discretization in the sense of mean-squared L2-norm. Numerical experiments are presented to confirm the strong convergence rates.

A discrete spectral method for time fractional fourth-order 2D diffusion-wave equation involving [formula omitted]-Caputo fractional derivative

In this work, the ψ-Caputo fractional derivative, as a generalization of the classical Caputo fractional derivative in which the fractional derivative of a sufficiently differentiable function is defined with respect to another strictly increasing function, ψ, is utilized to define the time fractional fourth-order 2D diffusion-wave equation. A Chebyshev–Gauss–Lobatto scheme is developed to solve this equation. In this way, a new operational matrix for the ψ-Riemann–Liouville fractional integration of the Chebyshev polynomials is derived. The solution of the equation is obtained by determining the solution of the algebraic system extracted from approximation the ψ-Caputo fractional derivative term via a finite discrete Chebyshev series and employing the expressed operational matrix. The validity of the established approach is examined by solving two examples.

DECAY RATES AND INITIAL VALUES FOR TIME-FRACTIONAL DIFFUSION-WAVE EQUATIONS

We consider a solution u(·,t) to an initial boundary value problem for time-fractional diffusion-wave equation with the order α ∈ (0,2)\ {1} where t is a time variable. We first prove that a suitable norm of u(·,t) is bounded by the rate t −α for 0 α 1 and t 1−α for 1 α 2 for all large t > 0. Second, we characterize initial values in the cases where the decay rates are faster than the above critical exponents. Differently from the classical diffusion equation α = 1, the decay rate can keep some local characterization of initial values. The proof is based on the eigenfunction expansions of solutions and the asymptotic expansions of the Mittag-Leffler functions for large time.

Temporal second-order difference schemes for the nonlinear time-fractional mixed sub-diffusion and diffusion-wave equation with delay

This paper investigates a nonlinear time-fractional mixed sub-diffusion and diffusion-wave equation with delay. The problem is particularly challenging due to its nonlinear nature, the presence of a time delay, and the incorporation of both fractional diffusion and fractional wave terms, introducing computational complexities for numerical analysis. To address this, we transform the model into a new generalized form involving Riemann–Liouville fractional integrals and a Caputo derivative of order α∈(0,1). Subsequently, a temporal second-order linearized difference scheme is presented to approximate the model, and its unconditional stability is rigorously proven based on discrete Gronwall’s inequality. To reduce computational and storage costs, we extend the discussion to a fast variant of the proposed difference scheme. The sum of exponentials approach is employed to approximate the kernel function in fractional-order operators, leading to fast difference analogs for the Riemann–Liouville fractional integral and the Caputo derivative. A fast variant of the developed direct method is introduced based on these analogs. Numerical results are provided to validate our theoretical findings and assess the accuracy and efficiency of the difference schemes.

Runoff Hydrograph Analysis of HEC-RAS 2D Flow Hydrodynamics Meteorological Rain-on-Grid on Observed Watershed: A Case Study of Wiroko Sub-Watershed

In contemporary hydrological analysis, numerical models have emerged as an alternative solution for water resources planning and management, particularly in forecasting extreme events and mitigating disaster risks. HEC-RAS 2D flow hydrodynamics enables the modeling of rainfall with input from meteorological boundary condition data. Precipitation within the studied watershed can be modeled as point rainfall or area-averaged rainfall in the rain-on-grid feature of HEC-RAS 2D flow hydrodynamics. In this study, the HEC-RAS 2D flow model utilizes the diffusion wave equations (DWE) for simulating unsteady flow routing. When modeling a complex watershed, the 2D hydrodynamics model has become a viable alternative, especially concerning the watershed's physical characteristics. This study aims to assess the reliability of HEC-RAS 2D flow hydrodynamics and compare its results with the hydrograph data from the observed watershed in the Wiroko Sub-Watershed of the upper Wonogiri Dam. In this research, extreme rainfall events of 100 mm in the first hour are simulated and compared with the hydrograph data from the observed watershed. Based on the numerical model results in the Wiroko Sub-Watershed, it was determined that conforming to the land use with a Manning roughness value of 0.12 (designated as developed area-medium intensity) resulted in a peak discharge (Qp) difference of 0.8%. Meanwhile, the time to peak (Tp) value exhibited a discrepancy of 1.93 hours longer between the numerical model and the observed hydrograph.

The Solution of Caputo Fractional Partial Differential Equations By Sumudu Transform

In this paper, we obtained an explicit solutions of the fractional diffusion-wave equations involving partial Caputo fractional derivative by using Sumudu transform. The solutions of the fractional diffusion-wave equations obtained in terms of Mittag-Leffler function and generalized Wright function.Some illustrative examples are also given.