Numerical analysis for viscoelastic fluid flow with distributed/variable order time fractional Maxwell constitutive models

Fractional calculus has been widely used to study the flow of viscoelastic fluids recently,and fractional differential equations have attracted a lot of attention.However,the research has shown that the fractional equation with constant order operators has certain limitations in characterizing some physical phenomena.In this paper,the viscoelastic fluid flow of generalized Maxwell fluids in an infinite straight pipe driven by a periodic pressure gradient is investigated systematically.Consider the complexity of the material structure and multi-scale effects in the viscoelastic fluid flow.The modified time fractional Maxwell models and the corresponding governing equations with distributed/variable order time fractional derivatives are proposed.Based on the L1-approximation formula of Caputo fractional derivatives,the implicit finite difference schemes for the distributed/variable order time fractional governing equations are presented,and the numerical solutions are derived.In order to test the correctness and availability of numerical schemes,two numerical examples are established to give the exact solutions.The comparisons between the numerical solutions and the exact solutions have been made,and their high consistency indicates that the present numerical methods are effective.Then,this paper analyzes the velocity distributions of the distributed/variable order fractional Maxwell governing equations under specific conditions,and discusses the effects of the weight coefficient(α)in distributed order time fractional derivatives,the orderα(r,t)in variable fractional order derivatives,the relaxation timeλ,and the frequencyωof the periodic pressure gradient on the fluid flow velocity.Finally,the flow rates of the distributed/variable order fractional Maxwell governing equations are also studied.

On Existence of Entropy Solution for a Doubly Nonlinear Differential Equation with L1 -Data

We consider a class of doubly nonlinear history-dependent problems having a convection term and a pseudomonotone nonlinear diffusion operator associated an equation of the type ?t(k * (b(v) - b(v0))) - div(a(x,Dv) + F(v)) = f where the right hand side belongs to L1. The kernel k belongs to the large class of PC kernels. In particular, the case of fractional time derivatives of order α ∈ (0,1) is included. Assuming b nondecreasing with L1-data, we prove existence in the framework of entropy solutions. The approach adopted for the proof is based on a several step approximation method and by using a result in the case of a strictly increasing b.

The Self-similar Solution to Some Nonlinear Integro-differential Equations Corresponding to Fractional Order Time Derivative

In this paper we study the self-similar solution to a class of nonlinear integro-differential equations which correspond to fractional order time derivative and interpolate nonlinear heat and wave equation. Using the space-time estimates which were established by Hirata and Miao in [1] we prove the global existence of self-similar solution of Cauchy problem for the nonlinear integro-differential equation in C＊（[0,∞];B^8pp,∞（R^n）.

Memory dependent anomalous diffusion in comb structure under distributed order time fractional dual-phase-lag model

This paper considers a novel distributed order time fractional dual-phase-lag model to analyze the anomalous diffusion in a comb structure,which has a widespread application in medicine and biology.The newly proposed constitution model is a generalization of the dual-phase-lag model,in which a spectrum of the time fractional derivatives with the memory characteristic governed by the weight coefficient is considered and the formulated governing equation contains both the diffusion and wave characteristics.With the L1-formula to discrete the time Caputo fractional derivatives,the finite difference method is used to discretize the model and the related numerical results are plotted graphically.By adding a source term,an exact solution is defined to verify the correctness of the numerical scheme and the convergence order of the error in spatial direction is presented.Finally,the dynamic characteristics of the particle distributions and the effects of involved parameters on the total number of particles in the x-direction are analyzed in detail.