
[Sponsors] 
[swak4Foam] funkySetFields: problem with processor boundary 

LinkBack  Thread Tools  Search this Thread  Display Modes 
May 14, 2015, 16:40 
funkySetFields: problem with processor boundary

#1 
New Member
Join Date: May 2015
Location: Russia
Posts: 3
Rep Power: 9 
Hi FOAMers!
I use funkySetFields (swakVersion: 0.3.1) for the decomposed case. Internal fields are set remarkably. But I am disturbed by the following problem with processor boundaries. For example, one can consider shockTube case: for x \in [5;0) T=348.432, p=1e5, U=0 for x \in [0;5] T=278.746, p=1e4, U=0 Mesh decomposition is done for 4 subdomains: proc0: [5;2.5] proc1: [2.5;0] proc2: [0;2.5] proc3: [2.5;5] funkySetFieldsDict contains expressions ( leftT { field T; expression "348.432"; condition "pos().x<0"; keepPatches true; } rightT { field T; expression "278.746"; condition "pos().x>=0"; keepPatches true; } leftP { field p; expression "100000"; condition "pos().x<0"; keepPatches true; } rightP { field p; expression "10000"; condition "pos().x>=0"; keepPatches true; } ); Run funkySetFields in parallel: mpirun np 4 funkySetFields time 0 parallel The execution is successful. There are no warnings and errors. So, we analyze result now. To start, processor0/0/T: FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 348.432; boundaryField { sides { type zeroGradient; } empty { type empty; } procBoundary0to1 { type processor; value uniform 278.746; } } But T is still 348.432 in subdomain 1... Similar, processor0/0/P: FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 1 2 0 0 0 0]; internalField uniform 100000; boundaryField { sides { type zeroGradient; } empty { type empty; } procBoundary0to1 { type processor; value uniform 10000; } } P.procBoundary0to1 must be 100000. Further, processor1/0/T: FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField nonuniform List<scalar> 3(348.432 348.432 278.746); boundaryField { sides { type zeroGradient; } empty { type empty; } procBoundary1to0 { type processor; value uniform 278.746; } procBoundary1to2 { type processor; value uniform 278.746; } } T.procBoundary1to0 must be 348.432. I hope, the problem is explained. Does somebody know how to solve this problem ? 

May 25, 2015, 17:33 

#2  
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 50 
Quote:
__________________
Note: I don't use "Friend"feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request 

May 25, 2015, 21:17 

#3 
New Member
Join Date: May 2015
Location: Russia
Posts: 3
Rep Power: 9 
Thank you for the advice Bernhard!
At the moment I have added calls correctBoundaryConditions for required fields in "createFields" of the solver I am using. I have maked the comparison of solutions for two variants: 1. ...> funkySetFields > decomposePar>... and 2. ..> decomposePar > funkySetFields >... Solutions are equal now! 

May 26, 2015, 09:13 

#4  
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 50 
Quote:
__________________
Note: I don't use "Friend"feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request 

May 26, 2015, 09:48 

#5 
New Member
Join Date: May 2015
Location: Russia
Posts: 3
Rep Power: 9 

Tags 
funkysetfields, problem, processor boundary 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
UDF compiling problem  Wouter  Fluent UDF and Scheme Programming  6  June 6, 2012 04:43 
[swak4Foam] Problem with compiling funkySetFields  Victor  OpenFOAM Community Contributions  4  July 19, 2011 16:34 
Gambit  meshing over airfoil wrapping (?) problem  JFDC  FLUENT  1  July 11, 2011 05:59 
natural convection problem for a CHT problem  SeHee  CFX  2  June 10, 2007 06:29 
Adiabatic and Rotating wall (Convection problem)  ParodDav  CFX  5  April 29, 2007 19:13 