[dokeun]

Originally Posted by Hisham

I got it. I post here for future readers.

1. In the blockMeshDict, define the cellZone names for required blocks:
e.g. (from the damBreak example)

blocks
(
hex (0 1 5 4 12 13 17 16) cellZoneName (23 8 1) simpleGrading (1 1 1)
hex (2 3 7 6 14 15 19 18) (19 8 1) simpleGrading (1 1 1)
hex (4 5 9 8 16 17 21 20) (23 42 1) simpleGrading (1 1 1)
hex (5 6 10 9 17 18 22 21) cellZoneName (4 42 1) simpleGrading (1 1 1)
hex (6 7 11 10 18 19 23 22) (19 42 1) simpleGrading (1 1 1)
);


2. In the setFieldsDict,

defaultFieldValues
(
volScalarFieldValue alpha1 0
);

regions
(
zoneToCell
{
name "cellZoneName";
fieldValues
(
volScalarFieldValue alpha1 1
);
}
);

Dear Hisham

Thank you for your post and comment for followers like me.

I'd also like to use non-box region initialization for the flow on cylinderial wall with multiphaseInterFoam.

The difference is that I import gmsh for mesh using gmshToFoam and there was no problem but I have a problem during setFields execution.

the code I have is as below.

deaultFieldsValues
(
    volScalarFieldValue alpha.water 1
    volScalarFieldValue alpha.air 0
    volVectorFieldValue U (0 0 0)
)

regions
(
  patchToFace
  {
   name inlet_water;
   fieldValues
   (
    ....
   );
  }
  patchToFace{...}
  patchToFace{...}
  
  zoneToCell
  {  
    name internal_water;
    fieldValues
    (
      volScalarFieldValue alpha.water 1
      volScalarFieldValue alpha.air 0
      volVectorFieldValue U (0 -0.5 0)
    );
  }
)

When I run setFields then I have a error message

-->  FOAM FATAL IO ERROR:

keyword name is undefined in dictionary ".zoneToCell"

file: .zoneToCell

    From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const in file db/dictionary/dictionary.C at line 566.

FOAM exiting

But when I check the Valid type by typing errata, it recommand the zoneToCell as one of the Valid type.

I guess zoneToCell is to be applied in setFields but I didn't get any idea from tutorials which use zoneToCell. They use zoneToCell in different utilit
y other than setFields.

I solved the problem with other region option, cylinderAnnulusToCell but it is little bit annoying because I have to know the coordinates.

Is there any fault that I have tried? And how can I fix it?

Thank you in advance.

[Owais Shabbir]

Hi Dokeun,


Have you found any solution to it? I am in a similar problem right now.

Thanks
OS

[dokeun]

Originally Posted by Owais Shabbir

Hi Dokeun,


Have you found any solution to it? I am in a similar problem right now.

Thanks
OS

Hi
You might need to check the coincidence of the name for ‘set’ in refineMeshDict file and the name in the topoSetDict.
Good luck.

[RANSES]

Hello to everybody. I'm profiting of this open thread for asking a simple question about setFields. By default, setFields is acting on the 0 folder. Do you know if you can use it also on other time-steps (i.e. time folders)?
If yes, how?
If no, do you know if there is an alternative for doing this operation?

thank you in advance.

[dokeun]

Originally Posted by RANSES

Hello to everybody. I'm profiting of this open thread for asking a simple question about setFields. By default, setFields is acting on the 0 folder. Do you know if you can use it also on other time-steps (i.e. time folders)?
If yes, how?
If no, do you know if there is an alternative for doing this operation?

thank you in advance.

Hello, Riccardo

When I need to run OpenFOAM case for non-zero time step, there are previous data for it. No need to use setField. Just one thin to do is changing the initial time from what you want to start in ControlDict.

If you want to start a case from other data set, then you need to use mapfields.

Good luck.

[RANSES]

Originally Posted by dokeun

Hello, Riccardo

When I need to run OpenFOAM case for non-zero time step, there are previous data for it. No need to use setField. Just one thin to do is changing the initial time from what you want to start in ControlDict.

If you want to start a case from other data set, then you need to use mapfields.

Good luck.

Hi Dokeun,

and thank you very much for your reply. The problem I'm facing is a dispersion study. I have to simulate a pollutant dispersion phenomenon starting from a converged state of the flow-field....that's why I want/need to apply setField on the concentration field in a time-step later than 0. I think I resolved this issue by copying the field I have modified with setField in the 0 folder in the one of the converged flow-field (kind of 2s). I was wondering if there was a "more elegant" way for doing this operation.

[gschaider]

Originally Posted by RANSES

Hi Dokeun,

and thank you very much for your reply. The problem I'm facing is a dispersion study. I have to simulate a pollutant dispersion phenomenon starting from a converged state of the flow-field....that's why I want/need to apply setField on the concentration field in a time-step later than 0. I think I resolved this issue by copying the field I have modified with setField in the 0 folder in the one of the converged flow-field (kind of 2s). I was wondering if there was a "more elegant" way for doing this operation.

I'd say "Use funkySetFields. It has a -time-option" But I'm extremely biased

[RANSES]

Originally Posted by gschaider

I'd say "Use funkySetFields. It has a -time-option" But I'm extremely biased

Thank you Bernhard, I will give a look to it.
Also, I'm having problem at defining setFields for a volVectorField (I've bee using it just for volScalarField). Do you know if setFields can be used also for vector quantities?

[dokeun]

Originally Posted by RANSES

Hi Dokeun,

and thank you very much for your reply. The problem I'm facing is a dispersion study. I have to simulate a pollutant dispersion phenomenon starting from a converged state of the flow-field....that's why I want/need to apply setField on the concentration field in a time-step later than 0. I think I resolved this issue by copying the field I have modified with setField in the 0 folder in the one of the converged flow-field (kind of 2s). I was wondering if there was a "more elegant" way for doing this operation.

Sorry for late reply.
I’m not so sure about my opinion because I’m newbie, too.
Anyway, as I think....
You might have a solution data without particles (and BC without particle definition) already from previous calculation for improving convergence.
And the field data folder for 2s are to be stored about the same calculation domain.

Then make a new case copying previous folders and files(system, control, 2.00, mesh file)
Make a mapfieldsPar file from examples or tutorials in control folder and change initial time as 2 sec in controlDict file.
Modify SetfieldPar so as to include boundary condition for particle with previous bc.
Then, just run SetfieldsPar, mapFieldsPar, Solver.

I hope you can crack the problem.