Integral Transforms and Dual Integral Equations to Solve Heat Equation with Mixed Conditions

Similar documents
Applied Mathematics Letters

System in Weibull Distribution

Quantum Particle Motion in Physical Space

Approximate Technique for Solving Class of Fractional Variational Problems

A new Approach for Solving Linear Ordinary Differential Equations

A Hybrid Variational Iteration Method for Blasius Equation

Several generation methods of multinomial distributed random number Tian Lei 1, a,linxihe 1,b,Zhigang Zhang 1,c

EXACT TRAVELLING WAVE SOLUTIONS FOR THREE NONLINEAR EVOLUTION EQUATIONS BY A BERNOULLI SUB-ODE METHOD

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

Quasi-Static transient Thermal Stresses in a Robin's thin Rectangular plate with internal moving heat source

Least Squares Fitting of Data

PROBABILITY AND STATISTICS Vol. III - Analysis of Variance and Analysis of Covariance - V. Nollau ANALYSIS OF VARIANCE AND ANALYSIS OF COVARIANCE

EXAMPLES of THEORETICAL PROBLEMS in the COURSE MMV031 HEAT TRANSFER, version 2017

On Pfaff s solution of the Pfaff problem

Determination of the Confidence Level of PSD Estimation with Given D.O.F. Based on WELCH Algorithm

Our focus will be on linear systems. A system is linear if it obeys the principle of superposition and homogenity, i.e.

The Parity of the Number of Irreducible Factors for Some Pentanomials

Special Relativity and Riemannian Geometry. Department of Mathematical Sciences

Denote the function derivatives f(x) in given points. x a b. Using relationships (1.2), polynomials (1.1) are written in the form

1. Statement of the problem

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

Least Squares Fitting of Data

14 The Postulates of Quantum mechanics

Lecture 12: Discrete Laplacian

ACTM State Calculus Competition Saturday April 30, 2011

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

Solving 2D-BKDV Equation by a Sub-ODE Method

Fermi-Dirac statistics

Chapter 1. Theory of Gravitation

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

Numerical Solution of Nonlinear Singular Ordinary Differential Equations Arising in Biology Via Operational Matrix of Shifted Legendre Polynomials

An Optimal Bound for Sum of Square Roots of Special Type of Integers

Study of the possibility of eliminating the Gibbs paradox within the framework of classical thermodynamics *

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Minimization of l 2 -Norm of the KSOR Operator

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

The First Integral Method to Nonlinear Partial Differential Equations

Research Article Green s Theorem for Sign Data

ON THE NUMBER OF PRIMITIVE PYTHAGOREAN QUINTUPLES

Convexity preserving interpolation by splines of arbitrary degree

Scattering by a perfectly conducting infinite cylinder

1 Definition of Rademacher Complexity

High resolution entropy stable scheme for shallow water equations

Lecture Note 3. Eshelby s Inclusion II

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Exact 3-D Solution for System with Rectangular Fins, Part 2

BAYESIAN CURVE FITTING USING PIECEWISE POLYNOMIALS. Dariusz Biskup

Appendix B. The Finite Difference Scheme

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Elastic Collisions. Definition: two point masses on which no external forces act collide without losing any energy.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

Atmospheric Radiation Fall 2008

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Perfect Fluid Cosmological Model in the Frame Work Lyra s Manifold

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Finite Vector Space Representations Ross Bannister Data Assimilation Research Centre, Reading, UK Last updated: 2nd August 2003

The Feynman path integral

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

Department of Economics, Niigata Sangyo University, Niigata, Japan

arxiv: v1 [math.co] 12 Sep 2014

The binomial transforms of the generalized (s, t )-Jacobsthal matrix sequence

Binomial transforms of the modified k-fibonacci-like sequence

Global Sensitivity. Tuesday 20 th February, 2018

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Modified parallel multisplitting iterative methods for non-hermitian positive definite systems

Fixed point method and its improvement for the system of Volterra-Fredholm integral equations of the second kind

PHYS 705: Classical Mechanics. Calculus of Variations II

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

THE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

AN ANALYSIS OF A FRACTAL KINETICS CURVE OF SAVAGEAU

A DISCONTINUOUS LEAST-SQUARES SPATIAL DISCRETIZATION FOR THE S N EQUATIONS. A Thesis LEI ZHU

One-sided finite-difference approximations suitable for use with Richardson extrapolation

Inductance Calculation for Conductors of Arbitrary Shape

Solving Fuzzy Linear Programming Problem With Fuzzy Relational Equation Constraint

A Computational Quadruple Laplace Transform for the Solution of Partial Differential Equations

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter

Digital Signal Processing

Haar Wavelet Collocation Method for the Numerical Solution of Nonlinear Volterra-Fredholm-Hammerstein Integral Equations

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

A FORMULA FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIAGONAL MATRIX

1 Introduction We consider a class of singularly perturbed two point singular boundary value problems of the form: k x with boundary conditions

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Calculus of Variations Basics

The Study of Teaching-learning-based Optimization Algorithm

,..., k N. , k 2. ,..., k i. The derivative with respect to temperature T is calculated by using the chain rule: & ( (5) dj j dt = "J j. k i.

Irreversible Work of Separation and Heat-Driven Separation

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation

International Conference on Advanced Computer Science and Electronics Information (ICACSEI 2013) equation. E. M. E. Zayed and S. A.

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

XII.3 The EM (Expectation-Maximization) Algorithm

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Optimal Pursuit Time in Differential Game for an Infinite System of Differential Equations

Transcription:

Int J Open Probles Copt Math, Vol 7, No 4, Deceber 214 ISSN 1998-6262; Copyrght ICSS Publcaton, 214 www-csrsorg Integral Transfors and Dual Integral Equatons to Solve Heat Equaton wth Mxed Condtons Naser Hoshan Departent of Matheatcs, Tafla Techncal Unversty, PO Box 179, Fax 22533, Tafla Jordan e-al: dr_naserh@yahooco eceved 28 Aprl 214; Accepted 2 October 214 Abstract The paper s devoted to deterne the soluton of the non-statonary heat equaton n a non axal syetrc cylndrcal coordnates subject to xed dscontnuous boundary condtons of the second knd and thrd knds, wth the ad of a fnte Fourer transfor and dual ntegral equatons ethod The soluton of the gven xed proble s ntroduced to a Fredhol ntegral equaton of the second knd Keywords: Integral transfors, dual ntegral equatons, xed boundary Condtons, heat equaton 1 Introducton Integral transfor ethod s wdely used to solve several probles n heat transfer theory wth dfferent coordnate systes for unxed boundary condtons [1,8] In onographs [3-6] Hankel and Laplace transfors were effectvely used to nvestgate exact solutons for Helholtz and heat equaton subject to xed boundary condtons of the frst,the second and of the thrd knds

Naser Hoshan 16 for cylnder In ths paper we propose the soluton of three-densonal nonstatonary heat equaton n a non axally syetrcal cylndrcal coordnates wth dscontnuous xed boundary condtons of the second and thrd knd on the level surface of a se-nfnte sold cylnder Exact soluton of the gven xed boundary value proble s obtaned wth the use of fnte Fourer, Hankel ntegral transfors separaton of varables and based on the applcaton of dual ntegral equatons ethod In ths paper we apply a fnte Fourer ntegral transfor and then Hankel ntegral transfors wth respect to coordnate varables and r, oreover, an ntal xed boundary value proble s transfored to a Helholtz boundary proble n cylndrcal coordnates, next,applcaton of xed boundary condtons yelds new for of a dual ntegral equaton wth Bessel functon of the frst knd of order n as a kernel, weght and unknown functons depend on paraeters The soluton of the obtaned dual ntegral equatons s ntroduced to a Fredhol ntegral equaton of the second knd wth kernel and free ter gven n for of proper ntegrals 2 Forulaton of the Proble The an goal n ths paper s to solve the non-statonary heat equaton for se-space n cylndrcal coordnates wth a non-axally syetry (21) 1 1 1 2 r r r r z a T T T T T where T T ( r, z,, ) s the teperature dstrbuton functon, r, z, 2, are the correspondng cylndrcal coordnates a s the teperature dffusvty coeffcent (constant) The boundary condtons (22) T (, z,, ) T (, z,, ) T ( r,,, ), T () T (2 ) Subject to xed dscontnuous boundary condtons of the second knd and of the thrd knd (23) T ( r,,, ) / z f 1( r,, ), r S (2,4) T ( r,,, ) / z T ( r,,, ) f 2( r,, ), r S where S (, ), S (, ), on a surface z,, constants The ntal condton s (2,5) T ( r, z,,) ( r, z,,) U

17 Integral Transfors and Dual Integral Where U s the ntal teperature (constant), the known functons f, 1,2 n (23),(24) contnuous and have the lted varaton wth respect of each varables r and,wth respect to r, and perodc wth perod 2 The physcal sgnfcance of the gven xed boundary value proble forulated such that, fnd the teperature dstrbuton functon T for a senfnte cylnder f the nsde dsk S a heat flux s gven accordng to Fourer low, whereas on the outsde dsk S a heat exchange obey Newton's low of coolng, on the lne of dscontnuty r no boundary condtons were gven As r z the teperature s vanshed; T 2 2 3 Soluton of the Proble Use a fnte Fourer transfor for a functon T wth respect to, to the ntal boundary value proble (21)-(25), we have T 2 TF d, K() K(2 ), Where cos, j 2 Fj sn, j 2 1 The nverse Fourer ntegral transfor s gven by the forula 1 T T ( T 2 cos T 21 sn ) 2 1 Where 2 T T () Td, 2 2 1 21 j j 2 T T (2 ) T F d T T (2 1) T F d 2 2 j Equaton (21) n a Fourer transfor range s T 1 T T 1 T (31) T r r r r z a To splfy the proble entoned above,use well known transforaton [1] (32) T ( r, z,, t ) exp( ) u( r, z, ) for (21)-(25)after applcaton Fourer transfor, s constant, we obtan a Helholtz equaton n cylndrcal coordnates u 1 u u (33) u u r r r r z a Separatng varables n (33), the general soluton of the Helholtz boundary value proble s obtaned n for of proper ntegral

Naser Hoshan 18 (34) u ( r, z, ) A ( p, )exp( h( p, )) J ( pr ) dp 2 where h( p, ) p / a, J ( pr ) s the Bessel functon of the frst knd of order, p s the paraeter of separaton of varables, A( p, ) unknown functon Applyng a xed boundary condtons (23) and (24) to (34), we obtan the dual ntegral equatons to deterne the unknown functon A( p, ) (35) A( p, ) J ( pr) h( p, ) dp f 1( r, ), r S (36) A( p, ) J ( pr )( h( p, ) ) dp f 2( r, ), r S 2 (, ) f (, ) j f r e r F d To solve the dual equatons (35),(36), rewrte the equatons n the standard for (37) B ( p, ) J ( pr) g ( p, ) dp f 1( r, ), r S (38) B ( p, ) J ( pr ) dp f 2( r, ), r S where h( p, ) w ( p, ), l w ( p, ) 1/, h( p, ) p B ( p, ) ( h( p, ) ) A( p, ), g ( p, ) w ( p, ) 1/ h( p, ) ewrte (38) n for (39) ( r, ), r S B ( p, ) J ( pr ) dp f 2( r, ), r S where ( r, ) s unknown contnuous functon defned outsde the dsk n S Applyng to (39) the nverse Hankel ntegral transfor [7] n the nterval S S we have (31) B ( p, ) ypj ( py ) ( y, ) dy ypj ( py ) f 2( y, ) dy Substtutng (31) nto (37), then nterchangng the order of ntegraton, we get a Fredhol ntegral equaton of the second knd for deternaton the unknown functon ( r, ) (311) ( r, ) ( y, ) K ( r, y, ) dy F ( r, ), r S

19 Integral Transfors and Dual Integral wth kernel (312) and free ter K ( r, y, ) pyj ( py ) J ( pr ) g ( p, ) dp (313) F ( r, ) f 2( r, ) p y f 2( y, ) J ( pr) J ( py ) g ( p, ) dp dy Integral equaton (311) should be solved wth the use of nuercal ethods for soe choces of f 1, f 2,,,,1,2, by usng soe software packages such atheatca or atlab The kernel gven n (312) contnuous and quadratc ntegrable n the square : r, y, for soe certan nuercal values of furtherore, the free ter (313) s ntegrable and bounded n the nterval y [2]Fnally the general soluton T ( r, z,, ) n the Fourer transfor doan s gven by the expresson 1 T ( r, z,, ) e h( p, ) (314) ypj ( py ) ( y, ) dy ypj ( py ) f 2( y, ) dy h( p, ) J ( pr) dp Put the value of the general soluton (314) nto the nverson forula of the nverse Fourer transfor, the general soluton of the ntal xed boundary value proble (21)-(25) T ( r, z,, ) (315) 1 T ( r, z,,) ( T ( r, z,,2 )cos T ( r, z,,2 1)sn ) 2 1 If, the soluton (315) reduced to the soluton of an axal syetry heat equaton wth xed condtons[6],oreover, f h( p, ) p, the above soluton s reduced to the soluton of the Laplace's' equaton wth xed condtons n dfferent engneerng and physcal applcatons[3-7] Theory entoned above of the applcatons of a ntegral transfors nvolvng exact soluton of the xed ntal boundary value proble can be used wdely to solve varous xed boundary probles for a non-statonary heat equaton n an nfnte or fnte cylnder, unsyetrcal cylndrcal coordnates for unbounded plate, sphercal coordnates and other xed probles 4 Concluson Fnally the above technque nvolvng applcaton of ntegral transfors for solvng xed boundary value probles can be used for nvestgatng several

Naser Hoshan 2 hoogeneous probles (heat equaton, Helholtz equaton and Laplace equaton) n dfferent coordnate syste and any areas of applcatons n techncal and physcal scences under xed boundary condtons of the frst, the second and of the thrd knds 5 Open Proble We consder an ntal xed boundary value proble (21) (22) for an nfnte plate of hgh h, n cylndrcal coordnates wth non axally syetry, subject to dscontnuous nhoogeneous xed boundary condtons of the second and of the thrd knd on a surface z (316) T ( r,,, ) / z f 1( r,, ), r S, (317) T ( r,,, ) / z T ( r,,, ) f 2( r,, ), r S On a level surface boundary condtons (318) T( r, h,, ) / z T ( r, h,, ) f3( r,, ) z h, located a thrd knd lnear nhoogeneous Where f, 1, 2,3 known functons, S (, ), S (, ),,,,, constants The above xed proble (316)-(318), ntroduced to soe type of dual ntegral equatons, however no one n the world solve ths proble, snce the boundary condton (318) coplcates soluton of the gven proble furtherore, known ethods concernng dual ntegral equatons ay be dffcult to use eferences [1] A Galtsyn, A Zhukovsk, Integral Transfors and Specal Functons n Heat Probles, Kev, Duka1976 [2] W Hackbusch, Integral Equaton, Theory and Nuercal Treatent Brkhäuser Verlag, Boston,1995 [3] N Hoshan, Dual Integral Equatons and Sngular Integral Equatons for Helholtz Equaton Internatonal Journal of Contep Math Scences, 29, V4, No 34,1695-1699,Hkar Ltd [4] N Hoshan, Integral Transfor Method n Soe Mxed Probles Internatonal Journal of Matheatcal Foru, 29, V4, No 4, 1977-198, Hkar Ltd

21 Integral Transfors and Dual Integral [5] N Hoshan, The dual ntegral equatons ethod nvolvng heat equaton wth xed boundary condtons, Engneerng Matheatcs Letters, No2 213, pp 137-142 [6] N Hoshan, The Dual Integral Equatons Method for Solvng Hlholtz Mxed Boundary Value Proble, Aercan Journal of Coputatonal and Appled Matheatcs, 213, 3(2), pp 138-142 [7] B Mandal, N Mandal, Advances n Dual Integral Equaton, London, CC1999 [8] M Ozsk, Heat Conducton, Wley & Sons NewYork, 22