An Efficient Legendre Spectral Tau Matrix Formulation for Solving Fractional Subdiffusion and Reaction Subdiffusion Equations

2015 ◽  
Vol 10 (2) ◽  
Author(s):  
E. H. Doha ◽  
A. H. Bhrawy ◽  
S. S. Ezz-Eldien
Keyword(s):  
Boundary Conditions ◽  
Spectral Method ◽  
Operational Matrix ◽  
Neumann Boundary ◽  
Variable Coefficients ◽  
Dirichlet Boundary Conditions ◽  
Fractional Integrals ◽  
Algebraic Equations ◽  
Matrix Formulation ◽  
Subdiffusion Equation

In this work, we discuss an operational matrix approach for introducing an approximate solution of the fractional subdiffusion equation (FSDE) with both Dirichlet boundary conditions (DBCs) and Neumann boundary conditions (NBCs). We propose a spectral method in both temporal and spatial discretizations for this equation. Our approach is based on the space-time shifted Legendre tau-spectral method combined with the operational matrix of fractional integrals, described in the Riemann–Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations in the unknown expansion coefficients of the sought-for spectral approximations. In addition, this approach is also investigated for solving the FSDE with the variable coefficients and the fractional reaction subdiffusion equation (FRSDE). For conforming the validity and accuracy of the numerical scheme proposed, four numerical examples with their approximate solutions are presented. Also, comparisons between our numerical results and those obtained by compact finite difference method (CFDM), Box-type scheme (B-TS), and FDM with Fourier analysis (FA) are introduced.

Numerical approximations for fractional diffusion equations via a Chebyshev spectral-tau method

Open Physics ◽  
2013 ◽  
Vol 11 (10) ◽  
Author(s):  
Eid Doha ◽  
Ali Bhrawy ◽  
Samer Ezz-Eldien
Keyword(s):  
Numerical Technique ◽  
Operational Matrix ◽  
Fractional Diffusion ◽  
High Accuracy ◽  
Variable Coefficients ◽  
Diffusion Equations ◽  
Algebraic Equations ◽  
Fractional Differentiation ◽  
Fractional Diffusion Equations ◽  
General Method

AbstractIn this paper, a class of fractional diffusion equations with variable coefficients is considered. An accurate and efficient spectral tau technique for solving the fractional diffusion equations numerically is proposed. This method is based upon Chebyshev tau approximation together with Chebyshev operational matrix of Caputo fractional differentiation. Such approach has the advantage of reducing the problem to the solution of a system of algebraic equations, which may then be solved by any standard numerical technique. We apply this general method to solve four specific examples. In each of the examples considered, the numerical results show that the proposed method is of high accuracy and is efficient for solving the time-dependent fractional diffusion equations.

On the vacuum stress induced by uniform acceleration or supporting the ether

1977 ◽  
Vol 354 (1676) ◽  
pp. 79-99 ◽  
Keyword(s):  
Scalar Field ◽  
Boundary Conditions ◽  
Neumann Boundary ◽  
Scalar Fields ◽  
Black Body ◽  
Dirichlet Boundary Conditions ◽  
Massless Scalar ◽  
Perfect Conductor ◽  
Infinite Plane ◽  
Uniform Acceleration

The stresses induced in the vacuum by the uniform acceleration of an infinite plane conductor are computed for the massless scalar and electromagnetic fields. Both Dirichlet and Neumann boundary conditions are considered for the scalar field; far from the conductor it is found, independently of the boundary condition, that the vacuum stress is ‘local’ and corresponds to the absence from the vacuum of black body radiation. Approaching the conductor, the energy density in the Dirichlet case is slightly lower than the ‘local’ term, and in the Neumann case slightly higher. At very small distances it again has the same asymptotic form for both scalar fields. For the electromagnetic field the results are similar to those for the scalar field with Dirichlet boundary conditions. Far from the conductor the spectrum is again black-body, though not Planckian. In all cases the acausal nature of ‘ perfect conductor ’ boundary conditions prevents the stress tensor from being finite on the conductor.

Asymptotic behaviour of some reaction-diffusion systems modelling complex combustion on bounded domains

1993 ◽  
Vol 123 (6) ◽  
pp. 1151-1163
Author(s):  
Joel D. Avrin
Keyword(s):  
Steady State ◽  
Boundary Conditions ◽  
Neumann Boundary ◽  
Reaction Diffusion ◽  
Convergence Result ◽  
Decay Rates ◽  
Dirichlet Boundary Conditions ◽  
Reaction Diffusion Systems ◽  
Convergence Results ◽  
Mass Fractions

SynopsisWe consider three models of multiple-step combustion processes on bounded spatial domains. Previously, steady-state convergence results have been established for these models with zero Neumann boundary conditions imposed on the temperature as well as the mass fractions. We retain here throughout the same boundary conditions on the mass fractions, but in our first set of results we establish steady-state convergence results with fixed Dirichlet boundary conditions on the temperature. Next, under certain physically reasonable assumptions, we develop, for two of the models, estimates on the decay rates of both mass fractions to zero, while for the remaining model we develop estimates on the decay rate of one concentration to zero and establish a positive lower bound on the other mass fraction. These results hold under either set of boundary conditions, but when the Dirichlet conditions are imposed on the temperature, we are able to obtain estimates on the rate of convergence of the temperature to its (generally nonconstant) steady-state. Finally, we improve the results of a previous paper by adding a temperature convergence result.

scholarly journals Parameters for the Collapse of Turbulence in the Stratified Plane Couette Flow

2018 ◽  
Vol 75 (9) ◽  
pp. 3211-3231 ◽  
Author(s):  
Ivo G. S. van Hooijdonk ◽  
Herman J. H. Clercx ◽  
Cedrick Ansorge ◽  
Arnold F. Moene ◽  
Bas J. H. van de Wiel
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Boundary Conditions ◽  
Couette Flow ◽  
Neumann Boundary ◽  
Surface Flux ◽  
Shear Capacity ◽  
Buoyancy Flux ◽  
Dirichlet Boundary Conditions ◽  
Obukhov Length

Abstract We perform direct numerical simulation of the Couette flow as a model for the stable boundary layer. The flow evolution is investigated for combinations of the (bulk) Reynolds number and the imposed surface buoyancy flux. First, we establish what the similarities and differences are between applying a fixed buoyancy difference (Dirichlet) and a fixed buoyancy flux (Neumann) as boundary conditions. Moreover, two distinct parameters were recently proposed for the turbulent-to-laminar transition: the Reynolds number based on the Obukhov length and the “shear capacity,” a velocity-scale ratio based on the buoyancy flux maximum. We study how these parameters relate to each other and to the atmospheric boundary layer. The results show that in a weakly stratified equilibrium state, the flow statistics are virtually the same between the different types of boundary conditions. However, at stronger stratification and, more generally, in nonequilibrium conditions, the flow statistics do depend on the type of boundary condition imposed. In the case of Neumann boundary conditions, a clear sensitivity to the initial stratification strength is observed because of the existence of multiple equilibriums, while for Dirichlet boundary conditions, only one statistically steady turbulent equilibrium exists for a particular set of boundary conditions. As in previous studies, we find that when the imposed surface flux is larger than the maximum buoyancy flux, no turbulent steady state occurs. Analytical investigation and simulation data indicate that this maximum buoyancy flux converges for increasing Reynolds numbers, which suggests a possible extrapolation to the atmospheric case.

Wavelets Galerkin Method for the Fractional Subdiffusion Equation

2016 ◽  
Vol 11 (6) ◽  
Author(s):  
M. H. Heydari
Keyword(s):  
Galerkin Method ◽  
Fractional Derivatives ◽  
Operational Matrix ◽  
Practical Importance ◽  
Time Derivative ◽  
Algebraic Equations ◽  
Legendre Wavelets ◽  
Equivalent Problem ◽  
Subdiffusion Equation ◽  
The Galerkin Method

The time fractional subdiffusion equation (FSDE) as a class of anomalous diffusive systems has obtained by replacing the time derivative in ordinary diffusion by a fractional derivative of order 0<α<1. Since analytically solving this problem is often impossible, proposing numerical methods for its solution has practical importance. In this paper, an efficient and accurate Galerkin method based on the Legendre wavelets (LWs) is proposed for solving this equation. The time fractional derivatives are described in the Riemann–Liouville sense. To do this, we first transform the original subdiffusion problem into an equivalent problem with fractional derivatives in the Caputo sense. The LWs and their fractional operational matrix (FOM) of integration together with the Galerkin method are used to transform the problem under consideration into the corresponding linear system of algebraic equations, which can be simply solved to achieve the solution of the problem. The proposed method is very convenient for solving such problems, since the initial and boundary conditions are taken into account, automatically. Furthermore, the efficiency of the proposed method is shown for some concrete examples. The results reveal that the proposed method is very accurate and efficient.

On the Number of Solutions to Semilinear Boundary Value Problems

2004 ◽  
Vol 4 (3) ◽  
Author(s):  
Markus Kunze ◽  
Rafael Ortega
Keyword(s):  
Boundary Conditions ◽  
Boundary Value Problems ◽  
Upper Bound ◽  
Dirichlet Boundary ◽  
Neumann Boundary ◽  
Neumann Boundary Conditions ◽  
Dirichlet Boundary Conditions ◽  
Number Of Solutions ◽  
Semilinear Elliptic ◽  
Semilinear Elliptic Problems

AbstractWe consider semilinear elliptic problems of the form Δu + g(u) = f(x) with Neumann boundary conditions or Δu+λ1u+g(u) = f(x) with Dirichlet boundary conditions, and we derive conditions on g and f under which an upper bound on the number of solutions can be obtained.

scholarly journals Synchronization for a class of generalized neural networks with interval time-varying delays and reaction-diffusion terms

2016 ◽  
Vol 21 (3) ◽  
pp. 379-399 ◽  
Author(s):  
Qintao Gan ◽  
Tielin Liu ◽  
Chang Liu ◽  
Tianshi Lv
Keyword(s):  
Neural Networks ◽  
Boundary Conditions ◽  
Neumann Boundary ◽  
Reaction Diffusion ◽  
Dirichlet Boundary Conditions ◽  
Linear Matrix ◽  
Time Varying ◽  
Interval Time ◽  
Generalized Neural Networks ◽  
Special Cases

In this paper, the synchronization problem for a class of generalized neural networks with interval time-varying delays and reaction-diffusion terms is investigated under Dirichlet boundary conditions and Neumann boundary conditions, respectively. Based on Lyapunov stability theory, both delay-derivative-dependent and delay-range-dependent conditions are derived in terms of linear matrix inequalities (LMIs), whose solvability heavily depends on the information of reaction-diffusion terms. The proposed generalized neural networks model includes reaction-diffusion local field neural networks and reaction-diffusion static neural networks as its special cases. The obtained synchronization results are easy to check and improve upon the existing ones. In our results, the assumptions for the differentiability and monotonicity on the activation functions are removed. It is assumed that the state delay belongs to a given interval, which means that the lower bound of delay is not restricted to be zero. Finally, the feasibility and effectiveness of the proposed methods is shown by simulation examples.

scholarly journals Non-universal Casimir Forces at Approach to Bose–Einstein Condensation of an Ideal Gas: Effect of Dirichlet Boundary Conditions

2020 ◽  
Vol 181 (3) ◽  
pp. 944-951
Author(s):  
M. Napiórkowski ◽  
J. Piasecki ◽  
J. W. Turner
Keyword(s):  
Boundary Conditions ◽  
Dirichlet Boundary ◽  
Neumann Boundary ◽  
Ideal Gas ◽  
Einstein Condensate ◽  
Dirichlet Boundary Conditions ◽  
Einstein Condensation ◽  
Casimir Forces ◽  
Leading Term ◽  
Bose Einstein

Abstract We analyze the Casimir forces for an ideal Bose gas enclosed between two infinite parallel walls separated by the distance D. The walls are characterized by the Dirichlet boundary conditions. We show that if the thermodynamic state with Bose–Einstein condensate present is correctly approached along the path pertinent to the Dirichlet b.c. then the leading term describing the large-distance decay of thermal Casimir force between the walls is $$\sim 1/D^{2}$$ ∼ 1 / D 2 with a non-universal amplitude. The next order correction is $$\sim \ln D/D^3$$ ∼ ln D / D 3 . These observations remain in contrast with the decay law for both the periodic and Neumann boundary conditions for which the leading term is $$\sim 1/D^3$$ ∼ 1 / D 3 with a universal amplitude. We associate this discrepancy with the D-dependent positive value of the one-particle ground state energy in the case of Dirichlet boundary conditions.

Sign in / Sign up

Export Citation Format

Share Document