It's interesting to see how the numerical diffusion manifests itself only in the direction of local convection, resulting first in an ellipse, and a circle after the second translation. Neat!
I am trying to simulate colliding droplets with multiphase VOF solver( interFoam) in OpenFOAM so, the problem I'm having is to give velocities of same magnitude but opposite in direction to droplets, using internal fields won't help same as with setFields which accepts only scalar fields so, any help in this would be highly appreciated
Hi Jozef thanks a lot for the informative tutorial. The velocity field you had mentioned is constant. How and we simulate in an accelerating field (for example U = t, where t in time)? If possible kindly made a video.
Create time Create mesh for time = 0 Reading setFieldsDict Setting field default values Setting internal values of volScalarField T Setting field region values --> FOAM FATAL IO ERROR: (openfoam-2012) Entry 'radius' has 6 excess tokens in stream 7(1 fieldValues ( volScalarFieldValue T 1 )) file: sphereToCell at line 29. FOAM exiting @József Nagy, any solution for this?
Thank you very much Jozsef Nagy sir. I tried by my self without looking this video. And I end up making 3 different cases . Lol. Then I saw this video and realised my mistake then I corrected it.... Thank you very much sir.
I've got it without explanation, it was not mentioned in previous tutorials, but one could infer the logic, how the simulation is run, and folders for each timestep are created
Hi professor, so even for" sphereToCell", it will generate a cylinder of radius 1 for this case? .....it seems using cylinderToCell is more intuitive. Based on the word "sphere", it seems spheretocell might generate a ball shape region at 0 0 0, and the circle shown on the front and back "empty plane" might not have the radius of 1.
Hi, How to know the required physical properties (eg. thermophysical properties ,transport properties..etc in the constant directory) for each solver? For example, if I want to solve a case using the icoFoam solver, may be I have to copy a basecase from an icoFoam tutorial and check out the required properties from that?
When I click the button in paraview to see the next step, it jumps over a period of time of 1.3333s rather than 1s. This causes paraview never to show the folder after 2 seconds since it jumps from 1.333 (approximately 1) to 2.6666 (approximately 3). The folders that are created in my file however only contain integers and I can also put these integers manually. Is there a reason for this strange behavior of paraview?
I am never able to see the initial conditions in paraview. Only after I run my simulation I am. I am talking about 12.09. This happened during the grid convergence tutorial as well. But if I use paraFoam instead of paraview & then I can see the initial values? why so?
@@OpenFOAMJozsefNagy I have got the same, I unchecked this box and refreshed, applied. I cannot see the initial value on the mesh which appears completely blue, right from the first second photogram the circle appears correctly. It is strange because paraFoam opens another version of paraview and there it is visible also the first photogram. I think this is a bug.
I fixed the problem, it was my fault. Just for knowledge, my mistake was this: in the 0 folder I kept the file names with the extension: T.org and U.org. I repeated changing the name in T, U (without .org) and I was able to see also the time 0 configuration. I am sorry for this.
how do we do it in one go ? l tried giving the if condition using runTime in the U file. But, I guess that was the wrong way to do. gave me so many errors which I failed to solve.
Hiii. In your youtube video "How to run your first openfoam case yourself" though I deleted phi, still the scalar transport is going out of the domain which is suppose to go in positive y direction. Please help me
@@OpenFOAMJozsefNagy in the first run, I ran for 4 seconds. After that I deleted the phi in the 4th second time step. And I changed the velocity in 0 folder to (0 1 0).in the control dict I changed endtime to 8 seconds. But still it's going out of domain
why things come to wrong ? Reading setFieldsDict Setting field default values Setting internal values of volScalarField T Setting field region values Adding cells with centre within sphere, with centre = (0 0 0) and radius = 1 --> FOAM Warning : From function Foam::autoPtr setCellField::iNew::operator()(Foam::Istream&) const in file setFields.C at line 174 field type volScalarFieldvalue not currently supported --> FOAM Warning : From function Foam::autoPtr setCellField::iNew::operator()(Foam::Istream&) const in file setFields.C at line 174 field type T not currently supported --> FOAM FATAL IO ERROR: wrong token type - expected word, found on line 27 the label 1 file: .sphereToCell.fieldValues at line 27. From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::word&) in file primitives/strings/word/wordIO.C at line 74. FOAM exiting
Create time Create mesh for time = 0 SIMPLE: no convergence criteria found. Calculations will run for 4 steps. Reading field T Reading field U Reading transportProperties Reading diffusivity DT Reading/calculating face flux field phi No finite volume options present Calculating scalar transport Courant Number mean: 0.1 max: 0.1 Time = 0.01 --> FOAM FATAL IO ERROR: (openfoam-2012) Entry 'T' not found in dictionary "" From const Foam::dictionary& Foam::dictionary::subDict(const Foam::word&, Foam::keyType::option) const in file db/dictionary/dictionary.C at line 537. FOAM exiting I got this error here 3 times, and still, I am unable to solve it. @József Nagy, can you tell me what it is?
You forgot to take some of the dictionaries. Check 0/T and system/fvSchemes and especially system/fvSolution, whether you have the correct entries about T.
It's interesting to see how the numerical diffusion manifests itself only in the direction of local convection, resulting first in an ellipse, and a circle after the second translation. Neat!
13:23 to go down to last line of the file in nano editor, press ctrl+end key combination.
Your videos are very useful for the beiginners like me. I have finished the basic and now moving forward to intermediate. Thanks a lot sir
Thanks a lot :)
Moving to the intermediate tutorials set.
I am trying to simulate colliding droplets with multiphase VOF solver( interFoam) in OpenFOAM so, the problem I'm having is to give velocities of same magnitude but opposite in direction to droplets, using internal fields won't help same as with setFields which accepts only scalar fields so, any help in this would be highly appreciated
Thank you for great tutorial for beginners
Hi Jozef thanks a lot for the informative tutorial. The velocity field you had mentioned is constant. How and we simulate in an accelerating field (for example U = t, where t in time)? If possible kindly made a video.
U = t means U = (t 0 0)
Well there are some function objects, where you can add a code. There you could overwrite the velocity each time step. This could work.
Create time
Create mesh for time = 0
Reading setFieldsDict
Setting field default values
Setting internal values of volScalarField T
Setting field region values
--> FOAM FATAL IO ERROR: (openfoam-2012)
Entry 'radius' has 6 excess tokens in stream
7(1 fieldValues ( volScalarFieldValue T 1 ))
file: sphereToCell at line 29.
FOAM exiting
@József Nagy, any solution for this?
I thank you Jozsef, you are the best!
Thank you very much Jozsef Nagy sir. I tried by my self without looking this video. And I end up making 3 different cases . Lol.
Then I saw this video and realised my mistake then I corrected it.... Thank you very much sir.
Great!
start from latest time in controldict some beginner couldn't catch this.
Thank you for life saver point!
I've got it without explanation, it was not mentioned in previous tutorials, but one could infer the logic, how the simulation is run, and folders for each timestep are created
Hi professor, so even for" sphereToCell", it will generate a cylinder of radius 1 for this case? .....it seems using cylinderToCell is more intuitive. Based on the word "sphere", it seems spheretocell might generate a ball shape region at 0 0 0, and the circle shown on the front and back "empty plane" might not have the radius of 1.
Yes, because we only have one cell over the thickness and they are inside both the sphere and the cylinder.
Hi, How to know the required physical properties (eg. thermophysical properties ,transport properties..etc in the constant directory) for each solver? For example, if I want to solve a case using the icoFoam solver, may be I have to copy a basecase from an icoFoam tutorial and check out the required properties from that?
That is a very good idea, yes.
Thanks a lot for this great tutorial
When I click the button in paraview to see the next step, it jumps over a period of time of 1.3333s rather than 1s. This causes paraview never to show the folder after 2 seconds since it jumps from 1.333 (approximately 1) to 2.6666 (approximately 3). The folders that are created in my file however only contain integers and I can also put these integers manually. Is there a reason for this strange behavior of paraview?
Reload your case. This is a problem of paraview
Happens with me everytime. I hate it.
sorry for the question, but what does the boundery inletOutlet means? I read the user and I still didnt get it
the user guide*
If the flow is pointing outwards it is zeroGradient if the flow points inward the fix value is used which you provide as inletValue.
I can not see the temperature at all in paraview same for the velocity, in my 0 folder there is the Temperature file how can this be?
On the panel there is an option says ignore the zero time. Make sure it is not selected
I am never able to see the initial conditions in paraview. Only after I run my simulation I am. I am talking about 12.09. This happened during the grid convergence tutorial as well. But if I use paraFoam instead of paraview & then I can see the initial values? why so?
In recent versions of Paraview, thee is an annoying check box with "Skip zero time". Uncheck it.
@@OpenFOAMJozsefNagy I have got the same, I unchecked this box and refreshed, applied. I cannot see the initial value on the mesh which appears completely blue, right from the first second photogram the circle appears correctly. It is strange because paraFoam opens another version of paraview and there it is visible also the first photogram. I think this is a bug.
I fixed the problem, it was my fault. Just for knowledge, my mistake was this: in the 0 folder I kept the file names with the extension: T.org and U.org. I repeated changing the name in T, U (without .org) and I was able to see also the time 0 configuration. I am sorry for this.
@@giacomodelbianco I am glad it works.
how do we do it in one go ? l tried giving the if condition using runTime in the U file. But, I guess that was the wrong way to do. gave me so many errors which I failed to solve.
You do it in separate steps.
Hello, in my
OpenFoam folder I'm missing the setFieldsDict file in the tutorials /
basic / scalarTransportFoam / location, how do I solve it?
You can copy it from \tutorials\compressible\sonicFoam\laminar\shockTube\system
After finishing the basic, which video should I watch? intermediate?btw, your video is really useful!
I am glad. Yes, intermediate is a good idea.
Hiii.
In your youtube video "How to run your first openfoam case yourself" though I deleted phi, still the scalar transport is going out of the domain which is suppose to go in positive y direction. Please help me
Did you run the case from the last time step?
@@OpenFOAMJozsefNagy in the first run, I ran for 4 seconds. After that I deleted the phi in the 4th second time step. And I changed the velocity in 0 folder to (0 1 0).in the control dict I changed endtime to 8 seconds. But still it's going out of domain
@@mmaruthiprasanna Change it in folder 4
@@OpenFOAMJozsefNagy yes. It worked. Thankyou very much
why things come to wrong ?
Reading setFieldsDict
Setting field default values
Setting internal values of volScalarField T
Setting field region values
Adding cells with centre within sphere, with centre = (0 0 0) and radius = 1
--> FOAM Warning :
From function Foam::autoPtr setCellField::iNew::operator()(Foam::Istream&) const
in file setFields.C at line 174
field type volScalarFieldvalue not currently supported
--> FOAM Warning :
From function Foam::autoPtr setCellField::iNew::operator()(Foam::Istream&) const
in file setFields.C at line 174
field type T not currently supported
--> FOAM FATAL IO ERROR:
wrong token type - expected word, found on line 27 the label 1
file: .sphereToCell.fieldValues at line 27.
From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::word&)
in file primitives/strings/word/wordIO.C at line 74.
FOAM exiting
volScalarFieldvalue --> volScalarFieldValue
I got the same error @jozsef Nazy, sir please clarify what you said in this
Create time
Create mesh for time = 0
SIMPLE: no convergence criteria found. Calculations will run for 4 steps.
Reading field T
Reading field U
Reading transportProperties
Reading diffusivity DT
Reading/calculating face flux field phi
No finite volume options present
Calculating scalar transport
Courant Number mean: 0.1 max: 0.1
Time = 0.01
--> FOAM FATAL IO ERROR: (openfoam-2012)
Entry 'T' not found in dictionary ""
From const Foam::dictionary& Foam::dictionary::subDict(const Foam::word&, Foam::keyType::option) const
in file db/dictionary/dictionary.C at line 537.
FOAM exiting
I got this error here 3 times, and still, I am unable to solve it.
@József Nagy, can you tell me what it is?
You forgot to take some of the dictionaries. Check 0/T and system/fvSchemes and especially system/fvSolution, whether you have the correct entries about T.