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

Oodles and more

Register Blogs Community New Posts Updated Threads Search

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   July 6, 2005, 11:14
Default Hi, I am trying to add an equ
  #1
Member
 
Radu Mustata
Join Date: Mar 2009
Location: Zaragoza, Spain
Posts: 99
Rep Power: 17
r2d2 is on a distinguished road
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
r2d2 is offline   Reply With Quote

 


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
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


All times are GMT -4. The time now is 07:45.