CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Post-Processing (http://www.cfd-online.com/Forums/openfoam-post-processing/)
-   -   foamLog not solving for Ux, Uy, Uz (http://www.cfd-online.com/Forums/openfoam-post-processing/99947-foamlog-not-solving-ux-uy-uz.html)

aerospain April 16, 2012 15:56

foamLog not solving for Ux, Uy, Uz
 
Hi all,

I recently re-installed OpenFOAM due to a major cock-up after some installation gone wrong, and I am not getting the same behaviour from OF as I had before. I also lost my home folder with all the cases I had run for the last two months.

Right now, when I run simple, pimple or piso (FOAM) it doesn't "solve for Ux" (or Uy or Uz) So, I cannot plot them either after running foamLog on the log file created during the run.

I have been reviewing the tutorials to check for any hints in the fvSchemes and/or fvSolutions files without any luck.

Could someone give me a hand, please?

Thanks!

wyldckat April 16, 2012 16:12

Hi aerospain,

Something very strange is going on there. Does it not work with any tutorial at all? Not even the very first we learn about on the User Guide, the "incompressible/icoFoam/cavity"?

The only files that I can remember that might be impairing functionality would be located in the semi-hidden folder "~/.OpenFOAM". If you remove or rename that folder, you might be able to make things work as intended once again.

Best regards,
Bruno

aerospain April 16, 2012 18:51

Hello Bruno,

That .OpenFOAM hidden directory doesn't exist in my user's home folder. I have double-checked the cavity tutorial and the velocity components convergence history is reported during the simulation.

I may have messed up when looking at different files in the tutorials to build up my case. Let me describe it in case you could provide some hints:

incompressible k-epsilon simulation with one inlet (-x), one outlet (+x), symmetry on y (for now) with 'empty' BC, and top/bottom z with slip wall, symmetry, cyclic or other BC. At z+-0.5 there is a cylinder aligned axially with the flow, therefore 'wall' condition.

The base of this cylinder is 1m and my velocity is 0.1m/s. I will be starting with Re number below.

If worse comes to worst, is there a robust way to uninstall OpenFOAM from my Ubuntu 11.10 completely before I re-install it again.

Regards,
Carlos

wyldckat April 17, 2012 03:25

Hi Carlos,

Quote:

Originally Posted by aerospain (Post 355002)
If worse comes to worst, is there a robust way to uninstall OpenFOAM from my Ubuntu 11.10 completely before I re-install it again.

Without knowing how you installed it the first place, I don't know what's the best suggestion and I'm not in the mood to detail every single possible uninstallation procedure :(

Best regards,
Bruno

aerospain April 17, 2012 10:41

Hi Bruno,

Sorry for my badly made question from last night. I wasn't expecting anyone to 'waste' time explaining that issue. I forgot to mention that I used the 'sudo apt-get install ..." procedure for Ubuntu. I will try using "sudo apt-get remove ..." and see what happens.

I can foresee a wonderful afternoon of debugging ;-) I'll get me a nice cup of tea :-D

cheers!
C.

aerospain April 18, 2012 10:01

Hi Bruno,

I have solved my problem by rewriting the blockMeshDict from scratch.

I realized something was going wrong during mesh generation because not all the boundaries defined in my dictionary appeared as generated. Besides, the fact that I was getting velocity vectors through my body's wall had been puzzling ;-) (with non-slip walls)

Anyway, my solver is reporting velocity components to the log file.

I am attaching the faulty blockMeshDict file if someone want's to investigate.

cheers!
C.

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

convertToMeters 1;

vertices
(
(-2 -0.1 0.5) //0
( 0 -0.1 0.5) //1
(27 -0.1 0.5) //2
(-2 -0.1 11.5) //3
( 0 -0.1 11.5) //4
(27 -0.1 11.5) //5
(-2 -0.1 -0.5) //6
( 0 -0.1 -0.5) //7
(27 -0.1 -0.5) //8
(-2 -0.1 -11.5) //9
( 0 -0.1 -11.5) //10
(27 -0.1 -11.5) //11

(-2 0.1 0.5) //12
( 0 0.1 0.5) //13
(27 0.1 0.5) //14
(-2 0.1 11.5) //15
( 0 0.1 11.5) //16
(27 0.1 11.5) //17
(-2 0.1 -0.5) //18
( 0 0.1 -0.5) //19
(27 0.1 -0.5) //20
(-2 0.1 -11.5) //21
( 0 0.1 -11.5) //22
(27 0.1 -11.5) //23

( 0 -0.1 0) //24
(27 -0.1 0) //25
( 0 0.1 0) //26
(27 0.1 0) //27
);

blocks // Normals along positive x,y,z
(
/* 0 */ hex (ls 0 1 13 12 3 4 16 15) (20 1 110) simpleGrading (1 1 1)
/* 1 */ hex ( 1 2 14 13 4 5 17 16) (270 1 110) simpleGrading (1 1 1)
/* 2 */ hex (24 25 27 26 1 2 14 13) (270 1 5) simpleGrading (1 1 1)
/* 3 */ hex ( 7 8 20 19 24 25 27 26) (270 1 5) simpleGrading (1 1 1)
/* 4 */ hex (10 11 23 22 7 8 20 19) (270 1 110) simpleGrading (1 1 1)
/* 5 */ hex ( 9 10 22 21 6 7 19 18) (20 1 110) simpleGrading (1 1 1)
);

edges
(
);

boundary // Normals pointing out of domain
(
inlet
{
type patch;
faces
(
( 0 3 15 12) // Block 0
( 9 6 18 21) // Block 5
);
}
outlet
{
type patch;
faces
(
( 5 2 14 17) // Block 1
( 2 25 27 14) // Block 2
(25 8 20 27) // Block 3
( 8 11 23 20) // Block 4
);
}
skinWall
{
type wall;
faces
(
( 1 0 12 13) // Block 0
( 6 7 19 18) // Block 5
)
}
baseWall
{
type wall;
faces
(
( 24 1 13 26) // Block 2
( 24 26 19 7) // Block 3
)
}
outerDomainUp
{
type wall;
faces
(
( 3 4 16 15) // Block 0
( 4 5 17 16) // Block 1
);
}
outerDomainDown
{
type wall;
faces
(
(10 9 21 22) // Block 5
(11 10 22 23) // Block 4
);
}
leftAndRight
{
type empty;
faces
(
( 0 1 4 3) // Block 0+
( 1 2 5 4) // Block 1+
(24 25 2 1) // Block 2+
( 7 8 25 24) // Block 3+
(10 11 8 7) // Block 4+
( 9 10 7 6) // Block 5+
(15 16 13 12) // Block 0-
(16 17 14 13) // Block 1-
(13 14 27 26) // Block 2-
(26 27 20 19) // Block 3-
(19 20 23 22) // Block 4-
(18 19 22 21) // Block 5-
);
}
);

mergePatchPairs
(
);

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


All times are GMT -4. The time now is 08:54.