[kk415]

Hi Foamers

I am trying to create a solver in which navier-stokes in one grid and temperature is solved in another grid. I developed this solver in a similar fashion to chtMultiRegion. I am stuck in the interpolation function of the solver. I am using the mapping utility of openFoam for this purpose.

I modified the MapConsistentMesh in mapfield.H for my purpose like this:

template<template<class> class CombineOp>
void MapConsistentMesh
(
const fvMesh& meshSource,
const fvMesh& meshTarget,
const meshToMesh0:rder& mapOrder,
IOobject& T,
IOobject& p,
IOobject& U,
IOobject& Ts,
IOobject& ps,
IOobject& Us
)
{
// Create the interpolation scheme
meshToMesh0 meshToMesh0Interp(meshSource, meshTarget);

Info<< nl
<< "Consistently creating and mapping fields for time "
<< meshSource.time().timeName() << nl << endl;

{
// Search for list of objects for this time
IOobjectList objects(meshSource, meshSource.time().timeName());

objects.add(U);
objects.add(p);
objects.add(T);
Info<< "object size" << objects.size()
<< endl;
// Map volFields
// ~~~~~~~~~~~~~
MapConsistentVolFields<scalar>
(
objects,
meshToMesh0Interp,
mapOrder,
CombineOp<scalar>()
);
MapConsistentVolFields<vector>
(
objects,
meshToMesh0Interp,
mapOrder,
CombineOp<vector>()
);
MapConsistentVolFields<sphericalTensor>
(
objects,
meshToMesh0Interp,
mapOrder,
CombineOp<sphericalTensor>()
);
MapConsistentVolFields<symmTensor>
(
objects,
meshToMesh0Interp,
mapOrder,
CombineOp<symmTensor>()
);
MapConsistentVolFields<tensor>
(
objects,
meshToMesh0Interp,
mapOrder,
CombineOp<tensor>()
);
Info<< " checking1 " << endl;
}

{
Info<< " checking " << endl;
// Search for list of target objects for this time
IOobjectList objects(meshTarget, meshTarget.time().timeName());
//objects.add(Us);
//objects.add(ps);
//objects.add(Ts);
// Mark surfaceFields as unmapped
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
UnMapped<surfaceScalarField>(objects);
UnMapped<surfaceVectorField>(objects);
UnMapped<surfaceSphericalTensorField>(objects);
UnMapped<surfaceSymmTensorField>(objects);
UnMapped<surfaceTensorField>(objects);

// Mark pointFields as unmapped
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~
UnMapped<pointScalarField>(objects);
UnMapped<pointVectorField>(objects);
UnMapped<pointSphericalTensorField>(objects);
UnMapped<pointSymmTensorField>(objects);
UnMapped<pointTensorField>(objects);
}
Info<< " checking " << endl;
// mapLagrangian(meshToMesh0Interp);
}

The output I am getting is like this :

Starting time loop

Time = 0.005


Solving for fluid region fluid

PISO: Operating solver in PISO mode

Reading transportProperties
fluid

Reading transportProperties
solid

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 8.90511e-06, No Iterations 19
smoothSolver: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.0492854, No Iterations 12
DICPCG: Solving for p, Initial residual = 0.590864, Final residual = 2.65225e-07, No Iterations 35

Consistently creating and mapping fields for time 0.005

object size3
interpolating T
checking2
interpolating p
checking2
interpolating U
checking2
checking1
*** glibc detected *** testFoam: double free or corruption (out): 0x00007ffc2f2b2fa0 ***
======= Backtrace: =========
/lib64/libc.so.6[0x3b0d275dee]
/lib64/libc.so.6[0x3b0d278c80]
/home/Rajesh/OpenFOAM/OpenFOAM-4.1/platforms/linux64GccDPInt64Opt/lib/libOpenFOAM.so(_ZN4Foam12HashPtrTableINS_8IOobject ENS_4wordENS_6string4hashEE5clearEv+0x44)[0x7ffadee37804]
/home/Rajesh/OpenFOAM/OpenFOAM-4.1/platforms/linux64GccDPInt64Opt/lib/libOpenFOAM.so(_ZN4Foam12IOobjectListD2Ev+0x9)[0x7ffadee35089]
testFoam[0x4f3564]
testFoam[0x44a2a9]
/lib64/libc.so.6(__libc_start_main+0xfd)[0x3b0d21ed1d]
testFoam[0x44d335]
======= Memory map: ========
00400000-0053c000 r-xp 00000000 08:11 60694501 /home/Rajesh/OpenFOAM/OpenFOAM-4.1/platforms/linux64GccDPInt64Opt/bin/testFoam
0073c000-00741000 rw-p 0013c000 08:11 60694501 /home/Rajesh/OpenFOAM/OpenFOAM-4.1/platforms/linux64GccDPInt64Opt/bin/testFoam
01e01000-02136000 rw-p 00000000 00:00 0 [heap]
3b0ce00000-3b0ce20000 r-xp 00000000 08:04 29884418 /lib64/ld-2.12.so
3b0d020000-3b0d021000 r--p 00020000 08:04 29884418 /lib64/ld-2.12.so
3b0d021000-3b0d022000 rw-p 00021000 08:04 29884418 /lib64/ld-2.12.so
3b0d022000-3b0d023000 rw-p 00000000 00:00 0
3b0d200000-3b0d38a000 r-xp 00000000 08:04 29884435 /lib64/libc-2.12.so
3b0d38a000-3b0d58a000 ---p 0018a000 08:04 29884435 /lib64/libc-2.12.so
3b0d58a000-3b0d58e000 r--p 0018a000 08:04 29884435 /lib64/libc-2.12.so
3b0d58e000-3b0d590000 rw-p 0018e000 08:04 29884435 /lib64/libc-2.12.so
3b0d590000-3b0d594000 rw-p 00000000 00:00 0
3b0d600000-3b0d617000 r-xp 00000000 08:04 29884442 /lib64/libpthread-2.12.so
3b0d617000-3b0d817000 ---p 00017000 08:04 29884442 /lib64/libpthread-2.12.so
3b0d817000-3b0d818000 r--p 00017000 08:04 29884442 /lib64/libpthread-2.12.so
3b0d818000-3b0d819000 rw-p 00018000 08:04 29884442 /lib64/libpthread-2.12.so
3b0d819000-3b0d81d000 rw-p 00000000 00:00 0
3b0da00000-3b0da02000 r-xp 00000000 08:04 29884533 /lib64/libdl-2.12.so
3b0da02000-3b0dc02000 ---p 00002000 08:04 29884533 /lib64/libdl-2.12.so
3b0dc02000-3b0dc03000 r--p 00002000 08:04 29884533 /lib64/libdl-2.12.so
3b0dc03000-3b0dc04000 rw-p 00003000 08:04 29884533 /lib64/libdl-2.12.so
3b0de00000-3b0de07000 r-xp 00000000 08:04 29884457 /lib64/librt-2.12.so
3b0de07000-3b0e006000 ---p 00007000 08:04 29884457 /lib64/librt-2.12.so


and so on ..... can anyone help me to know from what point this error is coming .... I mean I am getting the output as checking1 is printed. That means upto that point everything is correct. Then at what point the error is occuring????