dynamic mesh for multi-region problem
 

Hi Foamers,

I have a question on dynamic mesh for multi-region.
Many examples in OpenFOAM demonstrated it's capability of dealing with
body moving behavior in the computation domain. But there is only one computational domain.

I'm thinking it is quite possible to model multi-computation domain with dynamic mesh technique for OpenFOAM, for example based on chtMultiRegionFoam solver. Such that the fields (velocity, temperature, etc.) in each domain can be obtained.

As a test, I simply put this line to the solid region,
const scalar scaleFactor(0.0);
volPointInterpolation pInterp(mesh);
pointField zeroPoints(mesh.points());
pointField newPoints
(
zeroPoints
+ scaleFactor*pInterp.interpolate(U)().internalField ()
);

mesh.polyMesh::movePoints(newPoints);
many error thrown out after one successful iteration.

Any suggestions are highly appreciated!

==================== ERROR ======================


#0 Foam::error::printStack(Foam::Ostream&) in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::fv::gaussGrad<double>::gradf(Foam::Geometric Field<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::word const&) in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 Foam::fv::gaussGrad<double>::calcGrad(Foam::Geomet ricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::fv::gradScheme<double>::grad(Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#6 Foam::fv::correctedSnGrad<double>::fullGradCorrect ion(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fv::correctedSnGrad<double>::correction(Foam ::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8 Foam::fv::limitedSnGrad<double>::correction(Foam:: GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#9 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#10 Foam::fv::laplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#11 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/chtMultiRegionDyMeshFoam"
#12 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/chtMultiRegionDyMeshFoam"
#13
in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/chtMultiRegionDyMeshFoam"
#14 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#15
in "/share/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/chtMultiRegionDyMeshFoam"
Floating point exception (core dumped)

what was the outcome?
 

Hi Yijin -

I face a problem similar to this - when attempting to use different dynamic mesh libraries for coupled solid and fluid regions.

Did you ever get close to an answer on this issue?

Best, Eric

HI Eric:

what is your problem now?
Let's see what I known can help you.

Dear Yijin,

Have you been able to simulate multi region mesh motion in OF?
I am facing the same problem.

Yes. I did. You can email to ymao.mu@gmail.com for detail.

Hi,

Could you please post solution to the question of dynamic mesh in multi region cases? I am currently also working on this problem.
Thank you

Hi all,

Would anyone be able to post a solution for this, would highly appreciate.

Thanks,
Dasith