[atmcfd]

i`m going through anderson`s 'CFD and its applications",and i`m stuck on pg. 308
(analysis of quasi one dimensional subsonic-supersonic nozzle). here the initial condition of density,temperature and velocity is given by the equations

p = 1 - 0.3146x
T = 1 - 0.2314x . initial conditions at t = 0
V = (0.1 + 1.09x)T^1/2
where 'x' is the axis where the grid points are evenly spaced

forgive me if dis question`s stupid, but what is the origin of these eqns? there is no explanation offered in the text as to what are the governing principles behind these. are these just taken arbitrarily or what...
please help.

thanks in advance,
atm

[gocarts]

The subsonic-supersonic nozzle appears in Modern Compressible Flow, 2nd Edition, Anderson, pp. 365 (The Time-Marching Technique) with the following quote:

"...the guessed initial temperature distribution at time t=0. It is arbitrarily taken as linear variation."

I presume it's likely the same example in your book "CFD and its Application."

So I believe all the initial values are based on knowing the solution, so they have no theoretical foundation. They are probably chosen to be just close enough to the final (known) solution to make your solver work, but not far enough away to cause problems.

If you are following along and writing your own solver I would suggest substituting fixed values and see what converged solution you get.

[atmcfd]

thank u very much...
is it better if i write my computer code in C++ or MATLAB?
how do the 2 differ?
thanks in advance.
atm

[RenardP]

Since your problem is 1D you can use C++ or Matlab. For 2D and 3D problems Matlab is simply to slow.

Use the language you are more familiar C++ or Matlab for this simple problem.

[atmcfd]

oh... thank u very much.i guess C++ is better.