CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

custom inlet problem with parallel runs

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 10, 2013, 12:46
Default custom inlet problem with parallel runs
  #1
New Member
 
Hansjoerg Seybold
Join Date: Mar 2009
Posts: 16
Rep Power: 7
hansjoerg is on a distinguished road
Hi,
I created custom velocity inlets for
a) linear shear flow
b) logarithmic velocity profile
to use them with pisoFoam.

For testing I run the following commands:
1) ~>decomposePar
2) ~>mpirun -np 4 pisoFoam -parallel -case .
3) ~>pisoFoam -case .

In the case of the linear inlet (a): 1-3 run without error
In the case of the logarithmic inlet (b): only (1) and (3) run

running pisoFoam in parallel however fails with the error message below.
The corresponding section in the velocity 0/U:
...
inlet
{
type logVelocity;
n (1 0 0);
y (0 0 1);
uStar 1.0;
z0 0.1;
kappa 0.41;
//value uniform (0 0 0); // for paraview
}

I dont know where to search for the error, as my configs/custom
libs run fine in the serial case.

When checking the processor0 - processor3 directories I found that the kappa entry was missing in the U dictionary.
After adding the entry manually I succeeded running the case
Hope somebody has an idea why decomposePar omits the entry

Thanks a lot
hj

-----------------------------------------------------------------------
...
[0] --> FOAM FATAL IO ERROR:
[0] keyword kappa is undefined in dictionary "/net/segovia/home/hansjoerg/Projects/Particles_in_Flow/CFDEM/hansjoerg-PUBLIC-2.1.x/pisoShearFlowCFDEM1/CFD/processor0/0/U::boundaryField::inlet"
[0]
[0] file: /net/segovia/home/hansjoerg/Projects/Particles_in_Flow/CFDEM/hansjoerg-PUBLIC-2.1.x/pisoShearFlowCFDEM1/CFD/processor0/0/U::boundaryField::inlet from line 36 to line 41.
[0]
[0] From function dictionary::lookupEntry(const word&, bool, bool) const
[0] in file db/dictionary/dictionary.C at line 400.
[0]
FOAM parallel run exiting
....
-----------------------------------------------------------------------

Last edited by hansjoerg; June 10, 2013 at 13:51.
hansjoerg is offline   Reply With Quote

Old   June 10, 2013, 13:58
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 6,981
Blog Entries: 32
Rep Power: 68
wyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the rough
Greetings Hansjoerg,

Not much information to work with... this is one of those "which I had a crystal ball" kind of moments...

Nonetheless, I do believe I know what the problem is, because I fixed one of these problems not too long ago. Check the source code for "libMyPolynomVelocity" from the following two posts:
Compare them side-by-side and you might see what the problem is in your own boundary condition

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   June 10, 2013, 14:40
Default
  #3
New Member
 
Hansjoerg Seybold
Join Date: Mar 2009
Posts: 16
Rep Power: 7
hansjoerg is on a distinguished road
Hi,
Thanks a lot for he reply. Problem solved.

I went through your suggested posts and I think that i found the
problem:

At the end of the library definition in the function

void logVelocityFvPatchVectorField::write(Ostream& os) const

the entry for "kappa":

os.writeKeyword("kappa")<< kappa_ << token::END_STATEMENT << nl;

was missing. This seems to be ok in serial jobs but crucial in parallel calculations.



Cheers

hj

Last edited by hansjoerg; June 10, 2013 at 18:28.
hansjoerg is offline   Reply With Quote

Old   November 15, 2013, 10:05
Default
  #4
Member
 
Join Date: May 2012
Posts: 47
Rep Power: 4
styleworker is on a distinguished road
How did you read the entry kappa outside the logVelocityFvPatchVectorField?
styleworker 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
Problem with assigned inlet velocity profile as a boundary condition Ozgur_ FLUENT 2 August 3, 2012 05:54
Flux Problem at Inlet piyushsharma0487 Fluent UDF and Scheme Programming 9 April 19, 2012 02:57
How to assign the inlet B.C using a bunch of data set for an unsteady problem? ali8500 CFX 3 March 28, 2012 18:41
Parallel runs with sonicDyMFoam crashes (works fine with sonicFoam) jnilsson OpenFOAM Running, Solving & CFD 0 March 9, 2012 06:45
parallel mode - small problem? co2 FLUENT 2 June 1, 2004 23:47


All times are GMT -4. The time now is 14:22.