||May 9, 2010 23:50
Using mapFields with interDyMFoam
Has anybody else run into problems using mapFields and interDyMFoam?
I am running a model up to about 1.5 seconds, then stopping it in order to make a change to the boundary conditions. I use mapFields to map the last timestep of the initial run to the first timestep of the second stage run, and it appears to interpolate the fields correctly.
However, once I start the second stage of the run, it processes a few iterations and then fails with a floating point error.
#0 Foam::error: printStack(Foam::Ostream&) in "/usr/OpenFoam/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/usr/OpenFoam/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::Time::adjustDeltaT() in "/usr/OpenFoam/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 main in "/usr/OpenFoam/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/interDyMFoam"
#5 __libc_start_main in "/lib/libc.so.6"
#6 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
I have tried a number of different ways to do this, however all seem to meet with failure. I am assuming that it isn't carrying the data relating to the mesh refinement from the first stage run into the second stage, but is there anything that I can do to overcome this?