|
[Sponsors] |
March 6, 2013, 15:13 |
Patch to Patch mapping
|
#1 |
Member
Ronald McDonald
Join Date: Jul 2012
Posts: 38
Rep Power: 14 |
Hi OpenFoamers,
I have a really simple question and desire a simple, if possible, response. I have two boxes connected to each other. I have labelled those adjacent boundaries patch1 and patch2. I have one mesh throughout the whole domain. I want to be able to have two fields throughout the domain, where one ends at patch 1 and the values are transferred onto another field via patch 2. Is this possible in openFoam, if so how can I go about doing this? Best, Benjamin |
|
March 6, 2013, 17:32 |
|
#2 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 |
Hi Benjamin,
You could have told me in the PM you sent me that you had asked this question here on this thread Anyway, for other readers, the core of the answer I gave in the private message was: Quote:
Bruno
__________________
|
||
March 6, 2013, 19:40 |
Thanks, Bruno
|
#3 | |
Member
Ronald McDonald
Join Date: Jul 2012
Posts: 38
Rep Power: 14 |
Quote:
So I think I got the patch to patch mapping working with groovybc. The only thing is I'm getting a light warning message. HTML Code:
Time = 0.145 --> FOAM Warning : From function ExpressionResult::getUniformInternal(const label size,bool noWarn) in file ExpressionResult/ExpressionResultI.H at line 322 The minimum value 999.96 and the maximum 999.96 differ. I will use the average 999.96 DICPCG: Solving for H2O, Initial residual = 8.52782e-07, Final residual = 8.52782e-07, No Iterations 0 DICPCG: Solving for H2, Initial residual = 8.11531e-07, Final residual = 8.11531e-07, No Iterations 0 --> FOAM Warning : From function ExpressionResult::getUniformInternal(const label size,bool noWarn) in file ExpressionResult/ExpressionResultI.H at line 322 The minimum value 999.96 and the maximum 999.96 differ. I will use the average 999.96 DICPCG: Solving for H2O, Initial residual = 8.52782e-07, Final residual = 8.52782e-07, No Iterations 0 DICPCG: Solving for H2, Initial residual = 8.11531e-07, Final residual = 8.11531e-07, No Iterations 0 Benjamin |
||
March 7, 2013, 08:17 |
|
#4 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 |
Hi Benjamin,
That's one strange warning... there must be something beyond "999.96"... such as "99.96324231235"... but I guess it's perfectly OK. Best regards, Bruno
__________________
|
|
March 7, 2013, 11:51 |
|
#5 |
Member
Ronald McDonald
Join Date: Jul 2012
Posts: 38
Rep Power: 14 |
Hi Bruno,
So I ran my simulation with no boundary patches (the ones I were mapping) and found that my results there are linear. Which is theoretically correct. However, when I try to use groovyBC to combine the two patches I get linear up to the coupled patches then a parabolic after the patch, which is theoretically incorrect. Any ideas on how to solve this? I heard that pointtopointinterpolation.H is another way to go. I have no idea how to use the header file, however. Any thoughts, and/or direction would be greatly helpful. Benjamin |
|
March 7, 2013, 16:38 |
|
#6 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 |
Hi Benjamin,
I'm confused... didn't you want to map one field to another different field? Something like mapping "T" to "p" (temperature to pressure) or something like that? Or do you want to map a single field, e.g. "T", from one patch to another patch? Best regards, Bruno
__________________
|
|
March 7, 2013, 19:06 |
|
#7 | |
Member
Ronald McDonald
Join Date: Jul 2012
Posts: 38
Rep Power: 14 |
Quote:
At first I wanted to do the former, but now I want to do the latter, namely map within a single field one patch to another patch. I've been trying to use patchtopatchinterpolation.H but i haven't been able to get it to work. And groovyBC did work (there's an example of it in the link you provided) but the physics came out incorrectly. Any help would be great. Sincerely, Benjamin |
||
March 8, 2013, 18:18 |
|
#8 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 |
Hi Benjamin,
Right now I don't have time to do any tests, but here is a quick hint - run this command: Code:
grep -R map $FOAM_TUTORIALS Good luck! Bruno
__________________
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
y+ and u+ values with low-Re RANS turbulence models: utility + testcase | florian_krause | OpenFOAM | 114 | August 23, 2023 05:37 |
[GAMBIT] periodic faces not matching | Aadhavan | ANSYS Meshing & Geometry | 6 | August 31, 2013 11:25 |
CheckMeshbs errors | ivanyao | OpenFOAM Running, Solving & CFD | 2 | March 11, 2009 02:34 |
[Gmsh] Import gmsh msh to Foam | adorean | OpenFOAM Meshing & Mesh Conversion | 24 | April 27, 2005 08:19 |
Multicomponent fluid | Andrea | CFX | 2 | October 11, 2004 05:12 |