|
[Sponsors] |
BC which interpolate from an existing solution |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 5, 2013, 04:59 |
BC which interpolate from an existing solution
|
#1 |
New Member
Marcel Vonlanthen
Join Date: Nov 2012
Location: Zurich, Switzerland
Posts: 28
Rep Power: 13 |
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, 06:06 |
|
#2 |
Senior Member
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 168
Rep Power: 15 |
Dear Marcel,
I have the same problem, did you find a solution? Kind Regards Anne |
|
July 11, 2013, 05:14 |
|
#3 | |
New Member
Marcel Vonlanthen
Join Date: Nov 2012
Location: Zurich, Switzerland
Posts: 28
Rep Power: 13 |
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 09:26 |
grid dependancy | gueynard a. | Main CFD Forum | 19 | June 27, 2014 21:22 |
Why 3D solid-pore geometry showing diverged solution? | Sargam05 | OpenFOAM | 0 | December 3, 2012 15:45 |
Analytic solution for 2D steady Euler equations | jojo81 | Main CFD Forum | 0 | October 15, 2012 12:05 |
Wall functions | Abhijit Tilak | Main CFD Forum | 6 | February 5, 1999 01:16 |