[nakor]

Hey,
I am trying to simulate the heat transfer from a solid to a fluid.

I am using Gambit to generate my geometry, which is a tube with a solid inside.
Then I use fluent3DMeshToFoam :

Dimension of grid: 3
Number of points: 18170
PointGroup: 1 start: 0 end: 18169. Reading points...done.
Number of faces: 34860
FaceGroup: 4 start: 0 end: 19. Reading mixed faces...done.
FaceGroup: 5 start: 20 end: 17203. Reading mixed faces...done.
FaceGroup: 6 start: 17204 end: 17227. Reading mixed faces...done.
FaceGroup: 7 start: 17228 end: 17243. Reading mixed faces...done.
FaceGroup: 8 start: 17244 end: 18187. Reading mixed faces...done.
FaceGroup: 10 start: 18188 end: 34859. Reading mixed faces...done.
Number of cells: 8592
CellGroup: 2 start: 0 end: 8566 type: 1
CellGroup: 3 start: 8567 end: 8591 type: 1
Zone: 2 name: v_fluid type: fluid. Reading zone data...done.
Zone: 3 name: v_solid type: solid. Reading zone data...done.
Zone: 4 name: wall_internal type: wall. Reading zone data...done.
Zone: 5 name: sym type: symmetry. Reading zone data...done.
Zone: 6 name: outlet type: pressure-outlet. Reading zone data...done.
Zone: 7 name: inlet type: velocity-inlet. Reading zone data...done.
Zone: 8 name: wall type: wall. Reading zone data...done.
Zone: 10 name: default-interior type: interior. Reading zone data...done.

FINISHED LEXING

Creating patch 0 for zone: 4 name: wall_internal type: wall
Creating patch 1 for zone: 5 name: sym type: symmetry
Creating patch 2 for zone: 6 name: outlet type: pressure-outlet
Creating patch 3 for zone: 7 name: inlet type: velocity-inlet
Creating patch 4 for zone: 8 name: wall type: wall
Creating cellZone 0 name: v_fluid type: fluid
Creating cellZone 1 name: v_solid type: solid
Creating faceZone 0 name: default-interior type: interior
faceZone from Fluent indices: 18188 to: 34859 type: interior
patch 0 from Fluent indices: 0 to: 19 type: wall
patch 1 from Fluent indices: 20 to: 17203 type: symmetry
patch 2 from Fluent indices: 17204 to: 17227 type: pressure-outlet
patch 3 from Fluent indices: 17228 to: 17243 type: velocity-inlet
patch 4 from Fluent indices: 17244 to: 18187 type: wall

Then I have to split my mesh into two different part v_solid and v_fluid, but I don't know how to do that.
Should I use something like splitMesh ?

Thank you for your help !

[Aurelien Thinat]

Hello Nakor,

I don't remember right now the different steps you have to follow in order to set up a case for chtMultiRegionFoam, but there is already a big thread on this solver in the "running" part of this forum. I found all I needed in it.

Edit : Did you add a "regionProperties" file in the constant folder before using the fluent3DMeshToFoam ?

[nakor]

yes, I put a "region properties" file before using fluent3DMeshToFOam.

Then I try to use "splitMeshRegions -cellZones" and my output is :

Number of regions:2

Writing region per cell file (for manual decomposition) to "/home/users/romain/OpenFOAM/romain-1.7.x/run/backstep/V01-withcube/kw-SST-Groovy/constant/cellToRegion"

Writing region per cell as volScalarField to "/home/users/romain/OpenFOAM/romain-1.7.x/run/backstep/V01-withcube/kw-SST-Groovy/0/cellToRegion"

Region Cells
------ -----
0 8567
1 25

Region Zone Name
------ ---- ----
0 0 v_fluid
1 1 v_solid

Sizes inbetween regions:

Region Region Faces
------ ------ -----

Reading volScalarField cellToRegion
Reading volScalarField p
Reading volScalarField nut
Reading volScalarField k
Reading volScalarField omega
Reading volScalarField nuTilda
Reading volVectorField U
Reading volSymmTensorField R


Adding patches


Adding patches


Region 0
--------
Creating mesh for region 0 v_fluid
Testing:"/home/users/romain/OpenFOAM/romain-1.7.x/run/backstep/V01-withcube/kw-SST-Groovy/system/v_fluid/fvSchemes"
Writing dummy "v_fluid/fvSchemes"
Writing dummy "v_fluid/fvSolution"
Mapping fields
Mapping field cellToRegion
Mapping field k
Mapping field omega
Mapping field nuTilda
Mapping field nut
Mapping field p
Mapping field U
Mapping field R
Deleting empty patches
Writing new mesh
Writing addressing to base mesh
Writing map pointRegionAddressing from region0 points back to base mesh.
Writing map faceRegionAddressing from region0 faces back to base mesh.
Writing map cellRegionAddressing from region0 cells back to base mesh.

Region 1
--------
Creating mesh for region 1 v_solid
Testing:"/home/users/romain/OpenFOAM/romain-1.7.x/run/backstep/V01-withcube/kw-SST-Groovy/system/v_solid/fvSchemes"
Writing dummy "v_solid/fvSchemes"
Writing dummy "v_solid/fvSolution"
Mapping fields
Mapping field cellToRegion
Mapping field k
Mapping field omega
Mapping field nuTilda
Mapping field nut
Mapping field p
Mapping field U
Mapping field R
Deleting empty patches
Writing new mesh
Writing addressing to base mesh
Writing map pointRegionAddressing from region1 points back to base mesh.
Writing map faceRegionAddressing from region1 faces back to base mesh.
Writing map cellRegionAddressing from region1 cells back to base mesh.
End

But it did not create the folder v_solid and v_fluid in constant/ thus when I can not put my file thermophysicalproperties.
But it did create a new folder "mycase/1/" in which I have 2 folders "v_fluid" and "v_solid".
Should I put those two folders in "constant/" ?

[Aurelien Thinat]

I took a look at a case I have already setp up.

I also have the folder 0/ and 1/. It seems I used this command too and I copy-paste the folder 1/ files in the 0/ (fluid and wall on my case).

Then, in the constant/fluid/ folder, I have RASProperties, thermophisicalProperties...as in any constant folder for a single fluid region solver. So you can copy paste the the folders v_fluid v_solid in your constant/ folder and add the files needed by hand in each of the folders.

[nakor]

I did as you said, copy paste the the folders v_fluid v_solid in my constant/ folder and add the files needed by hand in each of the folders (cp, K,p_rgh,rho, T, g, thermophysical.., turbulenceproperties).

But it look like I have a problem with a library.

Create time

Create fluid mesh for region v_fluid for time = 0

Create solid mesh for region v_solid for time = 0

*** Reading fluid mesh thermophysical properties for region v_fluid

Adding to thermoFluid

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
#0 Foam::error:: printStack(Foam::Ostream&) in "/usr/local/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/usr/local/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::hPsiThermo<Foam:: pureMixture<Foam::constTransport<Foam::specieTherm o<Foam::hConstThermo<Foam:erfectGas> > > > >::calculate() in "/usr/local/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:: pureMixture<Foam::constTransport<Foam::specieTherm o<Foam::hConstThermo<Foam:erfectGas> > > > >::hPsiThermo(Foam::fvMesh const&) in "/usr/local/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#5 Foam::basicPsiThermo::addfvMeshConstructorToTable< Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > > >::New(Foam::fvMesh const&) in "/usr/local/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#6 Foam::basicPsiThermo::New(Foam::fvMesh const&) in "/usr/local/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#7 main in "/usr/local/OpenFOAM/OpenFOAM-1.7.x/applications/bin/linux64GccDPOpt/chtMultiRegionSimpleFoam"
#8 __libc_start_main in "/lib/libc.so.6"
#9 _start at /build/buildd/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Exception en point flottant

It is a bit stange because everything is working fine with the tutorial case (in "tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/"), but it looks like I am missing a file when I am running it in my case.

[Aurelien Thinat]

Which files do you have in your constant/v_fluid/ folders ?

[nakor]

I have the folder polymesh, and g, RASproperties,thermophysicalproperties, turbulenceporperties, cellToregions, (and k, nut, nutilda, omega, p, R, U, but I don't think they are really usefull here).

[Aurelien Thinat]

"(and k, nut, nutilda, omega, p, R, U, but I don't think they are really usefull here)."

Yes they are useless.


Can you copy paste the thermophysicalProperties file ?

[nakor]

FoamFile
{
version 2.0;
format ascii;
class dictionary;
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>;

mixture air 1 28.9 1000 0 1.8e-05 0.7;

I just took the one from the tutorials without modification.

[Aurelien Thinat]

Well, I don't know where the problem come from.

Do you add the v_fluid/ and v_solid/ in system/ ?
My case follows exactly what you have and it's runnable.

[nakor]

yes I have those two folders in system.
In each of them I have changeDictonaryDict, decomposeParDict, fvSchemes, and fvSolution.

[sixwp]

Hi nakor,

Did you manage to solve your issue?

Hi All

Any hint?

Best regards