|
[Sponsors] |
July 6, 2005, 11:14 |
Hi,
I am trying to add an equ
|
#1 |
Member
Radu Mustata
Join Date: Mar 2009
Location: Zaragoza, Spain
Posts: 99
Rep Power: 17 |
Hi,
I am trying to add an equation to oodles and then run it for a "cold jet" simulation. The equation would be for a non-reacting scalar and itīs done in a similar way as "ftEqn.H" of XiFoam. I will breafly describe what Iīve done ("mixf", it`s the name of my variable): firstly: mixfEqn.H reads: tmp<fv::convectionscheme<scalar> > mvConvection ( fv::convectionScheme<scalar>::New ( mesh, phi, mesh.divScheme("div(phi,mixf_2)") ) ); solve ( fvm::ddt(mixf) + mvConvection->fvmDiv(phi, mixf) - fvm::laplacian(sgsModel->nuEff(), mixf) ); I added in "createAverages.H" volScalarField* mixfMeanPtr; if (nAveragingStepsFile.good()) { ... Info<< "Creating field mixfMean\n" << endl; mixfMeanPtr = new volScalarField ( IOobject ( "mixfMean", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); } else { ... Info<< "Creating field mixfMean\n" << endl; mixfMeanPtr = new volScalarField ( IOobject ( "mixfMean", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), mixf ); } .... volScalarField& mixfMean = *mixfMeanPtr; Then in "calculateAverages.H" I added: mixfMean = nm1Coeff*mixfMean + nCoeff*mixf; and in "createFields.H" Info<< "Reading field mixf\n" << endl; volScalarField mixf ( IOobject ( "mixf", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); now, at runtime, in the 0 directory I added the "mixf" file with all the definitions of patch values and all that I run the case in parallel using LAM and it seems to run fine, for as long as I left it running Co did not go over 0.5. The trouble that I have is that at finalising the program I get a message like: --> FOAM FATAL ERROR : NO_READ specified for read-constructor of object Umean of class IOobject Function: regIOobject::readStream(const word&) in file: db/regIOobject/regIOobjectRead.C at line: 53. FOAM parallel run aborting ...and this Umean can be something else like mixf_0, or pmean or mixfmean ...with the same error. Then, of course (I guess) when I want to reconstructPar the simulation I get (e.g. in the aforementioned error case): Time = 0.0005 Reconstructing FV fields Reconstructing volScalarFields mixf_0 nuSgs epsilonMean k mixfMean --> FOAM FATAL IO ERROR : cannot open file file: /mnt/store1/radu/OpenFOAM/radu-1.1/run/oodles2/coldjet/processor9/0.0005/mixfMea n at line 0. Function: regIOobject::readStream(const word&) in file: db/regIOobject/regIOobjectRead.C at line: 68. FOAM exiting Am I missing something? Help...please! Cheers, Radu |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Coefficienct in the oodles | loneboard | OpenFOAM Running, Solving & CFD | 4 | March 12, 2009 21:18 |
A question on DES oodles | ivan_cozza | OpenFOAM Running, Solving & CFD | 2 | November 17, 2008 09:33 |
Oodles vs turbFoam | rolando | OpenFOAM Running, Solving & CFD | 9 | June 4, 2007 05:42 |
Parallel oodles | sriharsha | OpenFOAM Running, Solving & CFD | 0 | December 27, 2006 17:25 |
Parallel oodles | sriharsha | OpenFOAM Running, Solving & CFD | 0 | December 27, 2006 17:21 |