پروژه کدنویسی مقاله با فرترن

on MHD Flow Near a Stagnation Point on a Linearly Stretching

Sheet

چکیده گزارش

در این گزارش ابتدا مساله و معادلات حاکم بر آن تشریح شده و سپس با استفاده از حل تشابهی این معادلات دیفرانسیل جزئی به سری معادلات دیفرانسیل معمولی تبدیل شده است. با استفاده از روش رانج کوتای مرتبه چهار و شوتینگ غیر خطی این دستگاه معادلات حل شده­اند.

با استفاده از روش تشابهی معادلات فوق به معادلات دیفرانسیل معمولی تبدیل شده­اند.

روش حل:

طبق معادلات فوق با دستگاه معادلات دیفرانسیل غیر خطی سر و کار داریم اما برای حل به مقادیر  نیاز داریم. برای به دست آوردن آن­ها از روش شوتینگ غیر خطی استفاده می­کنیم. ابتدا برای آن­ها  یک مقدار اولیه حدس می­زنیم ، سپس با استفاده از روش شوتینگ نیوتن مقادیر را اصلاح می­کنیم.

در ضمن با نرم افزار گیت دیتا و نمودارهای مقاله و در نظر گرفتن  چند نقطه ،نمودار را با EXCELL رسم کرده و با نمودار مقاله کنار هم به نمایش گذاشته شده است.

function p1(j1) result(pj1)

real::pj1

real,intent(in)::j1

pj1=j1

end function p1

!!!!!!!!!!!!!!!!

function p2(j2) result(pj2)

real::pj2

real,intent(in)::j2

pj2=j2

end function p2

!!!!!!!!!!!!!!!!

function p3(j3,j4,j5) result(pj)

real::pj,m,la

real,intent(in)::j3,j4,j5

m=0.0;la=2.0

pj=-j3*j4+j5*j5+m**2*(j5-la)-la*la

end function p3

!!!!!!!!!!!!!!!!!

function p4(jjj4) result(pj4)

real::pj4

real,intent(in)::jjj4

pj4=jjj4

end function p4

!!!!!!!!!!!!!!!!

function p5(j8,j9,j10) result(p5j)

implicit none

real::p5j,e,S,Pr

real,intent(in)::j8,j9,j10

e=0;S=0;Pr=0.01

p5j=(-e*j8*j8-pr*j8*j9-pr*s*j10)/(1+e*j10)

end function p5

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!1shooting function

function q1(v) result(qv1)

real::qv1

real,intent(in)::v

qv1=-v

end function q1

!!!!!

function q2(v2) result(qv2)

......

جهت دانلود اینجا کلیک کنید