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

interFoam error with imported Gambit mesh

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 19, 2013, 00:13
Default interFoam error with imported Gambit mesh
  #1
Member
 
Yao Lu
Join Date: May 2013
Posts: 33
Rep Power: 12
shuoxue is on a distinguished road
I am solving the damBreak tutorial in user guide with imported Gambit mesh.
Boundary conditions in boundary file have been rewrited.

This is .msh file generated in Gambit.
https://skydrive.live.com/?cid=a9e4f...4756B9AE%21105

setFields is ok.
1.jpg

Error occurs when I solve the case with interFoam.
Code:
luyao@luyao:~/OpenFOAM/luyao-2.2.0/run/case/multiphase-interFoam/1/damBreak2$ interFoam
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.2.0-5be49240882f
Exec   : interFoam
Date   : Jun 19 2013
Time   : 10:26:21
Host   : "luyao"
PID    : 4020
Case   : /home/luyao/OpenFOAM/luyao-2.2.0/run/case/multiphase-interFoam/1/damBreak2
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar

Reading g
Calculating field g.h

No finite volume options present

time step continuity errors : sum local = 0, global = 0, cumulative = 0
DICPCG:  Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
deltaT = 0.00119048
Time = 0.00119048

MULES: Solving for alpha1
Phase-1 volume fraction = 0.130194  Min(alpha1) = 0  Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.130194  Min(alpha1) = 0  Max(alpha1) = 1
DICPCG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.0365699, No Iterations 3
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Uninterpreted: 
#3   at tensorField.C:0
#4  
 in "/opt/openfoam220/platforms/linuxGccDPOpt/bin/interFoam"
#5  
 in "/opt/openfoam220/platforms/linuxGccDPOpt/bin/interFoam"
#6  
 in "/opt/openfoam220/platforms/linuxGccDPOpt/bin/interFoam"
#7  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#8  
 in "/opt/openfoam220/platforms/linuxGccDPOpt/bin/interFoam"
Thanks in advance.
shuoxue is offline   Reply With Quote

Old   June 19, 2013, 08:08
Default
  #2
Member
 
Nicklas Linder
Join Date: Jul 2012
Location: Germany
Posts: 35
Rep Power: 13
nlinder is on a distinguished road
Hi!

I have similar problems, but only with a modified interFoam solver so far. Interesting that it also happens with the standard solver. What are your boundary conditions?

Greetings
nlinder is offline   Reply With Quote

Old   June 19, 2013, 08:28
Default
  #3
Member
 
Nicklas Linder
Join Date: Jul 2012
Location: Germany
Posts: 35
Rep Power: 13
nlinder is on a distinguished road
Update:

Can you confirm that commenting out the first two lines here (pEqn.H):
Code:
        if (pimple.finalNonOrthogonalIter())
        {
            //phi = phiHbyA - p_rghEqn.flux();

            //U = HbyA + rAU*fvc::reconstruct((phig - p_rghEqn.flux())/rAUf);
            U.correctBoundaryConditions();
            fvOptions.correct(U);
        }
"solves" the problem (at least the error message is gone - solution obviously not correct)?

The reason at this point is, that p_rghEqn.flux() returns a "-nan" field. Your first written timestep should be "-nan" everywhere in p_rgh, correct?

BTW: Serial or parallel?

greetings

Last edited by nlinder; June 19, 2013 at 09:49.
nlinder is offline   Reply With Quote

Old   June 19, 2013, 10:30
Default
  #4
Member
 
Nicklas Linder
Join Date: Jul 2012
Location: Germany
Posts: 35
Rep Power: 13
nlinder is on a distinguished road
Sorry for answering myself all the time but:
It seems as if
Code:
- ghf*fvc::snGrad(rho)
in "surfaceScalarField phig" produces an "-nan" on a boundary. Not in all cases, therefore I do not yet know where it comes from. Can you maybe confirm that in your case? Just comment it out and run again.

Greetings
Nicklas
nlinder is offline   Reply With Quote

Old   June 21, 2013, 00:42
Default
  #5
Member
 
Yao Lu
Join Date: May 2013
Posts: 33
Rep Power: 12
shuoxue is on a distinguished road
Quote:
Originally Posted by nlinder View Post
Sorry for answering myself all the time but:
It seems as if
Code:
- ghf*fvc::snGrad(rho)
in "surfaceScalarField phig" produces an "-nan" on a boundary. Not in all cases, therefore I do not yet know where it comes from. Can you maybe confirm that in your case? Just comment it out and run again.

Greetings
Nicklas
Hi, nlinder!
I am afraid that we got different problems although we received similar error.
Mine has been solved, see here.

Which line do you want me to comment?
Both
Code:
phi = phiHbyA - p_rghEqn.flux();
U = HbyA + rAU*fvc::reconstruct((phig - p_rghEqn.flux())/rAUf);
and
Code:
- ghf*fvc::snGrad(rho)
shuoxue is offline   Reply With Quote

Old   June 21, 2013, 04:48
Default
  #6
Member
 
Nicklas Linder
Join Date: Jul 2012
Location: Germany
Posts: 35
Rep Power: 13
nlinder is on a distinguished road
Hi!

Allright, thanks for the information. In my case it was enough to comment out
Code:
- ghf*fvc::snGrad(rho)
I now specified a boundary condition for rho and in my testCase I had some issues with the pressure and Velocity BC. I think I somehow "solved" that problem currently.
Thanks for discussing

Nicklas
nlinder 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
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
Question on InterFoam moving mesh capabilities ziv OpenFOAM Running, Solving & CFD 0 April 23, 2008 10:11
Convergence moving mesh lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 15:09
Successive Mesh in Gambit ashish FLUENT 0 April 28, 2006 01:42
How to mesh a circle by Gambit? Zhengcai Ye FLUENT 5 March 24, 2006 02:04


All times are GMT -4. The time now is 06:59.