[Saleh Abuhanieh]

Hi Foamers,

I am trying to develop a solver that can solve the following equations for different regions. I had to segregate the domain for different regions since a special BC shall be maintained between them.

Equations from the standard solver are:

solve
(
fvm::laplacian(phi) + rho/epsilon
);


rhoFlux = -k*mesh.magSf()*fvc::snGrad(phi);

solve
(
fvm::ddt(rho) + fvm::div(rhoFlux, rho)
);

Noting that before going to the solving step, all the associated fields and constants in the equations are replaced by the regions related ones. I used ptrLists to define the same names for fields and constants for each region.

The solver compiled without error, but when I run the case, it crashes once starts solving:

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::UPtrList<Foam::fvsPatchField<double> >:perator[](int) const in "/home/sfa/OpenFOAM/sfa-5.0/platforms/linux64GccDPInt32Opt/bin/electrostaticMultiRegionFoam"
#4 ? in "/home/sfa/OpenFOAM/sfa-5.0/platforms/linux64GccDPInt32Opt/bin/electrostaticMultiRegionFoam"
#5 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:?
#6 Foam::fv::gaussLaplacianScheme<double, Foam::SymmTensor<double> >::fvmLaplacianUncorrected(Foam::GeometricField<do uble, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#7 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#8 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double>(Foam::GeometricField< double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/sfa/OpenFOAM/sfa-5.0/platforms/linux64GccDPInt32Opt/bin/electrostaticMultiRegionFoam"
#9 ? in "/home/sfa/OpenFOAM/sfa-5.0/platforms/linux64GccDPInt32Opt/bin/electrostaticMultiRegionFoam"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? in "/home/sfa/OpenFOAM/sfa-5.0/platforms/linux64GccDPInt32Opt/bin/electrostaticMultiRegionFoam"
Segmentation fault (core dumped)

I thought maybe I need to define (explicitly) a fvMatrix for each region using ptrList also but i didn't know how since I've two equations, I tried one matrix for each equation just to see what will happen:

phiMatrix
(
fvm::laplacian(phi) + rho/epsilon
);

phiMatrix.solve();


rhoFlux = -k*mesh.magSf()*fvc::snGrad(phi);


rhoMatrix
(
fvm::ddt(rho) + fvm::div(rhoFlux, rho)
);

rhoMatrix.solve();

I received the the following compilation error:

error: no match for call to ‘(Foam::fvScalarMatrix {aka Foam::fvMatrix<double>}) (Foam::tmp<Foam::fvMatrix<double> >)’
);


Any help will be appreciated

Saleh

[Saleh Abuhanieh]

any clue pls.

[Saleh Abuhanieh]

please notify me if the provided information is not sufficient