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

conjugateHeatFoam error in new commits of 1.6-ext

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   June 22, 2011, 04:41
Default
  #21
Member
 
cosimo bianchini
Join Date: Mar 2009
Location: Florence, Tuscany, Italy
Posts: 88
Rep Power: 9
cosimobianchini is on a distinguished road
Send a message via Skype™ to cosimobianchini
I incurred in an analogous problem some time ago (same problem but in my own conjugate heat transfer solver: attached and detached modes are managed changing the specified boundary condition at the interface).
I solved it modifying the interface boundary condition in order to maintain a constant heat flux instead of a constant temperature gradient also in case of non-continuous conductivity.
Basically for the coupled interface both the internal and boundary coefficients were modified to guarantee continuity of heat flux and temperature on the interface itself.
1D conduction case was conformal to analytic solution to machine precision. Some more details of the implemented boundary condition may be find in the following paper presented at the open source CFD international conference 2008 http://brun.de.unifi.it/docpub/OF_Co..._C3X_final.doc
__________________
Cosimo Bianchini

Ergon Research s.r.l.
Via Panciatichi, 92
50127 Florence - ITALY
Tel: +39 055 0763716
Mob: +39 320 9460153
e-mail: cosimo.bianchini@ergonresearch.it
URL: www.ergonresearch.it
cosimobianchini is offline   Reply With Quote

Old   June 22, 2011, 06:11
Default new test
  #22
Senior Member
 
niaz's Avatar
 
A_R
Join Date: Jun 2009
Posts: 118
Rep Power: 9
niaz is on a distinguished road
Dear Hrvoje
1-I checked Benk test case. My left diffusivity is 10-5 and right is 10-3: on the interface has diffusivity 1.9802*10-5 for both sides. It shows that DT in boundary has updated.
2- For second test, I have used fixedGradient as BC. So, I should give the same Gradient value on the interface. But when I solved it, I faced to new problem. The value for flux on the interface is not similar!!!!
Attached Images
File Type: jpg babak.jpg (11.9 KB, 27 views)
niaz is offline   Reply With Quote

Old   June 22, 2011, 12:32
Default
  #23
Senior Member
 
Ben K
Join Date: Feb 2010
Location: Ottawa, Canada
Posts: 140
Rep Power: 11
benk is on a distinguished road
For the simple model that I made with a jump in diffusivity from 1e-5 to 1e-3, I calculated the fluxes by hand:
https://spreadsheets.google.com/spre...thkey=CMCd18IP

The fluxes match exactly if you use harmonic (or localMin/localMax) in laplaceSchemes but from my calculations, it looks like using Gauss linear corrected does have a problem with the fluxes through the interface.

It's also a good idea to refine the mesh close to the interface for large jumps in diffusivity like this though. I haven't tried this.

kamkari: To calc gradients in detatched state, you could just try creating another submesh which spans all regions and map the solution to this new submesh, then you can call snGrad on fields within this new mesh without needing to attachpatches.
benk is offline   Reply With Quote

Old   June 23, 2011, 03:38
Default
  #24
New Member
 
babak kamkari
Join Date: Dec 2010
Posts: 26
Rep Power: 7
kamkari is on a distinguished road
Dear Professor Jasak

Thank you for your explanation and sorry for annoying you. As you recommended I went to diffusivity filled. the values of updated diffusivities were the same on the both side of the interface as it was expected. Here the main problem is that the gradient of concentration (or temperature) on the sides of the interface are not correct using both linerMax and linear scheme for laplacian interpolation. I am working on your comments and Dr. Kenney results.
I also should add, it is my great honor to take part in your classes and you fail me because I will find another opportunity for learn more.

Thank for your time and patience
kamkari is offline   Reply With Quote

Old   July 11, 2011, 21:11
Default
  #25
New Member
 
Join Date: Oct 2010
Posts: 27
Rep Power: 7
albem is on a distinguished road
Hello

I am trying to run conjugateHeatFoam solver for the conjugateCavity problem, the program starts running but, after the first time calculation it crashed, I am using OpenFOAM 1.6-ext mounted in MacOSX 10.6.7

I followed the useful link from Hrv to instal OpenFoam

Extend Project 1.6-ext Binary Release for Mac OS X

Any help is welcome

thanks in advance


Starting time loop

Time = 0.01

Courant Number mean: 0 max: 0 velocity magnitude: 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.36284e-06, No Iterations 6
DILUPBiCG: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for p, Initial residual = 1, Final residual = 7.95294e-07, No Iterations 18
time step continuity errors : sum local = 1.41709e-08, global = -6.77626e-19, cumulative = -6.77626e-19
DICPCG: Solving for p, Initial residual = 0.431067, Final residual = 1.96388e-07, No Iterations 18
time step continuity errors : sum local = 6.11905e-09, global = 6.77626e-19, cumulative = 0


--> FOAM FATAL IO ERROR:
keyword preconditioner is undefined in dictionary ":reconditioner"

file: :reconditioner from line 39 to line 39.

From function dictionary::lookupEntry(const word&, bool, bool) const
in file db/dictionary/dictionary.C at line 396.

FOAM exiting




I don't get this error when using the binary mac version of 1.6-ext but I do when I use a newer commit of 1.6-ext. For the version that gives me a problem, git log gives:

commit 2f89ab8b718f44bda09af0ed421ec6751d7c52fa
Date: Fri May 27 20:53:20 2011 +0100

For the version that doesn't give me a problem, git log gives:

commit 3b1cb74951fe1630a2b3f39fc8d267622421f54d
Date: Fri Nov 26 11:08:01 2010 + 0000[/QUOTE]
albem is offline   Reply With Quote

Old   July 11, 2011, 22:00
Default
  #26
Senior Member
 
Ben K
Join Date: Feb 2010
Location: Ottawa, Canada
Posts: 140
Rep Power: 11
benk is on a distinguished road
Quote:
Originally Posted by albem View Post
--> FOAM FATAL IO ERROR:
keyword preconditioner is undefined in dictionary ":reconditioner"

file: :reconditioner from line 39 to line 39.
It sounds like a problem with your fvSolution file. I specify the preconditioner like this:

Code:
    T+T
    {
        solver           BiCGStab;
        preconditioner Cholesky;
        minIter      0;
        maxIter      1000;
        tolerance    1e-20;
        relTol       0.0;
    };
benk is offline   Reply With Quote

Old   September 20, 2011, 18:26
Default
  #27
New Member
 
Join Date: Oct 2010
Posts: 27
Rep Power: 7
albem is on a distinguished road
Quote:
Originally Posted by benk View Post
When I try to run the conjugateHeatFoam tutorial in 1.6-ext, I get the following error message in log.conjugateHeatFoam:



I don't get this error when using the binary mac version of 1.6-ext but I do when I use a newer commit of 1.6-ext. For the version that gives me a problem, git log gives:

commit 2f89ab8b718f44bda09af0ed421ec6751d7c52fa
Date: Fri May 27 20:53:20 2011 +0100

For the version that doesn't give me a problem, git log gives:

commit 3b1cb74951fe1630a2b3f39fc8d267622421f54d
Date: Fri Nov 26 11:08:01 2010 + 0000

Hi

I have installed OpenFOAM-1.6-ext on my Mac, I could compile conjugateHeatFoam solver but when I ran it I got

--> FOAM FATAL ERROR:
Attempt to cast type wall to type lduInterface

From function refCast<To>(From&)
in file /Users/hjasak/OpenFOAM/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude/typeInfo.H at line 115.


It is strange to me because the directory /Users/hjasak/OpenFOAM/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude

does not exist in my Mac

Any help is welcome

Alberto
albem is offline   Reply With Quote

Old   May 29, 2012, 05:51
Default
  #28
New Member
 
Join Date: May 2012
Posts: 2
Rep Power: 0
Rakenta is on a distinguished road
Hi to all!
I am new to OpenFOAM and I try to run the conjugateHeatFoam tutorial in 1.6-ext according to the tutorial of Johan Magnusson (1.5-dev).
Running the blockMesh for conjugateCavity and heatedBlock works fine, but whenever I try to run conjugateHeatFoam I get the following error message:

--> FOAM FATAL ERROR:
Cannot find file "points" in directory "constant/solid/polyMesh"

From function Time::findInstance(const fileName&, const word&, const IOobject::readOption)
in file db/Time/findInstance.C at line 148.


As I tried it for a long time to fix the problem without any success, I really would appreciate if you can give me any hints what I am doing wrong.
Thanks in advance
Rakenta is offline   Reply With Quote

Old   June 6, 2012, 06:09
Default
  #29
New Member
 
Join Date: May 2012
Posts: 2
Rep Power: 0
Rakenta is on a distinguished road
Well, I solved it. Stupid mistake.
You have to create the links of the folders 0, constant and system in the heated block, rename them to solid and shift them into the folders 0, constant and system of conjugateCavity.
Rakenta is offline   Reply With Quote

Old   June 18, 2013, 14:24
Default
  #30
Senior Member
 
mahdi abdollahzadeh
Join Date: Mar 2011
Location: Covilha,Portugal
Posts: 150
Rep Power: 7
mm.abdollahzadeh is on a distinguished road
Hi all

I want to modify the regioncouplying .

my boundary conditions is somehow"

k1grad(T)=k2grad(T)+sigma.

for this case I am already reading the the sigma field in regioncoupling. I think I need to access delta() of the coupled patch.

Any idea how to do it?

best
mahdi
mm.abdollahzadeh 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
Test directory missing in OpenFOAM 1.6 ext andrewryan OpenFOAM 2 March 20, 2011 16:41


All times are GMT -4. The time now is 15:58.