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

splitMeshRegions

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 3, 2010, 06:10
Default splitMeshRegions
  #1
VdG
New Member
 
Vincent de Graaf
Join Date: Jun 2010
Location: Germany
Posts: 12
Rep Power: 7
VdG is on a distinguished road
Hi all,

when running splitMeshRegions -cellzones in OF16x, boundaries of type "calculated" are not treated correctly. The values from the original solution are not mapped onto the split solution. I've attached an example based on the angledDuctImplicit tutorial. The original nut file (/10/nut) contains:

outlet
{
type calculated;
value nonuniform List<scalar>
400
(
0.000908662
0.000880609
....
)
;
}

In the split case, this becomes (in /11/outlet/nut):

outlet
{
type calculated;
value uniform 0;
}

A similar problem occurs for boundaries with nutWallFunction, kqRWallFunction etc.
Attached Files
File Type: zip splitMeshRegionsBugReport.zip (12.1 KB, 5 views)
VdG is offline   Reply With Quote

Old   June 3, 2010, 09:15
Default
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
splitMeshRegions was zeroing out all 'calculated' fvPatch fields since it assumed they were all on inter-region patches. I've pushed a fix to 1.6.x which only zeroes out the inter-region patches it has added. Can you let me know if it does the right thing?

(you will still have 'calculated' instead of nutWallFunction etc. since splitMeshRegions does not know about these boundary conditions. You can add a libs entry to your system/controlDict to force loading the libraries containing these boundary conditions)

Thanks for the nice bug report.
mattijs is offline   Reply With Quote

Old   June 3, 2010, 11:14
Default
  #3
VdG
New Member
 
Vincent de Graaf
Join Date: Jun 2010
Location: Germany
Posts: 12
Rep Power: 7
VdG is on a distinguished road
Hi Mattijs,

it works perfectly. Thanks!
VdG is offline   Reply With Quote

Old   June 9, 2010, 09:52
Default
  #4
VdG
New Member
 
Vincent de Graaf
Join Date: Jun 2010
Location: Germany
Posts: 12
Rep Power: 7
VdG is on a distinguished road
Another question: Why would you zero out the values on the inter-region patches? In my case, I need the values on these patches. Is it possible to write the values as a fixedValue, nonuniform BC?
VdG is offline   Reply With Quote

Old   June 9, 2010, 11:04
Default
  #5
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
The faces on these patches originate from internal faces. Hence there are no volFields on them (only possibly surfaceFields). You would have to do interpolation to get values on them. Usually these patches are to be used for coupling to they will have calculated values on them or some simple bc.

In your case maybe you can use mapFields to map a value onto them.
mattijs is offline   Reply With Quote

Old   June 14, 2010, 04:31
Default
  #6
VdG
New Member
 
Vincent de Graaf
Join Date: Jun 2010
Location: Germany
Posts: 12
Rep Power: 7
VdG is on a distinguished road
In my case, mapFields doesn't work, since I have insufficient memory. However, for now I can solve the problem by changing calculated to zeroGradient.
VdG is offline   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
splitMeshRegions naltang OpenFOAM 3 May 6, 2010 12:22
splitMeshRegions doesn't work. PetSul OpenFOAM Bugs 4 September 2, 2009 07:06


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