CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

chtIcoMultiRegionFoam - Incompressible version of chtMultiRegionFoam.

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

Like Tree5Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   November 5, 2010, 08:25
Default chtIcoMultiRegionFoam - Incompressible version of chtMultiRegionFoam.
  #1
Member
 
Fábio César Canesin
Join Date: Mar 2010
Location: Florianópolis
Posts: 67
Rep Power: 7
Canesin is on a distinguished road
Hi all,

I have been working with some modified chtMultiRegion codes for some time.. This version is the same I have send to Hrvoje Jasak for inclusion on the next version OpenFOAM-1.6-ext
To complie just enter the folder and type "wmake".

The solver was tested without turbulence, I will add an test case to this thread as soon as I have time to setup an easy to understand..

Please, keep the references in the source-code.. If you like to reference this solver you can:

CANESIN, F. C. : chtIcoMultiRegionFoam, incompressible multi-region fully segregated conjugated heat transfer - http://www.canesin.com/software

I hope to see good use of the code, I'm working in an fully coupled version and in some addons and a test case for my undergraduate thesis .. but don't hope for it so soon.

Best regards,
Fábio C. Canesin
Attached Files
File Type: gz chtIcoMultiRegionFoam.tar.gz (18.1 KB, 499 views)
awacs and ramon like this.
Canesin is offline   Reply With Quote

Old   November 5, 2010, 10:39
Default
  #2
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 6
stevenvanharen is on a distinguished road
Hi Fabio,

I have been working on the exact same thing.

It is on the shelf for now but beginning next year I will do some work with my solver.

Probably I will start with validating against work done by Tiselj (Conjugate heat transfer in channel flow) since I will do fully turbulent DNS.

What are you working on?

Regards,

Steven
stevenvanharen is offline   Reply With Quote

Old   November 5, 2010, 14:43
Default
  #3
Member
 
Fábio César Canesin
Join Date: Mar 2010
Location: Florianópolis
Posts: 67
Rep Power: 7
Canesin is on a distinguished road
Hi Steven,

I have been working with geometric optimization for active magnetic regenerators .. it is an setup for magnetic cooling.
Canesin is offline   Reply With Quote

Old   November 5, 2010, 17:24
Default Does it works for unstructured meshes?
  #4
Member
 
Pablo Caron
Join Date: Nov 2009
Location: Buenos Aires, Argentina
Posts: 73
Rep Power: 7
pcaron is on a distinguished road
Hi Fábio,

I've using chtMultiReginFoam, from 1.6.x, for some weeks. Apparently there is an issue with unstructured meshes, or I'm missing something . If you set up a case with uniform temperatures, spurious velocities arises.

Have you tried using unstructured meshes?

Best regards,

Pablo
pcaron is offline   Reply With Quote

Old   November 5, 2010, 20:27
Default
  #5
Member
 
Fábio César Canesin
Join Date: Mar 2010
Location: Florianópolis
Posts: 67
Rep Power: 7
Canesin is on a distinguished road
Hi Pablo,

I have not.. but, I know that in my solver will not have this issue, because it is a fully segregated approach, it solves velocity and them uses it to transport temperature..

The case in chtMultiRegionFoam cam be from the Boussinesq aproximation, where the temperature field is diverging the pressure.. or maybe some wrong boundary condition..

Best regards,
Fábio C. Canesin
Canesin is offline   Reply With Quote

Old   November 7, 2010, 16:30
Default
  #6
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 271
Rep Power: 9
phsieh2005 is on a distinguished road
Hi, Canesin,

I will try to see if I can modify the tutorial case in OF-1.7.x to run with your chtIcoMultiRegionFoam solver. Any suggestion before I start?

What is the reason of developing the chtIcoMultiRegionFoam over chtMultiRegionFoam (which uses compressible solver for fluid I believe)?

Pei
phsieh2005 is offline   Reply With Quote

Old   November 8, 2010, 05:26
Default
  #7
Member
 
Fábio César Canesin
Join Date: Mar 2010
Location: Florianópolis
Posts: 67
Rep Power: 7
Canesin is on a distinguished road
Quote:
Originally Posted by phsieh2005 View Post
Hi, Canesin,

I will try to see if I can modify the tutorial case in OF-1.7.x to run with your chtIcoMultiRegionFoam solver. Any suggestion before I start?
Yes, you can setup the case where the fluid has a configuration as like pimpleFoam and the solid as like the solid in chtMultiRegionFoam

Quote:
Originally Posted by phsieh2005 View Post
What is the reason of developing the chtIcoMultiRegionFoam over chtMultiRegionFoam (which uses compressible solver for fluid I believe)?

Pei
Besides using an compressible solver in chtMultiRegionFoam it can be used to simulate water, as is done in the new tutorial in 1.7.x, but it is like as little mess, you have to define and compressible::incompressible thermophysic property O.o.... Also it uses Boussinesq aproximation for natural(free) convection.. In that solver I propose there is no effects from the temperature in the fluid, like in forced convection.
Canesin is offline   Reply With Quote

Old   November 8, 2010, 13:40
Default
  #8
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 271
Rep Power: 9
phsieh2005 is on a distinguished road
Hi, Fábio,

Is there any reason why K is required for the fluid region? What value should I use for K if fluid is water?

Thanks!

Pei
phsieh2005 is offline   Reply With Quote

Old   November 8, 2010, 16:56
Default
  #9
Member
 
Fábio César Canesin
Join Date: Mar 2010
Location: Florianópolis
Posts: 67
Rep Power: 7
Canesin is on a distinguished road
Quote:
Originally Posted by phsieh2005 View Post
Hi, Fábio,

Is there any reason why K is required for the fluid region? What value should I use for K if fluid is water?
K is the thermal conductivity of the fluid.. It is needed for the temperature coupling with the solid region, to compute the heat flux..

A tipical value for water is 0.6 W/m*k

Best regards,

Fábio C. Canesin
Canesin is offline   Reply With Quote

Old   November 8, 2010, 21:08
Default
  #10
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 271
Rep Power: 9
phsieh2005 is on a distinguished road
Hi, Fábio,

Thanks for the explanation.

This is more involved that what I originally expected. I basically used the fluid properties from buoyantBoussinesqPimpleFoam case and also made some changes to the fvSchemes and fvSolution. But, I am getting strange error messages. I will have to look into your code in more detail.

Pei
--------------------
Solving for fluid region bottomAir
DILUPBiCG: Solving for T, Initial residual = 1, Final residual = 4.8181319e-05, No Iterations 1
max(T) [0 0 0 1 0 0 0] 304.7137


--> FOAM FATAL ERROR:

request for uniformDimensionedVectorField g from objectRegistry bottomAir failed
available objects of type uniformDimensionedVectorField are

0
(
)


From function objectRegistry::lookupObject<Type>(const word&) const
in file /home/phsieh/OpenFOAM/OpenFOAM-1.7.x/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 139.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/phsieh/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/phsieh/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::UniformDimensionedField<Foam::Vector<double> > const& Foam:bjectRegistry::lookupObject<Foam::UniformDi mensionedField<Foam::Vector<double> > >(Foam::word const&) const in "/home/phsieh/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libfiniteVolume.so"
#3 Foam::buoyantPressureFvPatchScalarField::updateCoe ffs() in "/home/phsieh/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::updateCoef fs() in "/home/phsieh/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#5 Foam::adjustPhi(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/home/phsieh/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libfiniteVolume.so"
#6
in "/home/phsieh/OpenFOAM/phsieh-1.7.x/applications/bin/linux64GccDPOpt/chtIcoMultiRegionFoam"
#7 __libc_start_main in "/lib64/libc.so.6"
#8
at /usr/src/packages/BUILD/glibc-2.11.2/csu/../sysdeps/x86_64/elf/start.S:116
Aborted
phsieh@rutgers:~/OpenFOAM/phsieh-1.7.x/run/snappyIcoMultiRegionHeater>
phsieh2005 is offline   Reply With Quote

Old   November 8, 2010, 21:52
Default
  #11
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 271
Rep Power: 9
phsieh2005 is on a distinguished road
Hi, Fábio,

Sorry to bother you again. I am wondering which fluid solver your chtIcoMultiRegionFoam was based on? PimpleFoam? buoyantBoyssinesqPimpleFoam?

At a quick glance, I did not find the location in the code that g is given.

Also, I compiled your chtIcoMultiRegionFoam on OpenFOAM-1.7.x. I am wondering if this makes any difference. However, compilation was successful.

Are you planning to write any paper or thesis on your work? Maybe I can read your thesis/paper to figure it out?

Pei
phsieh2005 is offline   Reply With Quote

Old   November 10, 2010, 09:43
Default
  #12
Member
 
Fábio César Canesin
Join Date: Mar 2010
Location: Florianópolis
Posts: 67
Rep Power: 7
Canesin is on a distinguished road
The fluid is based in pimpleFoam .. ... this is not the solver for my work, this is a first generation, as the final version needs the not-released block-matrix solver for do pressure coupling and become fully implicit ..

It should not as for g in the solver .. in the controlDict you have changed the solver ??? O.o ..

If you can hold until the weekend I can setup a very basic case for it, and you can them use it to make your case work.. sorry, but I'm doing 53hours/week..

Best regards,

Fábio C. Canesin
Canesin is offline   Reply With Quote

Old   November 10, 2010, 10:34
Default
  #13
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 271
Rep Power: 9
phsieh2005 is on a distinguished road
Hi, Fábio,

I made some more changes. Now the case is running. This is the same case in the OF-1.7.x/snappyMultiRegionHeater. I will check if the results are reasonable when the run completes.

So, this solver does not handle natural convection? Is there any reason why you did not pick buoyantBousinessqPimpleFoam as fluid solver?

Pei
phsieh2005 is offline   Reply With Quote

Old   November 11, 2010, 17:55
Default
  #14
Member
 
Fábio César Canesin
Join Date: Mar 2010
Location: Florianópolis
Posts: 67
Rep Power: 7
Canesin is on a distinguished road
Quote:
Originally Posted by phsieh2005 View Post
Hi, Fábio,

I made some more changes. Now the case is running. This is the same case in the OF-1.7.x/snappyMultiRegionHeater. I will check if the results are reasonable when the run completes.

So, this solver does not handle natural convection? Is there any reason why you did not pick buoyantBousinessqPimpleFoam as fluid solver?

Pei
Yes it does not handle natural convection.. it was developed for forced convection cases and small channels cooling.. consistence and easy to expand, that where the reasons.
Canesin is offline   Reply With Quote

Old   November 13, 2010, 16:05
Default
  #15
Member
 
Pablo Caron
Join Date: Nov 2009
Location: Buenos Aires, Argentina
Posts: 73
Rep Power: 7
pcaron is on a distinguished road
Quote:
Originally Posted by Canesin View Post
Hi Pablo,

I have not.. but, I know that in my solver will not have this issue, because it is a fully segregated approach, it solves velocity and them uses it to transport temperature..

The case in chtMultiRegionFoam cam be from the Boussinesq aproximation, where the temperature field is diverging the pressure.. or maybe some wrong boundary condition..

Best regards,
Fábio C. Canesin
Dear Fábio, sorry for the late reply. It was a hard week. The problem I had seems to be related with the pressure definitios previous to 1.7 version.

You can see the release notes for OF http://www.openfoam.com/archive/1.7....ease-notes.php
Quote:
Modifications to multiphase and buoyant solvers
...
Multiphase and buoyant flow solvers now solve for p_rgh=p-rho g\cdot x, rather than the static pressure p. This change is to avoid deficiencies in the handling of the pressure force / buoyant force balance on non-orthogonal and distorted meshes.
...
So, I'll start using the standard new cht version. I've been very busy last week.

Regards

Pablo
pcaron is offline   Reply With Quote

Old   January 14, 2011, 03:38
Default
  #16
Senior Member
 
maddalena's Avatar
 
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 12
maddalena is on a distinguished road
Hi Fabio,
and thanks for sharing your solver! Incompressible & forced convection & unstructured mesh: this is what I was looking for! Have you made any tests with the turbulence on? I need it in my case. Can you tell me something about that?
Regards

mad

Last edited by maddalena; January 14, 2011 at 03:57.
maddalena is offline   Reply With Quote

Old   January 14, 2011, 13:15
Default
  #17
Member
 
Fábio César Canesin
Join Date: Mar 2010
Location: Florianópolis
Posts: 67
Rep Power: 7
Canesin is on a distinguished road
I have not done many tests in the turbulence settings.. But it is "FOAM way" of doing turbulence.. If you look at the source code you will see that turbulence is added as an term in the equations using turbulent prandtl and nu ..

You should be able to use any turbulence model from OpenFOAM, but it will be tied with your study case.
Canesin is offline   Reply With Quote

Old   January 17, 2011, 04:26
Default
  #18
Senior Member
 
maddalena's Avatar
 
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 12
maddalena is on a distinguished road
Quote:
Originally Posted by Canesin View Post
I have not done many tests in the turbulence settings.. But it is "FOAM way" of doing turbulence.. [...] You should be able to use any turbulence model from OpenFOAM, but it will be tied with your study case.
Sorry but I cannot understand... This is what usually applies to turbulence, isn't it?
One more question, open to everybody: is there anyone that has a steady state version of this incompressible multiregion cht solver?
Thank you

mad
maddalena is offline   Reply With Quote

Old   January 25, 2011, 09:54
Default
  #19
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 6
stevenvanharen is on a distinguished road
Hi all,

does somebody have a simple tutorial for this solver?

I have been trying to create a case myself but it keeps generating errors, I guess I am making a mistake somewhere in setting up the case.

It would be great if one of you guys could share a simple case you have performed with this solver.

Thanks in advance.

Kind regards,

Steven
stevenvanharen is offline   Reply With Quote

Old   January 25, 2011, 11:37
Default
  #20
Member
 
Fábio César Canesin
Join Date: Mar 2010
Location: Florianópolis
Posts: 67
Rep Power: 7
Canesin is on a distinguished road
Quote:
Originally Posted by maddalena View Post
Sorry but I cannot understand... This is what usually applies to turbulence, isn't it?
One more question, open to everybody: is there anyone that has a steady state version of this incompressible multiregion cht solver?
Thank you

mad
Yes, it is what applies to turbulence...

I do not have a steady state version because my problem do not have steady state solution.

But, you could use larger times steps.. Run fist the potentialFoam to have good fluids fields... tham use something like GAMG to make it more tolerant to instabilities.
Canesin is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
OpenFOAM on cluster: version GLIBCXX_3.4.9 and GLIBCXX_3.4.11 not found ovie OpenFOAM 9 October 15, 2012 01:59
paraview installation woes vex OpenFOAM Installation 15 January 30, 2011 08:11
bubbleFoam validation case balkrishna OpenFOAM Running, Solving & CFD 24 August 30, 2010 04:37
Problem with paraFoam on a linux-64 bit bunni OpenFOAM Paraview & paraFoam 4 April 14, 2010 20:55
paraFoam reader for OpenFOAM 1.6 smart OpenFOAM Installation 13 November 16, 2009 22:41


All times are GMT -4. The time now is 10:39.