|
[Sponsors] |
March 15, 2013, 05:45 |
C_VOF(c,pt) value increases every time step
|
#1 |
Member
Indrajit
Join Date: Nov 2012
Posts: 33
Rep Power: 14 |
Hi everyone,
I am simulating a two phase flow case in which a waterfront advances and I am using an UDF to track the interface. While I detect the advancing front and display the x-value of the centroid and also the VOF value in the cell by: C_CENTROID(xc,c,t) and C_VOF(c,pt) resply,where ,pt= THREAD_SUB_THREAD(t,1), is the phase thread and t is the mixture thread. I have found that though the centroid does not change in two time steps the VOF value in the cell increases. Does the VOF value get added instead of being calculated afresh every time step??? Please help. Regards, Indrajit |
|
March 18, 2013, 03:41 |
|
#2 |
Senior Member
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 23 |
Shouldn't 'pt' in C_VOF(c,pt) have an index e.g. pt(1) for the phase in question?
pt[i], an element in the array, is a pointer to the corresponding phase-level thread for the ith phase, where i is the phase_domain_index. Source: UDF Manual. |
|
March 18, 2013, 05:27 |
C_VOF increases
|
#3 | |
Member
Indrajit
Join Date: Nov 2012
Posts: 33
Rep Power: 14 |
Quote:
Vasava, Thank you for your attention. However my THREAD_SUB_THREAD(t,1) specifies the required index and hence pt is not array of threads in this case but just thread pointer to the phase thread.So I guess that should not be a problem. Any thing else that can affect the VOF value to increase??? Regards, Indrajit |
||
March 18, 2013, 06:16 |
|
#4 |
Senior Member
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 23 |
I can't think of anything right now. Does this happen in between the simulation or from the beginning.
|
|
March 18, 2013, 09:08 |
|
#5 | |
Member
Indrajit
Join Date: Nov 2012
Posts: 33
Rep Power: 14 |
Quote:
When I start the simulation it shows some value and then every iteration the value keeps on increasing!!! I have not checked how long the value keeps on increasing but increase in the value at the same place with the interface not moving seems odd!! |
||
March 19, 2013, 03:33 |
|
#6 |
Senior Member
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 23 |
Did you use 'patch' after initialization?
If not then I suspect that this behavior is the effect of the initial condition or the initial state of the flow. You can run your simulation for bit longer and observe C_VOF. If it does not stabilizes to a reasonable range then there may be something wrong with the setup as well. |
|
March 19, 2013, 04:40 |
|
#7 | |
Member
Indrajit
Join Date: Nov 2012
Posts: 33
Rep Power: 14 |
Quote:
Yes I had used patch after initialization. Perhaps this behavior is an effect of the initial flow. I will see if I can observe stable values along some later time steps. Regards, Indrajit |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
High Courant Number @ icoFoam | Artex85 | OpenFOAM Running, Solving & CFD | 11 | February 16, 2017 14:40 |
Superlinear speedup in OpenFOAM 13 | msrinath80 | OpenFOAM Running, Solving & CFD | 18 | March 3, 2015 06:36 |
How to write k and epsilon before the abnormal end | xiuying | OpenFOAM Running, Solving & CFD | 8 | August 27, 2013 16:33 |
plot over time | fferroni | OpenFOAM Post-Processing | 7 | June 8, 2012 08:56 |
Problem with FloatingObject | Leech | OpenFOAM Running, Solving & CFD | 10 | March 29, 2012 16:24 |