CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

works in serial, but not parallel - GGDH in rhoSimpleFoam

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

LinkBack Thread Tools Search this Thread Display Modes
Old   April 30, 2017, 21:22
Default works in serial, but not parallel - GGDH in rhoSimpleFoam
Join Date: Aug 2012
Posts: 47
Rep Power: 12
JackW is on a distinguished road
Hi all,

I have a puzzling issue implementing a solver that incorporates Daly & Harlow's GGDH model in the energy Equation. I have updated the EEqn.H of rhoSimpleFoam to read:

volScalarField fluxTimeScale (
        0.3 * turbulence->mut() / (0.09 * rho * turbulence->k())

    volSymmTensorField alphaEff(

    fvScalarMatrix EEqn
        fvm::div(phi, he)
      + (
   == "e"
          ? fvc::div(phi, volScalarField("Ekp", 0.5*magSqr(U) + p/rho))
          : fvc::div(phi, volScalarField("K", 0.5*magSqr(U)))
        //   - fvm::laplacian(turbulence->alphaEff(), he) // linear
      - fvm::laplacian(alphaEff, he) // non-linear
      - fvm::laplacian(turbulence->alpha(), he)

        fvOptions(rho, he)
I am using the kOmegaSST model, so k(), R() and mut() all return well defined values.

In serial everything runs, but in parallel I get:

[1] Maximum number of iterations exceeded
[1]     From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
[1]     in file /home/jack/OpenFOAM/OpenFOAM-2.4.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.
FOAM parallel run aborting
I am using sensibleEnthalpy, perfectGas and const for transport. As I run the simulation I watch the h residuals rise and I get a crash after ~5-15 iterations. This is obviously the cause as h is proportional to T. Note I have also tried solving for e.

If I look at the flow in paraview, I see small regions (1 or 2 cells) of |T|>3e+8 (should be 297-330).

Note, this behaviour is reasonably fvScheme/fvSolution dependent and if I change the number of procs and decomposition. Also, I have tried from a converged solution from using the standard heat flux model and an initial condition.

Note that I have two timeVaryingMappedFixedValue boundary conditions.

if anyone can help - thanks! I am guessing it is something to do with the processor boundaries and using the turbulence member functions, but I do not know.


JackW is offline   Reply With Quote

Old   November 15, 2019, 06:49
Join Date: Dec 2012
Posts: 32
Rep Power: 11
mechkween is on a distinguished road
Hi Jack,

Did you manage to find a solution to this problem? I've also implemented a GGDH in the EEqn.H and my simulation crashes stating a floating point exception.

mechkween is offline   Reply With Quote


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
Explicitly filtered LES saeedi Main CFD Forum 16 October 14, 2015 11:58
Poisson eq w setReference works serial diverges in parallel tehache OpenFOAM Running, Solving & CFD 5 August 29, 2012 09:41
FLUENT in Windows Serial and Linux Parallel sham83 FLUENT 0 May 31, 2011 23:38
Differences solutions in Parallel and Serial wizard1301 OpenFOAM Running, Solving & CFD 2 January 30, 2009 05:38
Serial run OK parallel one fails r2d2 OpenFOAM Running, Solving & CFD 2 November 16, 2005 12:44

All times are GMT -4. The time now is 19:37.