LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER

Similar documents
LAPLACE TRANSFORM SOLUTION OF THE PROBLEM OF TIME-FRACTIONAL HEAT CONDUCTION IN A TWO-LAYERED SLAB

International Journal of Pure and Applied Sciences and Technology

Research Article On the Upper Bounds of Eigenvalues for a Class of Systems of Ordinary Differential Equations with Higher Order

Numerical Solution of Nonlinear Multi-order Fractional Differential Equations by Implementation of the Operational Matrix of Fractional Derivative

Lecture 36. Finite Element Methods

4. Eccentric axial loading, cross-section core

Lecture 4: Piecewise Cubic Interpolation

THE COMBINED SHEPARD ABEL GONCHAROV UNIVARIATE OPERATOR

Numerical solution of fractional elliptic PDE's by the collocation method

A Family of Multivariate Abel Series Distributions. of Order k

Two Coefficients of the Dyson Product

The Number of Rows which Equal Certain Row

Reactor Control Division BARC Mumbai India

90 S.S. Drgomr nd (t b)du(t) =u()(b ) u(t)dt: If we dd the bove two equltes, we get (.) u()(b ) u(t)dt = p(; t)du(t) where p(; t) := for ll ; t [; b]:

ON SIMPSON S INEQUALITY AND APPLICATIONS. 1. Introduction The following inequality is well known in the literature as Simpson s inequality : 2 1 f (4)

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Machine Learning Support Vector Machines SVM

Symmetries and Conservation Laws in Classical Mechanics

Numerical Solution of Linear Fredholm Fuzzy Integral Equations by Modified Homotopy Perturbation Method

arxiv: v1 [math.gm] 30 Dec 2015

ORDINARY DIFFERENTIAL EQUATIONS

6 Roots of Equations: Open Methods

Principle Component Analysis

Numerical Solution of Fredholm Integral Equations of the Second Kind by using 2-Point Explicit Group Successive Over-Relaxation Iterative Method

The Schur-Cohn Algorithm

Definition of Tracking

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Partially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism

NON-HOMOGENEOUS COMPOSITE BEAMS: ANALYTIC FORMULATION AND SOLUTION

Many-Body Calculations of the Isotope Shift

Strong Gravity and the BKL Conjecture

Review of linear algebra. Nuno Vasconcelos UCSD

CISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting

COMPLEX NUMBERS INDEX

Applied Statistics Qualifier Examination

A new Approach for Solving Linear Ordinary Differential Equations

Fractional Euler-Lagrange Equations of Order ( α, β ) for Lie Algebroids

Sequences of Intuitionistic Fuzzy Soft G-Modules

Katholieke Universiteit Leuven Department of Computer Science

Identification of Robot Arm s Joints Time-Varying Stiffness Under Loads

CENTROID (AĞIRLIK MERKEZİ )

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)

Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC

A Theoretical Study on the Rank of the Integral Operators for Large- Scale Electrodynamic Analysis

8. INVERSE Z-TRANSFORM

Study of Trapezoidal Fuzzy Linear System of Equations S. M. Bargir 1, *, M. S. Bapat 2, J. D. Yadav 3 1

Quiz: Experimental Physics Lab-I

Computation of Fifth Degree of Spline Function Model by Using C++ Programming

Mechanical resonance theory and applications

APPROXIMATE SOLUTION OF HIGHER ORDER LINEAR DIFFERENTIAL EQUATIONS BY MEANS OF A NEW RATIONAL CHEBYSHEV COLLOCATION METHOD

Dennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1

CHOVER-TYPE LAWS OF THE ITERATED LOGARITHM FOR WEIGHTED SUMS OF ρ -MIXING SEQUENCES

Activator-Inhibitor Model of a Dynamical System: Application to an Oscillating Chemical Reaction System

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

TWO INTERFACIAL COLLINEAR GRIFFITH CRACKS IN THERMO-ELASTIC COMPOSITE MEDIA

STATISTICAL MECHANICS OF THE INVERSE ISING MODEL

CHI-SQUARE DIVERGENCE AND MINIMIZATION PROBLEM

Decomposition of Boolean Function Sets for Boolean Neural Networks

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 9

NUMERICAL MODELLING OF A CILIUM USING AN INTEGRAL EQUATION

Soft Set Theoretic Approach for Dimensionality Reduction 1

Physics for Scientists and Engineers I

7.2 Volume. A cross section is the shape we get when cutting straight through an object.

Effect of Wind Speed on Reaction Coefficient of Different Building Height. Chunli Ren1, a, Yun Liu2,b

Digital Signal Processing

Remember: Project Proposals are due April 11.

Reproducing Kernel Hilbert Space for. Penalized Regression Multi-Predictors: Case in Longitudinal Data

Statistics and Probability Letters

Investigation phase in case of Bragg coupling

Multi-dimensional Central Limit Theorem

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method

Formulated Algorithm for Computing Dominant Eigenvalue. and the Corresponding Eigenvector

Jens Siebel (University of Applied Sciences Kaiserslautern) An Interactive Introduction to Complex Numbers

4. More general extremum principles and thermodynamic potentials

Linear and Nonlinear Optimization

INTRODUCTION TO COMPLEX NUMBERS

ME 501A Seminar in Engineering Analysis Page 1

An Ising model on 2-D image

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede

Solution of Tutorial 5 Drive dynamics & control

Lesson 2. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)

Advanced Machine Learning. An Ising model on 2-D image

Numerical Heat and Mass Transfer

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

Chapter Runge-Kutta 2nd Order Method for Ordinary Differential Equations

New Method for Solving Poisson Equation. on Irregular Domains

The areolar strain concept applied to elasticity

UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

Attribute reduction theory and approach to concept lattice

Electromagnetic modeling of a lightning rod

ESCI 342 Atmospheric Dynamics I Lesson 1 Vectors and Vector Calculus

523 P a g e. is measured through p. should be slower for lesser values of p and faster for greater values of p. If we set p*

ANALOG CIRCUIT SIMULATION BY STATE VARIABLE METHOD

VECTORS VECTORS VECTORS VECTORS. 2. Vector Representation. 1. Definition. 3. Types of Vectors. 5. Vector Operations I. 4. Equal and Opposite Vectors

829. An adaptive method for inertia force identification in cantilever under moving mass

Accurate Instantaneous Frequency Estimation with Iterated Hilbert Transform and Its Application

Work and Energy (Work Done by a Varying Force)

The Dirac equation in Rindler space: A pedagogical introduction

Physics 4B. A positive value is obtained, so the current is counterclockwise around the circuit.

Transcription:

Yn, S.-P.: Locl Frctonl Lplce Seres Expnson Method for Dffuson THERMAL SCIENCE, Yer 25, Vol. 9, Suppl., pp. S3-S35 S3 LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER b Sheng-Png YAN School of Mechncs nd Cvl Engneerng, Chn Unverst of Mnng nd Technolog, Xuzhou, Jngsu, Chn Orgnl scentfc pper DOI:.2298/TSCI463Y In ths pper, we frst propose the locl frctonl Lplce seres expnson method, whch s couplng method of seres expnson method nd Lplce trnsform v locl frctonl dfferentl opertor. An llustrtve exmple for hndlng the dffuson equton rsng n frctl het trnsfer s gven. Ke words: nltcl soluton, dffuson equton, het trnsfer, Lplce seres expnson method, Lplce trnsform Introducton Locl frctonl ntegrl trnsforms hve potentl pplctons for scence nd engneerng [-4]. The ws utlzed to fnd the solutons for dfferentl equtons n the mthemtcl modelng of complex sstems n engneerng to cpture the reltons n spce nd tme wth the kernels wthn non-dfferentblt nd rregulr sets lke frctls [5-]. The locl frctonl Lplce trnsform (LFLT ws ppled to couple other methods, such s decomposton method (DM [5] nd vrtonl terton method (VIM [-9]. Recentl, the locl frctonl seres expnson method (LFSEM ws suggested n [2] nd developed to solve the dfferentl equtons wthn locl frctonl dervtves (LFD [2, 22]. However, the couplng scheme of LFSEM wth LFLT s not consdered. The trget of ths pper s to present the locl frctonl Lplce seres expnson method to del wth the dffuson equton rsng n frctl het trnsfer [23-25]. Fundmentls The locl frctonl ntegrl opertor of ω ( x s defned s [-5, 23]: b j= N ( I ω ( = ω ( (d = lm ω ( ( Γ ( + ( Γ ( + b j = where t = tj+ tj, j =,..., N, t =, tn = b. As the nverse opertor of eq. (, the locl frctonl dervtve of Ω ( s defned s [-5, 23-25]: ( d Ω( [ Ω( Ω( ] Ω ( = lm = = (2 d ( wth [ Ω( Ω( ] Γ ( + Ω [ ( Ω ( ]. Author s e-ml: spn@cumt.edu.cn

S32 Yn, S.-P.: Locl Frctonl Lplce Seres Expnson Method for Dffuson THERMAL SCIENCE, Yer 25, Vol. 9, Suppl., pp. S3-S35 The LFLT of Ω ( s defned s [2-5]: Y, Y E ( { Ω ( } =Ω ( = ( Ω ( (d, < Γ + (3 The nverse LFLT of Ω ( s defned s [2-5]: β + Y, Y, Y E (2π β Ω ( = { Ω ( } = ( Ω ( (d (4 where = β +, nd Re( = β. Some propertes whch re ppled to ths mnuscrpt re [, 4]: Y { Ω ( + bω ( } = Y { Ω ( } + by { Ω ( } (5 2 2 n ( n n ( k ( n k Y{ Ω ( } = Y[ Ω( ] Ω ( (6 k = Y { E ( x } = (7 k Y = ( k + ( k Γ + (8 Anlss of the method We consder gven dfferentl equton n locl form: ( ψ = Κ ψ (9 where ( ψ ( x, /d = ψ nd K s lner locl opertor wth respect to x. We consder mult-term seprted functons of ndependent vrbles t nd x, nmel: ψ( x, = σ( ω ( x ( where σ ( nd ω ( x re two locl frctonl contnuous functons. Settng σ ( = / Γ ( +, we hve: = = Tkng the LFLT of eq. (, we obtn: ψ( x, = ω ( x ( Γ ( + Hence, we obtn: (, = ω ( (2 Γ ( + =

Yn, S.-P.: Locl Frctonl Lplce Seres Expnson Method for Dffuson THERMAL SCIENCE, Yer 25, Vol. 9, Suppl., pp. S3-S35 S33 Y ( { ( x t, } { ( } ( ( Y x ( t t t = ω + = ω + Γ + Γ + (3 = = Y { K (, } K ( ( ( ( = ω = ( K ω = Γ + = Γ + (4 Mkng use of eqs. (3 nd (4, from eq. (9 one obtn: ω+ = = = ( ( ( ( t t ( K ω Γ + Γ + (5 whch leds to the recurson: ω+ ( = ( Kω ( (6 Adoptng the recurson formul (6, we hve: where the convergent condton reds: (, = ω ( (7 Γ ( + = lm ω ( = Γ ( + Hence, the soluton of eq. (9 s determned b: (8 (, = Y { (, } = Y { ω ( } (9 Γ ( + An nltcl soluton for dffuson equton rsng n frctl het trnsfer = We now consder the dffuson equton rsng n frctl het trnsfer [23-25]: ( ( 2 ψ ( x, t ψ ( x, t =, < (2 We present ntl vlues s follows: Adoptng (6, we hve: such tht the recurrence terms re wrtten s: t x ψ ( x, = E ( x (2 ω+ ( = ( Kω( = ω( Y { ( x,} Y { E ( x } = = = (22

S34 Yn, S.-P.: Locl Frctonl Lplce Seres Expnson Method for Dffuson THERMAL SCIENCE, Yer 25, Vol. 9, Suppl., pp. S3-S35 ω ( = (23 2 ω ( = (24 3 ω ( = (25 nd so on. Hence, we get: Tkng nverse LFLT, the non-dfferentble soluton of dffuson equton rsng n frctl het trnsfer cn be wrtten s: ψ( x, = E ( x = Γ ( + = = E ( x E ( (27 nd ts grph s gven n fg.. Conclusons In ths work, we frst hd proposed the couplng scheme of LFSEM wth LFLT, whch clled locl frctonl Lplce seres expnson method (LFLSEM. Bsed on t, (, = (26 Γ ( + = we fnd the non-dfferentble soluton of dffuson equton rsng n frctl het trnsfer. The obtned result shows tht the presented technolog s es, smple, effcent nd ccurte. Fgure. The non-dfferentble soluton of dffuson equton rsng n frctl het trnsfer when the frctl dmenson s equl to ln2/ln3 Nomenclture x spce co-ordntes, [m] Y [ Ω ( ] LFLT of Ω (, [ ] Y, Y [ Ω ( ] nverse LFLT of Y, Ω (, [ ] Greek smbols tme frctl dmensonl order, [ ] tme, [s] ψ(x, concentrton, [ ] References [] Yng, X. J., Locl Frctonl Functonl Anlss & Its Applctons, Asn Acdemc Publsher Lmted, Hong Kong, 2 [2] Yng, X.-J., et l., Locl Frctonl Integrl Trnsforms nd Applctons, Elsever, 25 [3] Srvstv, H. M., et l., Specl Functons n Frctonl Clculus nd Relted Frctonl Dfferntegrl Equtons, World Scentfc, Sngpore, 25 [4] Yng, X. J. Locl Frctonl Integrl Trnsforms, Progress n Nonlner Scence, 4 (2,, pp. -225 [5] Cttn, C., et l., Frctonl Dnmcs, Emergng Scence Publshers, 25

Yn, S.-P.: Locl Frctonl Lplce Seres Expnson Method for Dffuson THERMAL SCIENCE, Yer 25, Vol. 9, Suppl., pp. S3-S35 S35 [6] Zhong, W. P., et l., Applctons of Yng-Fourer Trnsform to Locl Frctonl Equtons wth Locl Frctonl Dervtve nd Locl Frctonl Integrl, Advnced Mterls Reserch, 46 (22, Mrch, pp. 36-3 [7] Yng, A. M., et l., The Yng-Fourer Trnsforms to Het-Conducton n Sem-Infnte Frctl Br, Therml Scence, 7 (23, 3, pp. 77-73 [8] Yng, X. J., et l., A Novel Approch to Processng Frctl Sgnls Usng the Yng-Fourer Trnsforms, Proced Engneerng, 29 (22, Feb., pp. 295-2954 [9] Yng, X. J., et l., Mthemtcl Aspects of the Hesenberg Uncertnt Prncple wthn Locl Frctonl Fourer Anlss, Boundr Vlue Problems, 23 (23,, pp. -6 [] Wng, S. Q., et l., Locl Frctonl Functon Decomposton Method for Solvng Inhomogeneous Wve Equtons wth Locl Frctonl Dervtve, Abstrct nd Appled Anlss, 24 (24, ID 76395 [] Yng, X. J., Locl Frctonl Prtl Dfferentl Equtons wth Frctl Boundr Problems, Advnces n Computtonl Mthemtcs nd ts Applctons, (22,, pp. 6-63 [2] He, J.-H., A Tutorl Revew on Frctl Spcetme nd Frctonl Clculus, Interntonl Journl of Theoretcl Phscs, 53 (24,, pp. 3698-378 [3] Zhng, Y. Z., et l., Intl Boundr Vlue Problem for Frctl Het Equton n the Sem-Infnte Regon b Yng-Lplce Trnsform, Therml Scence, 8 (24, 2, pp. 677-68 [4] Zho, Y., et l., Mppngs for Specl Functons on Cntor Sets nd Specl Integrl Trnsforms v Locl Frctonl Opertors, Abstrct nd Appled Anlss, 23 (23, ID 36978 [5] Zho, C. G., et l., The Yng-Lplce Trnsform for Solvng the IVPs wth Locl Frctonl Dervtve, Abstrct nd Appled Anlss, 24 (24, ID 386459 [6] Lu, C. F., et l., Reconstructve Schemes for Vrtonl Iterton Method wthn Yng-Lplce Trnsform wth Applcton to Frctl Het Conducton Problem, Therml Scence, 7 (23, 3, pp. 75-72 [7] Yng, A. M., et l., Locl Frctonl Lplce Vrtonl Iterton Method for Solvng Lner Prtl Dfferentl Equtons wth Locl Frctonl Dervtve, Dscrete Dnmcs n Nture nd Socet, 24 (24, ID 36598 [8] L, Y., et l., Locl Frctonl Lplce Vrtonl Iterton Method for Frctl Vehculr Trffc Flow, Advnces n Mthemtcl Phscs, 24 (24, ID 64938 [9] Xu, S., et l., Locl Frctonl Lplce Vrtonl Iterton Method for Nonhomogeneous Het Equtons Arsng n Frctl Het Flow, Mthemtcl Problems n Engneerng, 24 (24, ID 94725 [2] Yng, A. M., et l., Locl Frctonl Seres Expnson Method for Solvng Wve nd Dffuson Equtons on Cntor Sets, Abstrct nd Appled Anlss, 23 (23, ID 3557 [2] Zho, Y., et l., Approxmton Solutons for Locl Frctonl Schrödnger Equton n the One-Dmensonl Cntorn Sstem, Advnces n Mthemtcl Phscs, 23 (23, ID 29386 [22] Yng, A. M., et l., Applcton of Locl Frctonl Seres Expnson Method to Solve Klen-Gordon Equtons on Cntor Sets, Abstrct nd Appled Anlss, 24 (24, ID 37274 [23] Yng, X. J., Advnced Locl Frctonl Clculus nd ts Applctons, World Scence, New York, USA, 22 [24] Yng, X. J., et l., Approxmte Solutons for Dffuson Equtons on Cntor Spce-Tme, Proceedngs of the Romnn Acdem, Seres A, 4 (23, 2, pp. 27-33 [25] Ho, Y. J., et l., Helmholtz nd Dffuson Equtons Assocted wth Locl Frctonl Dervtve Opertors Involvng the Cntorn nd Cntor-Tpe Clndrcl Coordntes, Advnces n Mthemtcl Phscs, 23 (23, ID 754248 Pper submtted: October, 24 Pper revsed: Jnur 2, 25 Pper ccepted: Februr 3, 25