CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

Continuity Residuals

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 29, 2013, 13:53
Default Continuity Residuals
  #1
Member
 
Shashank
Join Date: Apr 2011
Posts: 74
Rep Power: 14
shashank312 is on a distinguished road
I am solving a two-phase VOF model with phase change (both heat and mass transfer across the interface). Just wondering if the residual values of continuity could be greater than 1? I am observing residuals with values close to 1e+03 but they seem to be decreasing and eventually reaching to the convergence value (set at 1e-04).
shashank312 is offline   Reply With Quote

Old   May 29, 2013, 20:57
Default
  #2
Senior Member
 
Join Date: Aug 2011
Posts: 421
Blog Entries: 1
Rep Power: 21
blackmask will become famous soon enough
The residual is scaled by the maximum of first a few iterations (default 5). It the residual keep on growing after that then it could be larger than 1. It does not matter if it converges finally.
blackmask is offline   Reply With Quote

Old   May 30, 2013, 16:06
Default
  #3
Member
 
Shashank
Join Date: Apr 2011
Posts: 74
Rep Power: 14
shashank312 is on a distinguished road
Thanks again, Thomas. I kinda had an idea that it is converging but its taking a lot of iterations each time step. Since my evaporation rates and difference in densities of liquid and vapor are high at the liquid-vapor interface, my initial residuals might be high.

This is not related to this particular forum but do you know how to impose a constant temperature inside one of the phases (in a VOF model) throughout the simulation (every iteration of each time step)? I tried the following but it did not work.

DEFINE_ADJUST(somename,domain)
Thread *lt = pt[liquidphaseindex];

begin_c_loop_all(c,t)
{
if (C_VOF(c,lt) == 0.0)
{
C_T(c,t) = 329;
}
}
end_c_loop_all(c,t)

The temperature inside that phase is decreasing to a much lower value and after a few iterations reaching to zero Kelvin.
shashank312 is offline   Reply With Quote

Old   May 30, 2013, 20:49
Default
  #4
Senior Member
 
Join Date: Aug 2011
Posts: 421
Blog Entries: 1
Rep Power: 21
blackmask will become famous soon enough
One quick question, why
Code:
begin_c_loop_all(c,t)
instead of
Code:
begin_c_loop_all(c,lt)
Due to physical and numerical diffusion, the
Code:
C_VOF(c,lt)
does not necessary be zero even for a cell taken up by one specific phase. You should use
Code:
( fabs( C_VOF(c, lt) ) < 1.e-3 ) /*you might need to adjust the tolerance 1.e-3*/
instead.
blackmask is offline   Reply With Quote

Old   May 30, 2013, 20:58
Default
  #5
Member
 
Shashank
Join Date: Apr 2011
Posts: 74
Rep Power: 14
shashank312 is on a distinguished road
Could I loop over threads of a specific phase? I thought it generally gives a segmentation error. Let me try it anyway.

About the tolerance, I agree. The volume fraction may not necessarily be zero or one even though the cell is completely filled with a specific phase. Thanks for the extra tip.
shashank312 is offline   Reply With Quote

Old   May 30, 2013, 21:08
Default
  #6
Member
 
Shashank
Join Date: Apr 2011
Posts: 74
Rep Power: 14
shashank312 is on a distinguished road
BTW, if I use
begin_c_loop_all(c,lt)
instead of
begin_c_loop_all(c,t)
how can I access the mixture threads for temperature macro
C_T(c,t)
shashank312 is offline   Reply With Quote

Old   May 30, 2013, 21:10
Default
  #7
Member
 
Shashank
Join Date: Apr 2011
Posts: 74
Rep Power: 14
shashank312 is on a distinguished road
Could I use C_T(c,mt)? Here "mt" is a new thread different from "t", and specified as:

Thread *mt = THREAD_SUPER_THREAD(lt);
shashank312 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 16:33
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 05:03
Convergence moving mesh lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 15:09
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 03:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


All times are GMT -4. The time now is 05:45.