[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.

[wyldckat]

Greetings Vishal,

You're still using OpenFOAM 1.5? Well, the following details should still work... but some of the details I went to look for them here: http://repo.or.cz/w/OpenFOAM-1.5.x.git/tree

If I understand you correctly: since you want to use multiple domains, then you will have to use a solver similar chtMultiRegionFoam. And more specifically for the mesh, you will have to study the tutorial "chtMultiRegionFoam/multiRegionHeater" (http://repo.or.cz/w/OpenFOAM-1.5.x.g...tiRegionHeater) for more information on how to set-up the multiple regions.

This is because you first have to split the mesh into the various regions, only then can you assign patches to the faces that are between those regions.

Best regards,
Bruno

[nandiganavishal]

Thanks Bruno.

I have installed the latest OF version 2.1.1 now !! I read the chtmultiregionFoam and started creating the geometry. The mesh looks consistent now. However, I have a few queries regarding incorporating BCs. The questions are attached in the doc file. Please let me know your suggestions.

Thanks

Regards
Vishal

[wyldckat]

Hi Vishal,

Sorry, but this is beyond my level of expertise. All I can figure out is that you will need to take one of the "chtMultiRegion*Foam" solvers and modify the source code for solving your equations instead.

There is also a "chtIcoMultiRegionFoam" made by someone in the community, which you can use as a coding reference: http://www.cfd-online.com/Forums/ope...egionfoam.html

Best regards,
Bruno

[nandiganavishal]

Thanks Bruno.

1. Can you let me know what does the boundary condition "calculated" mean in chtMultiregionFoam!!

2. Further what does the boundary condition compressible::turbulentTemperatureCoupledBaffleMix ed; physically mean ? Can you give me the mathematical form of this equation !!

Thanks

Regards
Vishal

[wyldckat]

Hi Vishal,

1. I don't know much about the "calculated" patch. It rarely works for me and the user guide says:
   

   Originally Posted by http://www.openfoam.org/docs/user/boundaries.php

   calculated - Boundary field derived from other fields

2. All I know is where you can find more information about it. In the dedicated header file you will find some information on the "Description" part: https://github.com/OpenFOAM/OpenFOAM...hScalarField.H

Best regards,
Bruno

[nandiganavishal]

Thanks Bruno !!