# Blottner Coefficients

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

August 26, 2010, 07:03
Blottner Coefficients
#1
New Member

Freddie Sizer
Join Date: Apr 2009
Posts: 19
Rep Power: 10
I have been working with a CFD code that requires Blottner coefficients as an input. These are not easy to come by, but I've written a simple program in Scilab that calculates the coefficients for you.

The user inputs required are three temperatures and their corresponding viscosities, which are available on the NIST website.

http://webbook.nist.gov/chemistry/fluid/

The output is given in the x vector as the Blottner Coefficients A,B and C. (x(1)=A, x(2)=B etc.)

As an example I have included data for Helium at 1bar.

//Blottner
T1=100;
T2=800;
T3=1500;
u1=6.9644e-06;
u2=3.5932e-05;
u3=5.4056e-05;
//Calculating Blottner viscosity coefficients for He at 1bar
A=[(log(T1))^2 log(T1) 1 ; (log(T2))^2 log(T2) 1; (log(T3))^2 log(T3) 1];
u=[log(u1)-log(0.1); log(u2)-log(0.1); log(u3)-log(0.1)];
x=lusolve(sparse(A),u);
T=[10:10:1500]';
mu=0.1*exp((log(T)*x(1)+x(2)).*log(T)+x(3));
//Calculating Sutherland Viscosity
eta = ((mu0*(T0+C)).*((T./T0)^(3/2)))./(T+C);
//Import Helium data from NIST file
HEV=fscanfMat('HeliumNISTonebar.txt');
plot(T,mu,HEV(:,1),HEV(:,2))

The attached .txt file was taken from the NIST site by simply copying the html output option into excel, deleting not needed columns, and saving as .txt.

Matlab users will need to exchange the lusolve function for the gaussian elimination equivalent in Matlab.

I hope it helps someone!
Attached Files
 HeliumNISTonebar.txt (3.6 KB, 4 views)

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Xwang OpenFOAM 20 May 20, 2016 11:26 omurakbayir FLUENT 1 November 22, 2010 23:36 Aragon FLUENT 0 July 22, 2009 04:07 Brian FLUENT 0 September 8, 2004 09:04 Michael Main CFD Forum 2 December 8, 2003 13:15