[Arslan Arshad]

I am new to OpenFOAM so i dont know how to solve the laplacianfoam for the case of two materials present in a slab?
in that case how to mention two region properties and where?
Please can anyone help me.


Thanks in Advance

[mostafa R]

Hi
u can define new volscalarField of properties for your domain and use it in your code.

[Arslan Arshad]

Hi
Thanks for the reply, but tell me how i am new to OpenFoam.

[mostafa R]

define a new volscalarfield in file "arslan.H" for example and type :
//---------------------------------------------------------
#include "createMesh.H"

volScalarField dT
(
IOobject
(
"dT",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
//--------------------------------------------------------------

("dT" is such diffusion number)

if u want solve thermal or electrical field in your domain with different material this file has read with your code.

so u must create "dT" file such as "T" file in your "0" folder (0 - constant - system)

after this, for initializing the dT constant properties like u want append subsequent code end of write "volScalarField dT ......" :::::

//--------------------------------------------------------------------

dT.intenalField() [ mesh.cellZones()[0] ] = dT_0;
dT.intenalField() [ mesh.cellZones()[1] ] = dT_1;
.
.
.
//---------------------------------------------------------------------


because of error in "check dimension" set dimension dimensionless :
@ T , dT file
[0 0 0 0 0 0 0 0]

in body of main source code :

//--------------------------------------------------
#include "arslan.H"
slove(..
.
.
.
.
)
//--------------------------------------------------

[mostafa R]

of course in T file dimension set :

[0 2 0 0 0 0 0 0]

for eliminate Laplace operator dimension

[mostafa R]

plz che ck this solver:

chtMultiRegionFoam

perhaps good for your attend

[Arslan Arshad]

"dT.intenalField() [ mesh.cellZones()[0] ] = dT_0;
dT.intenalField() [ mesh.cellZones()[1] ] = dT_1;"


I want to make Slab (2-D) geometry in gambit and want to divide my length in 2 equal parts lets say 0.5m each.
Then how to mention cell zone in my case

[mostafa R]

in gambit create 2 rectangle volume (for example) that face between them is "interior"
remember geometry must be 3D;

set name of left volume be "left"
set name of right volume be "right"
(next tab u set bc in gambit : "fluid")

and type of bc must be wall with different name(this is for refuse unscheduled ERROR)
after using " fluent3DTo.... name.msh " command in boundary file change the type of bc from wall to everything u want...

when this is done u have :
cellzones - faceZones file and u can see inner cellZones file like this:

'''''''''''''''''''''''''''''''''''''''''
4
(
liq
{
type cellZone;
cellLabels List<label>
10620
(
0
1
2
3
4
.
.
.
''''''''''''''''''''''''''''''''''''''''''''''''

[Arslan Arshad]

Thank you for your help

I have tried with the geometry you told and coding but something is wrong.
I dont know what is wrong. Please tell me your email id. I will send the geometry files. Check it.

I shall be very thankful to you.

[Arslan Arshad]

Link to geometry file is given below:

https://www.mediafire.com/?52cc7cd7x2qwge8

[mostafa R]

hi
i check this and reply to u

[Arslan Arshad]

Hi
Thank You

[mostafa R]

hi again

I look inside file and check it in gambit:

you first create sides (out cubic) BC (interior - left - right - and sides that not applied)
and export *.msh file

create root folder for your running on linux (without any " "/space/ character like : "new folder")
it has an ERROR

inside this file create and put (constant - system) sub root folder
in constant folder : create empty folder polyMesh

finally run this :
fluent3DMeshToFoam *.msh

* is file name

after this ,go to polyMesh folder and look inside CellZone - faceZone files

[Arslan Arshad]

I tried the procedure u told me.
i made "arslan.H" file and include in my code "#include arslan.H"

but when i gave command in terminal "wmake" for compilation of dependency file i got the following error.
--------------------------------------------------------------------------------------------
In file included from laplacianFoam.C:71:0:
arslan.H: In function ‘int main(int, char**)’:
arslan.H:19:4: error: no match for call to ‘(Foam::volScalarField {aka Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}) ()’
arslan.H:20:4: error: no match for call to ‘(Foam::volScalarField {aka Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}) ()’
arslan.H:39:7: error: no match for call to ‘(Foam::volScalarField {aka Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}) ()’
arslan.H:40:7: error: no match for call to ‘(Foam::volScalarField {aka Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}) ()’
make: *** [Make/linux64GccDPOpt/laplacianFoam.o] Error 1
-------------------------------------------------------------------------------------------

what should i do?

[Arslan Arshad]

Also i am attaching the file "arslan.H" please check it too.
its urgent i have to submit this work in a week.

[Arslan Arshad]

Thanks in advance.

[mostafa R]

i see.
//---------------------------------------------------------//
#include "createMesh.H"

Info<< "Reading field DT\n" << endl;

volScalarField DT
(
IOobject
(
"DT",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

// edit :

double dt[2] = {0.0087, 0.0084};

forAll(mesh.cellZones(),i)
{
forAll(mesh.cellZones()[i],j)
{
int index = mesh.cellZoens()[i][j];
DT.internalField()[index] = dt[i]; // (i don't remember - maybe DT.internalFields() )
}
}

////////////////////////////////

do this for k

and then if the BC not set u can set them
by this variables:

DT.boundaryField()[][]

mesh.faceZones() // like zones this is for interface BC

enjoy it

[Arslan Arshad]

i am still facing errors i will send you my solver with geometry.......i hope u will check it.

[mostafa R]

send your errors

i do this and works good

[Arslan Arshad]

I am attaching the solver
please u have the geometry link
and there is nothing in "polyMesh" folder it is empty.
(please use the command fluent3DMeshToFoam HSlab.msh)

and when you use the command "wmake" for compilation of dependency file.
error appears
please check i think i have some errors in declaration

Thank You