COMPUTATIONAL FLUID DYNAMICS - Help with SIMPLE algorithm implementation

I have the following program (by some standards, it might be called
long). I want any advice as how to better understand it. For those in
the know, it is meant to implement the SIMPLE algorithm for horizontal
channel flow. It uses a mesh/grid system for this. I will try to
offer explanations to anyone who is interested in offering help.

Here's the program. I would very much like to know how to modify it
for horizontal and vertical fluid flow, and how to deal with a channel
that looks like what I am about to draw (the third channel):

Horizontal Channel:
<---------------------10m----------------------->
-------------------------------------------------/\

1m

-------------------------------------------------\/

Boundary Conditions/Initial Conditions.

On the LHS of the channel, the velocity is given by:
u=6(0.25-y^2)
v=0

On the walls of the channel, the velocity is given by:
u=v=0

At the RHS of the channel, the velocity derivative is given by:
du/dx=0, dv/dx=0

The pressure is given by p = -0.12x

Vertical Channel:

<1m>
| | /\
| |
| |
| |
| | 10m
| |
| |
| |
| | \/


Here, the outlet is the top part and we have the analogous boundary
conditions.


Channel With 'Bump':



/\-------------------------------------------------------------------------


1m


\/---------------------
/\ |
|
1m |
|
\/
|-----------------------------------------------------
< 2m > < 15m
>

Here, we have the same boundary conditions (inlet, outlet and wall
conditions).

I know what the is given by roughly, but would like to know what
adaptations it would be wise to make to the below program to get it to
work.......ESPECIALLY for the third case.

Thanks for any help!

.

Relevant Pages

  • Re: OT: G-20
    ... But I think a six outlet device will fit behind my stereo/TV much ... wire retaining bales for standard IEC plugs. ... Throw that two channel crap away and get a real "stereo"! ...
    (sci.electronics.design)
  • Re: Help stop dog fighting
    ... also another outlet to plug earphones into.>> ... for the left channel, and one for the right channel. ... be white & red for audio, ...
    (soc.retirement)
  • Re: Velocity measurement
    ... quadrature, which measures the time between transition events such as ... channel A rising edges. ... If I were to use PIC for measuring the velocity, ... an incline and require some simple geometry to correct. ...
    (comp.arch.embedded)
  • Re: Role of velocity in note_off event ?
    ... ('note_off', dtime, channel, note, velocity) ... so that future developments could be incorporated easily. ...
    (comp.music.midi)
  • Re: Role of velocity in note_off event ?
    ... ('note_off', dtime, channel, note, velocity) ... how do synthesisers interpret "velocity" in a ...
    (comp.music.midi)