[nandiganavishal]

Dear Foamers,

I would like to solve the following equation on a multiple domains as shown in the attachment !!

Equation: div.(A1 grad(phi)) = B

The numbers mentioned in the figure indicate the boundaries where the necessary boundary conditions have to be specified.

1. While incorporating the geometry in the blockMesh, I found that the red and the yellow regions were getting defined as one domain and I was not able to define the boundaries 7,8,9,15,16,17,18,19. I got the following error when the blockMesh is executed.

Trying to specify a boundary face 4(34 35 91 90) on the face on cell 13 which is either an internal face or already belongs to some other patch. This is face 0 of patch 8 named alumina_nanopore_upperwall.#0 Foam::error:rintStack(Foam::O stream&) in "/home/cmn/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/cmn/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt /libOpenFOAM.so"
#2 Foam:olyMesh:olyMesh(Foam::IOobject const&, Foam::Field<Foam::Vector<dou ble> > const&, Foam::List<Foam::cellShape> const&, Foam::List<Foam::List<Foam::f ace> > const&, Foam::List<Foam::word> const&, Foam::List<Foam::word> const&, Foa m::word const&, Foam::word const&, Foam::List<Foam::word> const&, bool) in "/hom e/cmn/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#3
in "/home/cmn/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/blockMesh"
#4
in "/home/cmn/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/blockMesh"
#5
in "/home/cmn/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/blockMesh"
#6 __libc_start_main in "/lib/libc.so.6"
#7 __gxx_personality_v0 in "/home/cmn/OpenFOAM/OpenFOAM-1.5/applications/bin/li nux64GccDPOpt/blockMesh"


Please let me know, how to modify the blockMesh so that both all the boundaries are consistently defined. Attached please find the blockMesh file.