[alimea]

Hi
I made some changes on "laplacian solver" and I wanted to solve a simple problem. but I saw this error:

stream failure while reading header on line 55 of file "/home/drnili/Desktop/laplacianFoam/example/laplacianFoam/test/0/sai" for essential objectsai

Of course the file "sai" in "0" folder doesn't have 55 lines!!
Could you please help me?

[floquation]

If you open any of OpenFoam's case files (e.g. constant/transportProperties, but any file will do), you'll find a header that looks somewhat like:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

As you can see, there is information in this header. OpenFoam needs that information to read the file.

Therefore, you should include the above header (but then adapted for your file, of course). Also mind the OpenFoam version, so do not simply copy my header above, but rather take a header from one of your other case files.

[alimea]

Originally Posted by floquation

If you open any of OpenFoam's case files (e.g. constant/transportProperties, but any file will do), you'll find a header that looks somewhat like:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

As you can see, there is information in this header. OpenFoam needs that information to read the file.

Therefore, you should include the above header (but then adapted for your file, of course). Also mind the OpenFoam version, so do not simply copy my header above, but rather take a header from one of your other case files.

Dear Kevin
Thanks for your answer.
But I copy a similar header (e.g. the header of T in "0" folder) for the "sai" parameter. so what is the problem?

[alexeym]

Hi,

Surely game called "Guess what I have done wrong!" is fascinating and people here like to play it VERY much. Yet, it would be MUCH easier if:

1. You tell a little bit more about `some changes on "laplacian solver"`.
2. Post your sai file.

[alimea]

Originally Posted by alexeym

Hi,

Surely game called "Guess what I have done wrong!" is fascinating and people here like to play it VERY much. Yet, it would be MUCH easier if:

1. You tell a little bit more about `some changes on "laplacian solver"`.
2. Post your sai file.

I sent it to you on other question.
Regards

[alexeym]

Hi,

In other question you have posted your modifications. This question is about problem with sai file. And there is no sai file in this thread.

[alimea]

Originally Posted by alexeym

Hi,

In other question you have posted your modifications. This question is about problem with sai file. And there is no sai file in this thread.

This is the file of mylaplacian.C :

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------

\*---------------------------------------------------------------------------*/

#include "fvCFD.H"
#include "simpleControl.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
#include "setRootCase.H"

#include "createTime.H"
#include "createMesh.H"
#include "createFields.H"

simpleControl simple(mesh);

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Info<< "\nCalculating temperature distribution\n" << endl;

while (simple.loop())
{
Info<< "Time = " << runTime.timeName() << nl << endl;

while (simple.correctNonOrthogonal())
{
solve(fvm::laplacian(sai) == omega);

Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
}

Info<< "End\n" << endl;

return 0;
}


// ************************************************** *********************** //


This is the "createFields.H" :



Info<< "Reading field sai\n" << endl;

volScalarField sai
(
IOobject
(
"sai",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);


Info<< "Reading field omega\n" << endl;

volScalarField omega
(
IOobject
(
"omega",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

Info<< "Reading field gradU\n" << endl;

volTensorField gradU
(
IOobject
(
"gradU",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

forAll(omega, counter)
{
omega[counter] = gradU[counter].component(tensor::XY) - gradU[counter].component(tensor::YX);
}




Info<< "Reading transportProperties\n" << endl;

IOdictionary transportProperties
(
IOobject
(
"transportProperties",
runTime.constant(),
mesh,
IOobject::MUST_READ_IF_MODIFIED,
IOobject::NO_WRITE
)
);


Info<< "Reading diffusivity DT\n" << endl;

dimensionedScalar DT
(
transportProperties.lookup("DT")
);

[alimea]

Originally Posted by alexeym

Hi,

In other question you have posted your modifications. This question is about problem with sai file. And there is no sai file in this thread.

The last problem is solved.
Thank you.

I got the "gradU" file from the "100" folder of other solution (it means the results of that solution at 100 second) and paste it in the "0" folder of present problem as an initial condition. The problem was run but no result folder is written in the folder of problem!!
why?
Regards

[alexeym]

Cause you do not call runTime.write(). And it seems, you are playing another game called "Randomly put curly brackets and see if your code can be compiled".

[alimea]

Originally Posted by alexeym

Cause you do not call runTime.write(). And it seems, you are playing another game called "Randomly put curly brackets and see if your code can be compiled".

yes. Thank you