CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Multi-block

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   May 2, 2002, 17:03
Default Multi-block
  #1
Arthur
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 inter-block surface? Thanks
  Reply With Quote

Old   May 3, 2002, 03:26
Default Re: Multi-block
  #2
jing
Guest
 
Posts: n/a
in the multo-block 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 inter-block surfaces can therefore be transfered with each other.

  Reply With Quote

Old   May 5, 2002, 17:39
Default Re: Multi-block
  #3
Arthur
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 multi-block), 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.
  Reply With Quote

Old   May 8, 2002, 03:24
Default Re: Multi-block
  #4
jing
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!
  Reply With Quote

Old   May 9, 2002, 09:17
Default Re: Multi-block
  #5
Jonas Holdeman
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 multi-blocks 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 Sherman-Morrison 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.

  Reply With Quote

Old   May 12, 2002, 22:25
Default Re: Multi-block
  #6
Arthur
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.
  Reply With Quote

Old   May 15, 2002, 11:02
Default Re: Multi-block
  #7
Jing
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.

  Reply With Quote

Old   May 15, 2002, 11:12
Default Re: Multi-block
  #8
Jing
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.
  Reply With Quote

Old   May 20, 2002, 16:51
Default Re: Multi-block
  #9
Arthur
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
  Reply With Quote

Old   May 27, 2002, 05:44
Default Re: Multi-block
  #10
Jing
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 "multi-block" or "Message passing library". You will find enough useful information.

Good Luck!

Jing
  Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
coupling and multi block meshing abdullahkarimi CD-adapco 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


All times are GMT -4. The time now is 23:53.