
[Sponsors] 
May 2, 2002, 17:03 
Multiblock

#1 
Guest
Posts: n/a

Is it possible to treat the multi block using an iterative manner? Just like feeding the outflow of one block to the other block as an inflow. If it is possible, what is the boundary conditions on the interblock surface? Thanks


May 3, 2002, 03:26 
Re: Multiblock

#2 
Guest
Posts: n/a

in the multoblock technique, two layer cells are extended out around every block. They are named as virtual cells and dummy cells. They are correspondent with the cells in the neigbouring blocks, which means the virtual cells and dummy cells have the same flow field parameters with the correspondent cells in other block. In the iterative process, the boundary conditions on the interblock surfaces can therefore be transfered with each other.


May 5, 2002, 17:39 
Re: Multiblock

#3 
Guest
Posts: n/a

Thanks for your help. But can you give me more details . For example, for a very simple case, if I devide a backstep into two zones (though this is some kind of unnecessary to use multiblock), and use dummy cells as you said. When I solve for the first zone, what is the boundary condition that I can use for the dummy cells, outflow? And if using an iteritive method, I can use that outflow from zone1 for my second zone and solve the zone2. But what is the next step, how the iterative way continue? Thanks very much.


May 8, 2002, 03:24 
Re: Multiblock

#4 
Guest
Posts: n/a

When you devide a backstep into two zones, the two zones have a common interface. The dummy cells and the virtual cells should be developed around the whole block, including the common interface and other boundaries. At the beginning of the calculation, the flow field in the solution domain should be initialized. The flow field distribiton (or boundary condition) for the dummy cells in the first clok on the common interface should be valued by the distribution in the corresponding cells in the block 2. The flow field distribution for the dummy cells in the first clock on the other boundary are taken from normal boundary conditions. For the block 2, it is same. After the first iteration, revalue the flow field disribution for the dummy cells in the first block on the common interface according to the flow field you obtains in the block 2. For the block 2, it is same. It means you should exchange the data for the dummy cells on the interface in each iteration.
You can also find some documents in the internet which introduce how to exchang data in the multi blocks and paralle calculation. Hope it will be helpful! 

May 9, 2002, 09:17 
Re: Multiblock

#5 
Guest
Posts: n/a

Has anyone tried the following procedure?
I am assuming a finite element method using sparce matrix techniques. Assume the blocks have been independently assembled into a global matrix. This global matrix will be block diagonal. Matrix blocks corresponding to two ajoining multiblocks will each contain partial contributions to degrees of freedom on common the edges. Assume one has a fast way to solve the individual uncoupled blocks, say by a parallel computation with a separate processor handling each block. The blocks can be coupled by replacing one column of the set of common DOFs (a pair if the DOF is on an edge, more if on a common corner of several blocks) by the sum of the columns from different blocks and the rows by the sum of the rows from different blocks. The remaining rows and columns are replaced by zeros except for a 1 on the diagonals and a 1 in the summed column (a simple change can make this symmetrical). The RHS would be modified as well. This procedure completes the assembly of one of the common DOFs and constrains the remaining common DOFs to be equal. This modification of the global matrix is of the form so that the Woodbury formula (a generalization of the ShermanMorrison formula for updating the solution or inverse of a matrix if one matrix element is changed see "Numerical Recipes in C/C++/Fortran/etc") can be used to find the (correct) coupled solution from uncoupled solutions using the fast (parallel) block solvers. I suppose this could be used iteratively, turning on and off the coupling between pairs of blocks. I would appreciate a thoughtful critique or references if anyone has seen this done. 

May 12, 2002, 22:25 
Re: Multiblock

#6 
Guest
Posts: n/a

Can I do it in this way: I get the velocity at the interface, and use it as inflow for the second zone. After generating the pressure distribution in second zone, I put the pressure on the common boundary and compute the flow in first zone again, so iteration keeps on. It seems ok since the result is very close to that of FLUENT, but is there problem for such scheme? Thanks very much.


May 15, 2002, 11:02 
Re: Multiblock

#7 
Guest
Posts: n/a

It seems Ok to exhange the data in this way. If the data contained in the boundary conditions are enough for the calculation in the second block, it should work too.


May 15, 2002, 11:12 
Re: Multiblock

#8 
Guest
Posts: n/a

Sorry,I forget one point. With your method, as the flowfield on the interface in the first block is taken as the inflow of the second block, the corresponding cells near the interfaces in the two blocks have the same flowfield. Theoretically, it is not right. Since continuity property is expected on the interface, not the equal flow field. But you can taken it as a coarse numerical algorithm, no problem.


May 20, 2002, 16:51 
Re: Multiblock

#9 
Guest
Posts: n/a

Thanks. Since I am using finite volume method and matched grid, the same velocity distribution can satisfy continuity, do you think so? And I am not clear some points in your reply in 8th, May. 1. Why dummy cells are all around the domain, are they computed besides the cells around interfaces? 2. How do you start compute zone2 for the back step case? Because the inflow face for this zone is the interface, how do you decide the inflow for this zone? And can you recommend me a more detailed paper for such kind of multiblock manipulations? Thanks very much


May 27, 2002, 05:44 
Re: Multiblock

#10 
Guest
Posts: n/a

It is not so easy to explain it only with text. And I don't have a right paper just on hand. Please search in the Internet with the key words "multiblock" or "Message passing library". You will find enough useful information.
Good Luck! Jing 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
coupling and multi block meshing  abdullahkarimi  CDadapco  1  November 11, 2010 03:20 
Icem Mesh to Foam  jphandrigan  OpenFOAM Mesh Utilities  4  March 9, 2010 03:58 
Multi block geometry  Samir  FLUENT  0  January 14, 2005 16:55 
Multi Block case.  Constantinos Constantinou  FLUENT  1  October 26, 2001 02:00 
F90, modules and multi block  Frank Muldoon  Main CFD Forum  0  October 13, 1999 22:37 