CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

The Blasius equation in Python

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 3, 2011, 15:49
Default The Blasius equation in Python
  #1
New Member
 
Join Date: May 2011
Posts: 1
Rep Power: 0
javiergra24 is on a distinguished road
Hi everybody

I'm writing a script in Python to solve the Blasius equation but it does not work, numerical results does not match with data I've seen in fluid mechanics books.

Please anyone could help me? Thank you very much

Script is

import sys, pylab, numpy
from pylab import *
from numpy import *
import matplotlib.pyplot as plt


deta=0.0001
total=10


e = []
f = []
g = []
h = []


fvec=0; gvec=0; hvec=0.3219

etavec=0.0

while etavec <= 10:

fvec = fvec + gvec * deta
gvec = gvec + hvec * deta
hvec = hvec -(1/2) * fvec * hvec * deta
etavec = etavec + deta
e.append(etavec)
f.append(fvec)
g.append(gvec)
h.append(hvec)
print etavec,fvec,gvec

numpy.savetxt("blasius.dat", transpose((e,f,g,h)),delimiter=' ')

pylab.figure()
pylab.plot(e,g)
pylab.plot(e,h)
xticklines = getp(gca(), 'xticklines')
yticklines = getp(gca(), 'yticklines')
xgridlines = getp(gca(), 'xgridlines')
ygridlines = getp(gca(), 'ygridlines')
xticklabels = getp(gca(), 'xticklabels')
yticklabels = getp(gca(), 'yticklabels')
#lines1 = pylab.plot(e,f)
lines2 = pylab.plot(e,g)
lines3 = pylab.plot(e,h)
#pylab.setp(lines1, color='b', linewidth=3.0)
pylab.setp(lines2, color='r', linewidth=3.0)
pylab.setp(lines3, color='k', linewidth=3.0)
setp(xticklines, 'linewidth', 3)
setp(yticklines, 'linewidth', 3)
setp(xgridlines, 'linestyle', '--')
setp(ygridlines, 'linestyle', '--')
setp(yticklabels, fontsize='xx-large')
setp(xticklabels, fontsize='xx-large')
xlabel('$\eta$',fontsize=30)
ylabel('$f(\eta) \quad \partial_{\eta} f(\eta) \quad \partial_{\eta \eta}f(\eta)$',fontsize=30)
axis([0.0,10,-0.9,20.9])
#legend( (lines1, lines2, lines3), ('$\f(\eta)$', '$\partial_{\eta} f(\eta)$', '$\partial_{\eta \eta}f(\eta)$'))


grid(True)
pylab.show()
javiergra24 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 06:21
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 07:27
blasius equation varunjain89 Main CFD Forum 15 June 25, 2012 10:42
Blasius equation help :( johnt447 Main CFD Forum 0 March 18, 2011 18:57
Constant velocity of the material Sas CFX 15 July 13, 2010 09:56


All times are GMT -4. The time now is 19:17.