Error in thermophysical properties (chtMultiRegionFoam)
 

Hello foamers!

I have modified multiRegionHeater Tutorial of chtMultiRegionFoam to following geometry, as I want to simulate plasma actuator induced flow:
http://i40.tinypic.com/2j61uag.jpg

In tutorials, following regions were created by toposetDict:

topAir
bottomAir
heater
leftSolid
rightSolid

Then I have modified the geometry of this tutorial like as above image in which I have defined following regions by topoSetDict

topAir (same as tutorial but differ in dimension)
bottomAir (same as tutorial but differ in dimension)
heater (as a dielectric material which is kapton film)
leftSolid (as a top electrode which is copper)
rightSolid (as a bottom electrode which is copper)
innerelec (as a inner electrode which is also copper)

I have set themophysical properties for each region according to their material properties except topAir and bottomAir as they are same air regions in my case like tutorials.

I have run following commands in the terminal:

But I found following error after executing
mpirun -np 4 chtMultiRegionFoam -parallel


How can I get rid of it?

Thanks in advance.

Best regards,
Mukut

Did you make your own library?
libfluidThermophysicalModels.so?
The standard should just be basicThermophysicalModels.
If you did make your own library, you have to make sure that everything is connected and pointed to properly.

The error is referring to the calculate function in the sensible enthalpy section so it may also be a typo there.

If all you did was change the geometry by adding more solids, you shouldn't be getting any of those kinds of errors since they relate to the solver itself not the problem.
Granted that is from 2.1.1, but I doubt they would add fluid just like that.

Thank you so much. the problem was solved.. :)

foam::error::printStack compressibleInterFoam parallel run
 

Hi mukut

How did you solve the problem with the foam::error::printStack? I am facing the same problem like you did in OpenFOAM 2.2.2. I perform simulations with compressibleInterFoam on a case with round about 2e6 cells. In serial the simulation runs perfect but would take weeks to finish. In parallel run I get the following error massage (I cleaned the full error output to show only errors from one calculation node):


The error occurs when OpenFOAM tries to initialize the thermophysical model of the air phase. I only changed the geometry from the tutorial case. Moreover, the error shows a problem of the thermophysical model, not the case itself. Thanks for your help.


Cheers


Fabian

twoPhaseMixtureThermo bug?
 

Hi

I started from scratch to find the trigger of my error in compressibleInterFoam case. I took the tutorial case and ran it serial and parallel. Everything is fine. When I change the temperatures to the range I have to simulate my specific case (round about 900 °C), then the serial case runs flawless. However the parallel run stops with printStack error after constructing the twoPhaseMixtureThermo. So the only problem is the increase in temperature? I also tried with janaf thermo model without success. Could this be a bug (OF 2.2.2)? Help is still welcome!

Cheers

Fabian

hello,

Did you check if after decomposePar, you did not get in the proc../0/T files some "value (non)uniform 0;" ?
Because some tools (merge / stitch patch, maybe other ...) seem to do that (i don"t know why).

regards,
olivier

Hi

Sorry for this huge delay. I think you have a good point here. I used snappy to generate the mesh. To not have to reconstruct I copied the 0 time folder into the processor directories after meshing. Then I ran setFields. The important point know: I use

in my fields. So I do not set the processor boundaries myself. This could lead to the error. I will check this an confirm the solution here asap.

Cheers

Fabian

Dear Mukut or anyone else who might be able to help,

I am trying to use chtMultiRegionFoam to model transient conduction between two solids (no fluid). I am having a very difficult time and I hope you don't mind that I trouble you with such a basic problem. Is there an possibility that you could have a look at my two-solid heat transfer case and see where I have gone wrong?

https://www.dropbox.com/sh/1fxb42g1b...92kKyZgha?dl=0

I currently do not have any fluids in the region properties. I suppose this could cause trouble for chtMultiRegion.

Thank you,
Cliff

Cliff I would suggest changing the region properties file. You haven't changed it. You are telling it that there are fluid regions involved, but then not providing them.

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object regionProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

fluidRegionNames ( );

solidRegionNames ( solid1 solid2);

// ************************************************** *********************** //

is a good place to start. Then you have to make sure that you have the appropriate folders in your system folder for the regions, including the changeDictionaryDict files in their respective system/regionName folder

Thank you very much. It's running now. Much joy here.

You wouldn't happen to know why in the controlDict I am not able to successfully run with

adjustTimeStep yes;

when I remove this, it runs, albeit slowly. The adjustable time step works well for the multiRegionHeater tutorial case, but not in my simple case.

Cliff

Hello everyone,
I am trying to solve conjugate heat transfer problem in chtMultiRegionFoam with geometry of cavity with some insulation. That means two region...solid and air. But it is showing some problem when it read thermophysical property..please suggest something..

https://www.dropbox.com/s/kdt1prx29c...heck1.PNG?dl=0

-Baran

nonuniform value 0
 

Hi

I usually get those errors when the initial fields and boundary conditions of pressure p and temperature T have nonphysical values. Sometimes, for example when you decompose a case for a parallel run, the fields get initialized with nonuniform values. So when you find a pressure or temperature boundary with value 0 this causes the error. Most obvious, the thermodynamic model can't handle a fluid with 0 K temperature and 0 Pa pressure. Best way to get around this is using changeDictionary and overwrite all fields before running.

Cheers

Fabian

Very very thanks for your guidance...Its working for this now. But for my case solver is not reading "changeDictionary"... I manually change bc in 0 directory...but it is showing problem regarding boundary condition of T..
https://www.dropbox.com/s/s8faeevvgd...heck3.PNG?dl=0

After that I change bc manually in 0 directory to
After that it is showing following error..
https://www.dropbox.com/s/3vs6c1h64y...heck2.PNG?dl=0

Maybe this is very basic thing..but...as I am new to OF...so struggling to understand this thing...please suggest ...

-baran

Greetings to all!

@baran: I saw the PM you sent me and came looking for public posts you had. I'm glad the first issue is solved.

But the new problem is because you're using the boundary condition type "calculated" where it should not be used.
Keep in mind that the "calculated" type of boundary condition means that the result on that patch for that specific field, is a result of the flow on the rest of the domain and that it can be calculated based on the other fields or from the flow itself. In other words, the error message you're getting is implying that you cannot use the "calculated" type for your patch.

Beyond this, I suggest that you take a few steps back and test ideas using a simpler test case. For example, you can use the case given here: http://openfoamwiki.net/index.php/Ge..._-_planeWall2D - since it's a simpler case than yours, and since it's one where it's easier to calculate the analytical solution for comparison, you can then make small changes to make it closer and closer to your problem at hand. For example, this way you can test various thermodynamic properties with this test case, to assess what you can and cannot do.
In other words: isolate and conquer one small problem at a time ;)

Best regards,
Bruno

Greeting everyone,

@bruno.....thanks for your guidance....I went through the tutorial of multiRegionHeater of chtMultiRegionFoam. According to this I modified my case as well as boundary conditions which was implemented in wrong way initially. But still there is some problem in turbulence model which is observed in log.chtMultiRegionFoam though I used same thing like in tutorial .

I am unable to sort it out this issue.....!!! Please suggest something...

Regards,
Baran

Hi baren

probably you found out already that changeDictionary has to be called before starting the solution process with the solver. You have to prepare a changeDictionaryDict for all your regions. And then execute the change with the command

in the shell. The error you get from the solver concerning you turbulence model could be due to a wrong thermoType. Check this thread:

http://www.cfd-online.com/Forums/ope...tml#post381177

Took me 1 minute to Google and read :(

Regards Fabian

Hello everyone,
@ fabian , thanks...for your suggestion...i sort it out the problem......of thermophysical properties...

but there is some mpi running error is coming for my case...its not running..
then i go through this...
http://webcache.googleusercontent.co...&ct=clnk&gl=in
and change my deltaT to 10^-8 ....but still it showing same error....please suggest something...

Regards,
Baran

Hello everyone,
my problem is solved....I did the following mistakes...

1. I convert one gambit mesh file into openfoam ,so its scaling is forgot to change..
2. Adjustable run time is on...so deltaT is changed to some different value...which create the floating point exception issue...

Regards,
Baran

Error in chtMultiRegionFOAM
 

Hi everyone!

I have this error during the execution of my simulation:

Solving for fluid region water
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.167835e-009, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 5.9374103e-008, No Iterations 4


--> FOAM FATAL ERROR:
Kappa defined to employ fluidThermo method, but thermo package not available

From function temperatureCoupledBase::kappa(const scalarField&) const
in file derivedFvPatchFields/temperatureCoupledBase/temperatureCoupledBase.C at line 116.

FOAM exiting

The solver is ChtmultiRegionFOAM.

Anyone can help me?

Thank you so much!

Quick questions @fkika: Please provide more details, such as:

1. Which OpenFOAM version are you using?
2. Which installation instructions did you follow for installing OpenFOAM?
3. Can you reproduce this error with one of OpenFOAM's tutorials?
4. How is the content of the files "thermophysicalProperties" in your case? If you don't know where they are, run:
   
   Code:

   ---

   find . -name thermophysicalProperties

   ---