|
[Sponsors] |
BC which interpolate from an existing solution |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 5, 2013, 05:59 |
BC which interpolate from an existing solution
|
#1 |
New Member
Marcel Vonlanthen
Join Date: Nov 2012
Location: Zurich, Switzerland
Posts: 28
Rep Power: 14 |
Dear foamer,
I try to create a boundary condition which interpolate the volumic velocity field of an existing case (sourceCase) to the inlet patch of a running case (runCase). To do this, I use the interpolate static method include in surfaceInterpolation.H. Here is a piece of code included in the "updateCoeff" member function of my boundary: Code:
fileName sourceRootDir = "path/to/source"; fileName sourceCaseDir = "sourceCase"; fileName sourceTime = "1000"; Time runTimeSource ( Time::controlDictName, sourceRootDir, sourceCaseDir ); fvMesh meshSource ( IOobject ( sourceRegion, runTimeSource.timeName(), runTimeSource ) ); volVectorField Usource ( IOobject ( "U", sourceTime, meshSource, IOobject::MUST_READ, IOobject::NO_WRITE ), meshSource ); // do the interpolation operator==(fvc::interpolate(Usource)); fixedValueFvPatchField<vector>::updateCoeffs(); Code:
--> FOAM FATAL IO ERROR: keyword interpolate(U) is undefined in dictionary "/path/to/source/sourceCase/system/fvSchemes.interpolationSchemes" file: /path/to/source/sourceCase/system/fvSchemes.interpolationSchemes From function dictionary::lookupEntry(const word&, bool, bool) const in file db/dictionary/dictionary.C at line 402. FOAM exiting Code:
interpolationSchemes { interpolate(U) linear; } |
|
July 10, 2013, 07:06 |
|
#2 |
Senior Member
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 168
Rep Power: 16 |
Dear Marcel,
I have the same problem, did you find a solution? Kind Regards Anne |
|
July 11, 2013, 06:14 |
|
#3 | |
New Member
Marcel Vonlanthen
Join Date: Nov 2012
Location: Zurich, Switzerland
Posts: 28
Rep Power: 14 |
Quote:
I fixed the error "interpolate(U)" by just adding one line in the IOobject of meshSource Code:
fvMesh meshSource ( IOobject ( sourceRegion, runTimeSource.timeName(), runTimeSource, IOobject::MUST_READ //this makes all the difference... ) ); In the meantime, I used a different approach, but let me know if you have a working boundary! Cheers, Marcel |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Asymmetrical solution using RANS | Marli | FLUENT | 2 | January 9, 2015 10:26 |
grid dependancy | gueynard a. | Main CFD Forum | 19 | June 27, 2014 22:22 |
Why 3D solid-pore geometry showing diverged solution? | Sargam05 | OpenFOAM | 0 | December 3, 2012 16:45 |
Analytic solution for 2D steady Euler equations | jojo81 | Main CFD Forum | 0 | October 15, 2012 13:05 |
Wall functions | Abhijit Tilak | Main CFD Forum | 6 | February 5, 1999 02:16 |