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

Flow in a Diverging Nozzle

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 25, 2013, 11:47
Default Flow in a Diverging Nozzle
  #1
New Member
 
Join Date: Jan 2013
Posts: 17
Rep Power: 4
kamyar is on a distinguished road
Hi
I'm Writing a code in matlab to solve the euler equation for flow in a diverging nozzle with supersonic inlet and subsonic outlet (sectioan 12-5 of Computatioanl fluid dynamics Vol.2 Hoffman). I'm using Steger and warming Flux Vector Splitting explicit formulation. the code doesn't work. Here is the code:
% Allocating variables:
dx = 0.1;
dt = 0.00001;
x = 0:dx:1;
s = 1.398 + 0.374 * tanh(0.8*x - 4);
ds = 0.2776 * ( 1 - tanh(0.8*x - 4).^2 );
gama = 1.4;
n = size(x,2);
ro = zeros(1,n);
u = zeros(1,n);
P = zeros(1,n);
et = zeros(1,n);
Q = zeros(1,n);
H = zeros(1,n);
E_plus = zeros(3,n);
E_minus = zeros(3,n);
DELP = 1;
Iteration = 0;
% Setting initial and boundary conditions for Supersonic inlet and
% subsonic outlet:
u(1:29) = 1676.55;
u(30:n) = 572.76;
ro( : ) = 0.002241;
P( : ) = 2000;
et= 0.5 * u.^2 + P ./ ( ro .* (gama-1) );
a = sqrt(gama.*P./ro);
Q(1,: ) = ro;
Q(2,: ) = ro .* u;
Q(3,: )= ro .* et;
% Sloving the euler equation using Steger and Warming Flux Vector
% Splitting:
while DELP > 0.1
P0 = P;
H(2,: ) = ds .* P;
E_plus(1,: ) = s .* ro .* ( 2 * gama * u + a - u ) / 2 / gama;
E_plus(2,: ) = s .* ro .* ( ( 2 * ( gama - 1 ) * u.^2 ) + ( u + a ).^2 ) / 2 / gama;
E_plus(3,: ) = s .* ro .* ( ( ( gama - 1 ) * u.^3 ) + ( ( u + a ).^3 / 2 )...
+ ( ( ( 3 - gama ) .* ( u + a ) .* a.^2 ) / ( 2 * (gama - 1) ) )) / 2 / gama;
E_minus(1,: ) = s .* ro .* ( u - a ) / 2 / gama;
E_minus(2,: ) = s .* ro .* ( ( u - a ).^2 ) / 2 / gama;
E_minus(3,: ) = s .* ro .* ( ( ( ( u - a ).^3) / 2 ) + ...
( ( 3 - gama ) * ( u - a ) .* a.^2 ) / ( 2 * ( gama - 1 ) ) ) / 2 / gama;
for ii = 2:n-1
if u(ii) > a(ii)
E_plus(:,ii-1:ii) = E_plus(:,ii-1:ii) + E_minus(:,ii-1:ii);
E_minus(:,ii:ii+1) = 0;
end
Q(:,ii) = Q(:,ii) + dt * ( H(:,ii) - ( E_plus(:,ii) - E_plus(:,ii-1)...
+ E_minus(:,ii+1) - E_minus(:,ii) ) / dx ) / s(ii);
end
Q(1,n) = Q(1,n-1);
Q(3,n) = Q(3,n-1);
Q(2,n) = u(n) * ro(n);
ro = Q(1,: );
u( : ) = Q(2,: ) ./ Q(1,: );
et( : ) = Q(3,: ) ./ Q(1,: );
P = ro .* ( gama - 1 ) .* ( et - 0.5 * u.^2);
a = sqrt(gama * P./ro);
DELP = sum( abs(P - P0) );
Iteration = Iteration + 1;
end

What is wrong?
thanks.
kamyar is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Supersonic Inlet for a Compressible Flow in a Nozzle Daco FLUENT 8 November 8, 2012 04:34
help - Running Converging Diverging Nozzle with a shock wave Daco FLUENT 0 September 3, 2012 08:23
Review: Reversed flow CRT FLUENT 0 July 20, 2012 13:03
3D-Simulation of Water Flow From Nozzle into Air Navin Sampath FLOW-3D 2 February 26, 2009 11:46
Inviscid Drag at subsonic, subcritical Mach # Axel Rohde Main CFD Forum 1 November 19, 2001 13:19


All times are GMT -4. The time now is 13:25.