CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

solid to fluid heattransfer with chtMultiRegionFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 24, 2010, 13:07
Default solid to fluid heattransfer with chtMultiRegionFoam
  #1
New Member
 
Romain
Join Date: Jun 2010
Location: Lyon
Posts: 28
Rep Power: 16
nakor is on a distinguished road
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 :

Quote:
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 !
nakor is offline   Reply With Quote

Old   August 25, 2010, 04:37
Default
  #2
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 16
Aurelien Thinat is on a distinguished road
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 ?
Aurelien Thinat is offline   Reply With Quote

Old   August 25, 2010, 05:35
Default
  #3
New Member
 
Romain
Join Date: Jun 2010
Location: Lyon
Posts: 28
Rep Power: 16
nakor is on a distinguished road
yes, I put a "region properties" file before using fluent3DMeshToFOam.

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

Quote:
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/" ?
nakor is offline   Reply With Quote

Old   August 25, 2010, 05:50
Default
  #4
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 16
Aurelien Thinat is on a distinguished road
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.
Aurelien Thinat is offline   Reply With Quote

Old   August 25, 2010, 06:51
Default
  #5
New Member
 
Romain
Join Date: Jun 2010
Location: Lyon
Posts: 28
Rep Power: 16
nakor is on a distinguished road
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.
Quote:
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.
nakor is offline   Reply With Quote

Old   August 25, 2010, 06:53
Default
  #6
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 16
Aurelien Thinat is on a distinguished road
Which files do you have in your constant/v_fluid/ folders ?
Aurelien Thinat is offline   Reply With Quote

Old   August 25, 2010, 06:56
Default
  #7
New Member
 
Romain
Join Date: Jun 2010
Location: Lyon
Posts: 28
Rep Power: 16
nakor is on a distinguished road
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).
nakor is offline   Reply With Quote

Old   August 25, 2010, 07:30
Default
  #8
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 16
Aurelien Thinat is on a distinguished road
"(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 ?
Aurelien Thinat is offline   Reply With Quote

Old   August 25, 2010, 08:39
Default
  #9
New Member
 
Romain
Join Date: Jun 2010
Location: Lyon
Posts: 28
Rep Power: 16
nakor is on a distinguished road
Quote:
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.
nakor is offline   Reply With Quote

Old   August 25, 2010, 08:49
Default
  #10
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 16
Aurelien Thinat is on a distinguished road
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.
Aurelien Thinat is offline   Reply With Quote

Old   August 25, 2010, 08:53
Default
  #11
New Member
 
Romain
Join Date: Jun 2010
Location: Lyon
Posts: 28
Rep Power: 16
nakor is on a distinguished road
yes I have those two folders in system.
In each of them I have changeDictonaryDict, decomposeParDict, fvSchemes, and fvSolution.
nakor is offline   Reply With Quote

Old   March 21, 2011, 09:28
Default
  #12
Member
 
Samuel ARNAUD
Join Date: Feb 2011
Location: Grenoble, FRANCE
Posts: 39
Rep Power: 15
sixwp is on a distinguished road
Hi nakor,

Did you manage to solve your issue?

Hi All

Any hint?

Best regards
__________________
Sam
sixwp is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Can CFX do CHT simulations with a solid domain rotating in a stationary fluid domain? acro CFX 15 September 23, 2016 12:16
fluid to solid heattransfer with chtMultiRegionFoam schteff OpenFOAM 5 August 20, 2010 08:45
Theory behind the solid domain rotation in fluid using Bets property Riken CFX 0 October 28, 2009 10:51
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11
CFX4.3 -build analysis form Chie Min CFX 5 July 13, 2001 00:19


All times are GMT -4. The time now is 00:48.