[irishdave]

Hello,

I have been looking at the timeVaryingMappedFixedValue fvPatch. I am currently making a program that will be used to create inlet data. This patch looks to be the perfect solution in order to use the created data within a case. My question is where can I find more documentation on the patch? I want to be able to create the constant/boundaryData/inlet files.

What is contained in the points file and is there an easy way to create one?

Why for the tutorial case simpleFoam/pitzDailyExptInlet/ which uses the BC do the files contain 70 values when the fvPatch faceCells() function returns 30 cells for the case?

I am currently trying to implement a inflow data generator for OpenFOAM which upon completion will be made available open-source to the OpenFOAM community.

Thanks
David

[mattijs]

The points are the locations of the data. The data at those points is in the time directory.

The patch interpolates the point data onto the face centre.

The point data comes from e.g. measurements.

[irishdave]

Thank you Mattijs, I played around with it a bit and started to realize what what was going on.

[jerome]

Hello,

I would like to use timeVaryingMappedFixedValue for my inlet velocities and I want to set different velocity values for each time step of my simulation (in this case, deltaT = 0.0002 s).

I created the directories in constant/boundaryData/inlet/:

constant/boundaryData/inlet/points
constant/boundaryData/inlet/0/U
constant/boundaryData/inlet/0.0002/U
constant/boundaryData/inlet/0.0004/U
constant/boundaryData/inlet/0.0006/U
(…)

However, when I run the case, I have a problem at the first time step:

Starting time loop

Time = 0.0002

Courant Number mean: 0 max: 0.241854


--> FOAM FATAL IO ERROR : cannot open file

file: …/run/ExpeCylinder1uP7/constant/boundaryData/inlet/0.0004 at line0.

From function regIOobject::readStream(const word&)
in file db/regIOobject/regIOobjectRead.C at line 66.

FOAM exiting


It runs nevertheless if I only provide the directory "0" (constant inlet velocity values).
Could someone tell me the correct way to use this patch?

Thank you in advance,
Jerome

[mattijs]

Switch on the debug flag for timeVaryingMappedFixedValue. It should print

timeVaryingMappedFixedValueFvPatchField : In directory XXXX found times YYYYY. See if it includes your 0.0004 directory.

Does the $FOAM_TUTORIALS/simpleFoam/pitzDailyExptInlet case work if you add another time directory?

[jerome]

Hello,

Thank you for your answer,

Yes, I have the same problem with pitzDailyExptInlet: the simulation stopped at the first time step and said:

--> FOAM FATAL IO ERROR : cannot open file

file: .../pitzDailyExptInlet/constant/boundaryData/inlet/2 at line 0.

From function regIOobject::readStream(const word&)
in file db/regIOobject/regIOobjectRead.C at line 66.

FOAM exiting


I switched on the debug flag for timeVaryingMappedFixedValue. I obtain:

Reading field p

Reading field U

timeVaryingMappedFixedValue : construct from dictionary

timeVaryingMappedFixedValueFvPatchField : Read 100 sample points from ".../ExpeCylinder1uP7/constant/boundaryData/inlet/points"
timeVaryingMappedFixedValueFvPatchField : Used points (0.00138321 0.00138321 0) (0.00154865 0.00111043 0) (0.00167642 0.000812023 0) to define coordinate system with normal (0 0 -1)

readSamplePoints : Dumping triangulated surface to triangulation.stl
readSamplePoints : Dumping face centres to ".../ExpeCylinder1uP7/localFaceCentres.obj"

timeVaryingMappedFixedValueFvPatchField : In directory ".../ExpeCylinder1uP7/constant/boundaryData/inlet" found times

9
(
0
0.0002
0.0004
0.0006
0.0008
0.001
0.0012
0.0014
0.0016
)

findTime : Found time 0 inbetween index:0 time:0 and index:1 time:0.0002

checkTable : Reading startValues from "boundaryData/inlet/0"
checkTable : Reading endValues from "boundaryData/inlet/0.0002"

updateCoeffs : Sampled, interpolated values between start time:0 and end time:0.0002 with weight:0
updateCoeffs : set fixedValue to min0 0 0.0134146) max0 0 0.302318)

timeVaryingMappedFixedValue : copy construct, resetting internal field
Reading/calculating face flux field phi


Starting time loop

Time = 0.0002

Courant Number mean: 0 max: 0.241854

findTime : Found time 0.0002 inbetween index:1 time:0.0002 and index:2 time:0.0004
checkTable : Setting startValues to (already read) "boundaryData/inlet/0.0002"
checkTable : Reading endValues from "boundaryData/inlet/0.0004"


--> FOAM FATAL IO ERROR : cannot open file

file: .../ExpeCylinder1uP7/constant/boundaryData/inlet/0.0004 at line0.

From function regIOobject::readStream(const word&)
in file db/regIOobject/regIOobjectRead.C at line 66.

FOAM exiting




Do you know what the problem could be?

Thank you in advance,

Jerome

[mattijs]

Sounds like a bug. Can you report it in the bug-reports section and we'll have a look at it.

[johndeas]

I have a similar problem using timeVaryingMappedFixedValue:

Time = 0.2

Courant Number mean: 0.0215601 max: 0.1955


--> FOAM FATAL IO ERROR : cannot open file

file: /home/vinz/OpenFOAM/vinz-1.4.1/run/05-TVMFV/constant/boundaryData/entree/0.4 at line 0.

From function regIOobject::readStream(const word&)
in file db/regIOobject/regIOobjectRead.C at line 66.

FOAM exiting

It does not appear when I have two time directories in boundaryData, but when a third is added (even just by copying/renaming another one), the previous error message appears.

[johndeas]

And it appears when the second directory has to be read.

For example, if I have the folders 0 0.1 and 0.2 in boundaryData, the problem will appear at Time = 0.1

[johndeas]

Here is the report using debug flag, with directories 0 0.1 and 0.2:

Time = 0.1

Courant Number mean: 0.0215667 max: 0.191
findTime : Found time 0.1 inbetween index:1 time:0.1 and index:2 time:0.2
checkTable : Setting startValues to (already read) "boundaryData/entree/0.1"
checkTable : Reading endValues from "boundaryData/entree/0.2"


--> FOAM FATAL IO ERROR : cannot open file

file: /home/vinz/OpenFOAM/vinz-1.4.1/run/05-TVMFV/constant/boundaryData/entree/0.2 at line 0.

From function regIOobject::readStream(const word&)
in file db/regIOobject/regIOobjectRead.C at line 66.

FOAM exiting

[mattijs]

I don't see this problem. Running the pitzDailyExptInlet tutorial case now with 4 time files. (0,0.2,0.4 and 0.9). Do you have a 0.2 file? Readable?

[johndeas]

Yes, even if it is a direct copy of the 0 it doesn't work.

[johndeas]

Here is how you can reproduce my errors, considering three different cases:

take the case pitzDailyExptinlet

change controlDict/endTime to 20 (just to make it quick)

several cases now:

case 1: if I copy boundaryData/0 to boundaryData/25, it gives me:
------

Reading field U

timeVaryingMappedFixedValue : construct from dictionary
timeVaryingMappedFixedValueFvPatchField : Read 70 sample points from "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/constant/boundaryData/ inlet/points"
timeVaryingMappedFixedValueFvPatchField : Used points (-0.0206 0.000125 0) (-0.0206 0.000375 0) (-0.0206 0.000125 0.173) to define coordinate system with normal (1 0 0)
readSamplePoints : Dumping triangulated surface to triangulation.stl
readSamplePoints : Dumping face centres to "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/localFaceCentres.obj"
timeVaryingMappedFixedValueFvPatchField : In directory "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/constant/boundaryData/ inlet" found times
2
(
25
0
)



--> FOAM FATAL ERROR : Cannot find starting sampling values for current time 0
Have sampling values for times
2
(
25
0
)

In directory "constant/boundaryData/inlet"
on patch inlet of field U in file "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/0/U"

From function findTime
in file fields/fvPatchFields/derived/timeVaryingMappedFixedValue/timeVaryingMappedFixedV alueFvPatchField.C at line 469.

FOAM exiting

case 2: if I copy boundaryData/0 to boundaryData/21, -or- boundaryData/26 it works fine, I do not understand why, since 21, 25 and 26 are above the upper limit of the calculation, it might not be a problem of interpolation, why the bug with 25 ?
-----

case 3: if I copy boundaryData/0 to boundaryData/3 -and- boundaryData/26, the bug appear at Time=2 with the following message:
-----

Time = 2

findTime : Found time 2 inbetween index:1 time:2 and index:2 time:26
checkTable : Setting startValues to (already read) "boundaryData/inlet/2"
checkTable : Reading endValues from "boundaryData/inlet/26"


--> FOAM FATAL IO ERROR : cannot open file

file: /home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/constant/boundaryData/i nlet/26 at line 0.

From function regIOobject::readStream(const word&)
in file db/regIOobject/regIOobjectRead.C at line 66.

FOAM exiting

I hope those reproducible errors might help to point at a misnapulation of my own/bug in the patch.

[mattijs]

The problem seems to be the sorting of the times. I get:

timeVaryingMappedFixedValueFvPatchField : In directory "../pitzDailyExptInlet/constant/boundaryData/inlet" found times
2
(
0
25
)

which is in the correct, incremental, order. (findTimes should return times in sorted order). Don't directly see why you should get them differently.

[johndeas]

I tried again today, copying boundaryData/0 to 5 and 30, and get the following message:

(in this case, findTime provide them in the right order...)

Exec : simpleFoam . 07-pitzDailyExptInlet
Date : May 14 2008
Time : 12:50:05
Host : pcaea1
PID : 9819
Root : /home/vinz/OpenFOAM/vinz-1.4.1/run
Case : 07-pitzDailyExptInlet
Nprocs : 1
Create time

Create mesh for time = 0

Reading field p

Reading field U

timeVaryingMappedFixedValue : construct from dictionary
timeVaryingMappedFixedValueFvPatchField : Read 70 sample points from "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyE
xptInlet/constant/boundaryData/inlet/points"
timeVaryingMappedFixedValueFvPatchField : Used points (-0.0206 0.000125 0) (-0.0206 0.000375 0) (-0.0206 0.000125 0.17
3) to define coordinate system with normal (1 0 0)
readSamplePoints : Dumping triangulated surface to triangulation.stl
readSamplePoints : Dumping face centres to "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/localFaceCentres.
obj"
timeVaryingMappedFixedValueFvPatchField : In directory "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/const
ant/boundaryData/inlet" found times
3
(
0
5
30
)

findTime : Found time 0 inbetween index:0 time:0 and index:1 time:5
checkTable : Reading startValues from "boundaryData/inlet/0"
checkTable : Reading endValues from "boundaryData/inlet/5"
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0
updateCoeffs : set fixedValue to min1.51773 0 0) max14.3576 0 0)
timeVaryingMappedFixedValue : copy construct, resetting internal field
Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model kEpsilon
timeVaryingMappedFixedValue : construct from dictionary
timeVaryingMappedFixedValueFvPatchField : Read 70 sample points from "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyE
xptInlet/constant/boundaryData/inlet/points"
timeVaryingMappedFixedValueFvPatchField : Used points (-0.0206 0.000125 0) (-0.0206 0.000375 0) (-0.0206 0.000125 0.17
3) to define coordinate system with normal (1 0 0)
readSamplePoints : Dumping triangulated surface to triangulation.stl
readSamplePoints : Dumping face centres to "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/localFaceCentres.
obj"
timeVaryingMappedFixedValueFvPatchField : In directory "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/const
ant/boundaryData/inlet" found times
3
(
0
5
30
)

findTime : Found time 0 inbetween index:0 time:0 and index:1 time:5
checkTable : Reading startValues from "boundaryData/inlet/0"
checkTable : Reading endValues from "boundaryData/inlet/5"
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0
updateCoeffs : set fixedValue to min:0.0828453 max:3.12421
timeVaryingMappedFixedValue : copy construct, resetting internal field
timeVaryingMappedFixedValue : construct from dictionary
timeVaryingMappedFixedValueFvPatchField : Read 70 sample points from "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyE
xptInlet/constant/boundaryData/inlet/points"
timeVaryingMappedFixedValueFvPatchField : Used points (-0.0206 0.000125 0) (-0.0206 0.000375 0) (-0.0206 0.000125 0.17
3) to define coordinate system with normal (1 0 0)
readSamplePoints : Dumping triangulated surface to triangulation.stl
readSamplePoints : Dumping face centres to "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/localFaceCentres.
obj"
timeVaryingMappedFixedValueFvPatchField : In directory "/home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/const
ant/boundaryData/inlet" found times
3
(
0
5
30
)

findTime : Found time 0 inbetween index:0 time:0 and index:1 time:5
checkTable : Reading startValues from "boundaryData/inlet/0"
checkTable : Reading endValues from "boundaryData/inlet/5"
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0
updateCoeffs : set fixedValue to min:2.27183 max:9813.84
timeVaryingMappedFixedValue : copy construct, resetting internal field

Starting time loop

Time = 1

findTime : Found time 1 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.2
updateCoeffs : set fixedValue to min1.51773 0 0) max14.3576 0 0)
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.0880712, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0415937, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00901764, No Iterations 158
time step continuity errors : sum local = 1.23171, global = -0.000635655, cumulative = -0.000635655
findTime : Found time 1 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.2
updateCoeffs : set fixedValue to min1.51773 0 0) max14.3576 0 0)
findTime : Found time 1 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.2
updateCoeffs : set fixedValue to min:2.27183 max:9813.84
DILUPBiCG: Solving for epsilon, Initial residual = 0.487709, Final residual = 0.0374107, No Iterations 1
findTime : Found time 1 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.2
updateCoeffs : set fixedValue to min:0.0828453 max:3.12421
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.0202643, No Iterations 1
ExecutionTime = 0.54 s ClockTime = 0 s

Time = 2

findTime : Found time 2 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.4
updateCoeffs : set fixedValue to min1.51773 0 0) max14.3576 0 0)
DILUPBiCG: Solving for Ux, Initial residual = 0.554546, Final residual = 0.0041229, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.358215, Final residual = 0.0337988, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.257793, Final residual = 0.00244823, No Iterations 167
time step continuity errors : sum local = 1.26227, global = -0.00138772, cumulative = -0.00202337
findTime : Found time 2 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.4
updateCoeffs : set fixedValue to min1.51773 0 0) max14.3576 0 0)
findTime : Found time 2 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.4
updateCoeffs : set fixedValue to min:2.27183 max:9813.84
DILUPBiCG: Solving for epsilon, Initial residual = 0.798777, Final residual = 0.00811118, No Iterations 2
findTime : Found time 2 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.4
updateCoeffs : set fixedValue to min:0.0828453 max:3.12421
DILUPBiCG: Solving for k, Initial residual = 0.759915, Final residual = 0.00607932, No Iterations 2
ExecutionTime = 0.76 s ClockTime = 1 s

Time = 3

findTime : Found time 3 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.6
updateCoeffs : set fixedValue to min1.51773 0 0) max14.3576 0 0)
DILUPBiCG: Solving for Ux, Initial residual = 0.442564, Final residual = 0.0389076, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.236774, Final residual = 0.0204506, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.133109, Final residual = 0.00123854, No Iterations 165
time step continuity errors : sum local = 0.808968, global = 0.00450947, cumulative = 0.00248609
findTime : Found time 3 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.6
updateCoeffs : set fixedValue to min1.51773 0 0) max14.3576 0 0)
findTime : Found time 3 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.6
updateCoeffs : set fixedValue to min:2.27183 max:9813.84
DILUPBiCG: Solving for epsilon, Initial residual = 0.335955, Final residual = 0.00300697, No Iterations 2
findTime : Found time 3 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.6
updateCoeffs : set fixedValue to min:0.0828453 max:3.12421
DILUPBiCG: Solving for k, Initial residual = 0.228812, Final residual = 0.00174147, No Iterations 2
ExecutionTime = 0.97 s ClockTime = 1 s

Time = 4

findTime : Found time 4 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.8
updateCoeffs : set fixedValue to min1.51773 0 0) max14.3576 0 0)
DILUPBiCG: Solving for Ux, Initial residual = 0.310706, Final residual = 0.00289445, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.14289, Final residual = 0.00143542, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.208651, Final residual = 0.00205794, No Iterations 171
time step continuity errors : sum local = 0.521486, global = -0.000293979, cumulative = 0.00219212
findTime : Found time 4 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.8
updateCoeffs : set fixedValue to min1.51773 0 0) max14.3576 0 0)
findTime : Found time 4 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.8
updateCoeffs : set fixedValue to min:2.27183 max:9813.84
DILUPBiCG: Solving for epsilon, Initial residual = 0.139524, Final residual = 0.0130092, No Iterations 1
findTime : Found time 4 inbetween index:0 time:0 and index:1 time:5
updateCoeffs : Sampled, interpolated values between start time:0 and end time:5 with weight:0.8
updateCoeffs : set fixedValue to min:0.0828453 max:3.12421
DILUPBiCG: Solving for k, Initial residual = 0.119173, Final residual = 0.0105309, No Iterations 1
ExecutionTime = 1.23 s ClockTime = 1 s

Time = 5

findTime : Found time 5 inbetween index:1 time:5 and index:2 time:30
checkTable : Setting startValues to (already read) "boundaryData/inlet/5"
checkTable : Reading endValues from "boundaryData/inlet/30"


--> FOAM FATAL IO ERROR : cannot open file

file: /home/vinz/OpenFOAM/vinz-1.4.1/run/07-pitzDailyExptInlet/constant/boundaryData/i nlet/30 at line 0.

From function regIOobject::readStream(const word&)
in file db/regIOobject/regIOobjectRead.C at line 66.

[mattijs]

Hi John,

there is a space in the 'i nlet' in the message above. Is that just from editing?

- any difference if you up the writePrecision to 15?
- anything special about your architecture/disk (nfs?)
- is this release 1.4.1?

All works fine for me. At time 5:

Time = 5

findTime : Found time 5 inbetween index:1 time:5 and index:2 time:30
checkTable : Setting startValues to (already read) "boundaryData/inlet/5"
checkTable : Reading endValues from "boundaryData/inlet/30"
updateCoeffs : Sampled, interpolated values between start time:5 and end time:30 with weight:0
updateCoeffs : set fixedValue to min1.51773 0 0) max14.3576 0 0)

[ivanwhlau]

hi,
I am trying to learn to use timeVaryingUniformFixedValue
However, I got the following error.
In fact, I am trying the sample from the tutorial. Also, I try to use the directMapped value from 30 points instead. Is it possible?
Can anyone give me some advise?
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field p

Reading field U



Cannot find points that make valid normal.
Need at least three sample points which are not in a line.
on patch inlet of field U in file "/home/ilau/OpenFOAM/ilau-1.5/run/test/simpleFoam/pitzDaily-timeVaryingMappedFix edValue/0/U"

From function timeVaryingMappedFixedValueFvPatchField<type>::rea dSamplePoints()
in file fields/fvPatchFields/derived/timeVaryingMappedFixedValue/timeVaryingMappedFixedV alueFvPatchField.C at line 330.

FOAM exiting

Regards,
Ivan

[miku]

Hi,

I have used timeVaryingMappedFixedValue for creating boundaryData. Now I also want to have fixed values for my internal field. Does anybody know if I can use timeVaryingMappedFixedValue for this problem? Perhaps by creating an folder for internalData instead of boundaryData? I want to get fixed values for pressure and velocity for all timesteps.

Thanks in advance.

[sega]

Hello again.

Is there any hint on how to use the timeVaryingMappedFixedValue boundary condition?

Unfortunately there is not much commend in the code and I'm not sure how to use the bc.

[johndeas]

Hi,

I am having some trouble setting up the timeVarryingMappedFixedValue boundary condition for a direct numerical simulation case. I thought it might be appropriate to advertise to my fellow tVMFV users, for my post located at http://www.cfd-online.com/Forums/ope...tml#post264064 .