how to set a block with phoenics3.4
 

I have been working hard to the numerical simulation of combustion with BFC. Lately, to my painful, I cann¡¯t set a block in the above of nozzle inlet. When I set a block, the iterations(sweeps) end with only 2 sweeps. However,If I substitute 4 plates for block, the solving process is very satisfied. I don¡¯t know why this is? In fact, the subsititute is not accord with the physical model. How to set a block with BFC.? How to modify Q1file?(considering vpor,epor,hpor etc.)

Re: how to set a block with phoenics3.4
 

zhujianguo,

I am not sure if this will answer your question, since I run v3.2 (also with BFC grid), but let's try.

It is recommended to use solid region rather than VPOR etc. therefore you should first add in GROUP 7:

store(prps)

Now you set IC in GROUP 11 - initialize prps and then set the blocked region with a solid material:

fiinit(prps) = 0.

patch(block, inival, ixf, ixl, iyf, iyl, izf, izl, 1, 1) ; init(block, prps, 0., 198.)

That should be all. For any blocked region surface that is a no-slip wall, you should add WALL command in GROUP 13, e.g.:

wall (esurf, east, ixfe, ixle, iyfe, iyle, izfe, izle, 1, Lstep)

However, due to a bug in v3.2 - "wall" command of the boundaries of blocked regions is treated as IC patch and should therefore replaced by "patch" and "coval", manually adding all affected variables and especially adding the parallel velocity component and the turbulence solved variables (if any). For the above wall, this should be:

patch(esurf, ewall, ixfe, ixle, iyfe, iyle, izfe, izle, 1, Lstep) ; coval(esurf, v1, grnd2, 0.) ; coval(esurf, ke, grnd2, grnd2) ; coval(esurf, ep, grnd2, grnd2) ; ...

I hope it helps, Rami