Blockage in CFD code
I am currently playing with the CFD code in that comes with the book by Peric and Ferziger.
I would very much like to implement a method for setting blockage internally in the domain.
Can any one give me a hint on how to do this ?
What is the best way ?
Thanks in advance.
Re: Blockage in CFD code
I am not sure how the code u are referring to is written. but the general way is (Patankar's funda). for any grid point u write the equation as
Ap_phi_p = Sigma A_i Phi_i +Su
Where Ap= Sigma A_i+Sp
Phi is a general variable (u,v,w,t,p,k,e ....)
Make Sp very large (say 1e25) and similarly Su= 1e25 (must be equal to value of Sp) for the grid point for which u want to freeze the value. U can initialise the value required and it will remain unchanged throuthout the iterations. if u put Phi_p =0 for u,v,w u will create islands in the flow path. Specify these values for all the grid-points that u want to freeze.
Note: for equations having large value of variables this may blow up due to overflow (typically k-e) Some care is required in grid spacing at blocked areas.
Hope this helps
|All times are GMT -4. The time now is 23:18.|