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

interFoam: No finite volume options present

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 10, 2014, 15:59
Default interFoam: No finite volume options present
  #1
New Member
 
CFDGuy
Join Date: Nov 2013
Posts: 5
Rep Power: 3
mehdinekooee is on a distinguished road
Hi Foamers,

When I run interFoam on a cluster I got the following output:
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.0592593, global = -0.0592593, cumulative = -0.0592593
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 7.78695e-11, No Iterations 221
time step continuity errors : sum local = 5.71578e-12, global = 4.65865e-14, cumulative = -0.0592593
Courant Number mean: 12.7285 max: 41.8982

Starting time loop
-----
"No finite volume options present" what does this phrase mean?
Does it have any effect on cpu time ?

Thanks
mehdinekooee is offline   Reply With Quote

Old   January 12, 2014, 08:37
Default
  #2
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 919
Rep Power: 17
akidess will become famous soon enough
I'm guessing you are running a case without an fvOptions file. Shouldn't be a problem, just ignore the warning.
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Help define the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oam-technology
akidess is offline   Reply With Quote

Old   July 14, 2014, 03:57
Default
  #3
Member
 
Vignesh
Join Date: Oct 2012
Location: Darmstadt, Germany
Posts: 33
Rep Power: 4
vigneshTG is on a distinguished road
Dear Friends !!

I recently modified the surface tension force implementation in the interfoam solver and successfully compiled it. But when i run a case using this solver i get error as shown in the snapshot.

The modification i made was F_{\sigma}=\sigma* \nabla^{2}\alpha* \hat{n} instead of F_{\sigma}=\sigma*\kappa* \nabla \alpha

I made the changes in PEqn.H file in interfoam directory and the modified code is shown below (modifications are written in bold)
Code:
{
    
    /* New one to calculate laplacian of alpha*/
    volScalarField lapalpha(fvc::laplacian(alpha1));
     /* New one to calculate laplacian of alpha*/

    volVectorField HbyA("HbyA", U);
    HbyA = rAU*UEqn.H();

    surfaceScalarField phiHbyA
    (
        "phiHbyA",
        (fvc::interpolate(HbyA) & mesh.Sf())
      + fvc::ddtPhiCorr(rAU, rho, U, phi)
    );

    adjustPhi(phiHbyA, U, p_rgh);
    phi = phiHbyA;

    surfaceScalarField phig
    (
        (
 
                fvc::interpolate(lapalpha)*interface.sigma()*interface.nHatf()
           - ghf*fvc::snGrad(rho)
        )*rAUf*mesh.magSf()
    );

    phiHbyA += phig;

    while (pimple.correctNonOrthogonal())
    {
        fvScalarMatrix p_rghEqn
        (
            fvm::laplacian(rAUf, p_rgh) == fvc::div(phiHbyA)
        );

        p_rghEqn.setReference(pRefCell, getRefCellValue(p_rgh, pRefCell));

        p_rghEqn.solve(mesh.solver(p_rgh.select(pimple.finalInnerIter())));

        if (pimple.finalNonOrthogonalIter())
        {
            phi = phiHbyA - p_rghEqn.flux();

            U = HbyA + rAU*fvc::reconstruct((phig - p_rghEqn.flux())/rAUf);
            U.correctBoundaryConditions();
            fvOptions.correct(U);
        }
    }

    #include "continuityErrs.H"

    p == p_rgh + rho*gh;

    if (p_rgh.needReference())
    {
        p += dimensionedScalar
        (
            "p",
            p.dimensions(),
            pRefValue - getRefCellValue(p, pRefCell)
        );
        p_rgh = p - rho*gh;
    }
}
Can someone tell me what is this error and how to solve this ?



Attached Images
File Type: jpg Screenshot from 2014-07-14 09:02:55.jpg (83.0 KB, 58 views)
__________________
Thanks and Regards

Vignesh TG
vigneshTG is offline   Reply With Quote

Old   July 15, 2014, 04:57
Default
  #4
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 185
Rep Power: 7
vonboett is on a distinguished road
I bet you devide by zero in some cells, yust give it a try with lapalpha(alpha1) with alpha1 > 0 in all cells to try it out.
vonboett is offline   Reply With Quote

Old   July 15, 2014, 11:05
Default
  #5
Member
 
Vignesh
Join Date: Oct 2012
Location: Darmstadt, Germany
Posts: 33
Rep Power: 4
vigneshTG is on a distinguished road
Thank you Albrecht !! I found the error and corrected it
__________________
Thanks and Regards

Vignesh TG
vigneshTG 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
dynamic Mesh is faster than MRF???? sharonyue OpenFOAM Running, Solving & CFD 14 August 26, 2013 07:47
BlockMesh FOAM warning gaottino OpenFOAM Native Meshers: blockMesh 7 July 19, 2010 14:11
About finite volume methods for N-S Lionel BRS Main CFD Forum 6 July 16, 2006 03:15
Finite Element vs. Finite Volume Codes Salman Main CFD Forum 9 January 6, 2006 23:59
control volume based finite element araz Main CFD Forum 1 November 15, 2003 18:50


All times are GMT -4. The time now is 03:23.