[LUIS1717]

Hi everyone,

I'm trying to program a solver and its need a volscalarfield that initially must be zero everywhere in the mesh and I don't want to read it from any initial condition or to write it down. Anyone know how to do it?

thanks in advance Luis

[alexeym]

Hi,

Just an example:

$ sol
$ grep -r volScalarField *
...lots of output...
stressAnalysis/solidDisplacementFoam/readThermalProperties.H
...
$ less stressAnalysis/solidDisplacementFoam/readThermalProperties.H

Go to line 16.

Another:

$ src
$ cd turbulenceModels/incompressible/turbulenceModel/laminar
$ less laminar.C

Go to line 90 (this one maybe a little bit too complicated, since involves certain things in addition to field creation).

[LUIS1717]

Thanks alexeym,

I was wondering if there were some way to avoid to create the IOobject, or at least dont have to write it down.

[rapierrz]

Hi Luis,

For example for vol scalar field alpha

volScalarField Temp
(
IOobject
(
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
dimensionedScalar("Temp",dimensionSet(0,0,0,0,0,0, 0),0)
);

[alexeym]

@LUIS1717

If you take a look at GeometricField class documentation (http://foam.sourceforge.net/docs/cpp/a00912.html, volScalarField is just a typedef GeometricField<scalar, fvPatchField, volMesh>), it can be created either as a copy of another GeometricField, or on the base of IOobject.

So in brief, no, you can not avoid IOobject constructor.

Surely if you really-really-really hate IOobject constructor, you can create wrapper taking field name, dimensions, and mesh as parameters and returning autoPtr.

[ziad]

Originally Posted by LUIS1717

Hi everyone,

I'm trying to program a solver and its need a volscalarfield that initially must be zero everywhere in the mesh and I don't want to read it from any initial condition or to write it down. Anyone know how to do it?

thanks in advance Luis

Here's a REALLY ugly hack, but it's simple, it works and you won't need IOobject.

volScalarField toto ("use any combination of variables with correct dimensions");
toto *= scalar(0);