 |
 |
 |
|
|
|
[Sponsors] | ||||
June 19, 2010, 11:28
|
twoLiquidMixingFoam
|
#1 |
|
Member
Maolong LIU
Join Date: Apr 2010
Location: USA
Posts: 31
Rep Power: 16  |
Hi, because the twoLiquidMixingFoam in OpenFoam is a transent solver. I was trying to change it to a steady state solver according to simpleFoam and twoLiquidMixingFoam.
Now the new solver can run successfully and it seems that the calculation result is correct. I compared u with experiment result and previous transient result. But the problem is that there is a converge problem. After run a long time, pressure, k, epsilon etc. are all converged, but onlu u does not converge. Time = 30575 DILUPBiCG: Solving for alpha1, Initial residual = 5.30999e-08, Final residual = 5.30999e-08, No Iterations 0 Phase 1 volume fraction = 0.500002 Min(alpha1) = 0 Max(alpha1) = 1 DILUPBiCG: Solving for Ux, Initial residual = 0.135853, Final residual = 1.4931e-08, No Iterations 6 DILUPBiCG: Solving for Uy, Initial residual = 0.0184276, Final residual = 1.39218e-08, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 3.33005e-07, Final residual = 3.33005e-07, No Iterations 0 time step continuity errors : sum local = 6.94663e-09, global = -1.76227e-11, cumulative = 2.20847e-08 DILUPBiCG: Solving for epsilon, Initial residual = 9.27893e-08, Final residual = 9.27893e-08, No Iterations 0 DILUPBiCG: Solving for k, Initial residual = 3.85187e-07, Final residual = 1.09717e-08, No Iterations 1 ExecutionTime = 22.01 s ClockTime = 22 s Can anyone give me any advice? Thanks. |
|
|
|
|
|
August 9, 2010, 06:20
|
|
#2 |
|
Senior Member
Balkrishna Patankar
Join Date: Mar 2009
Location: Pune
Posts: 123
Rep Power: 17 |
Do try changing the relaxation values ....
|
|
|
|
|
|
|
August 9, 2010, 06:31
|
|
#3 | |
|
Member
Maolong LIU
Join Date: Apr 2010
Location: USA
Posts: 31
Rep Power: 16 |
Thank you for your reply. I modified the U equation
UEqn== -fvc::grad(p) to UEqn== fvc::reconstruct ( ...... ) according to the twoLiquidMixingFoam, and now this problem has been solved. Now I am trying to speed up the convergence because my mesh number is really huge. Now the procedure of this solver is first solver U eqn and then p eqn and then alpha1 eqn. Also I change the relax method of alpha eqn just like p eqn (restore alpha1 and relax alpha1 instead alpha1 eqn). How do you think of it? Thank you! Quote:
|
||
|
|
|
||
|
August 9, 2010, 06:45
|
|
#4 |
|
Senior Member
Balkrishna Patankar
Join Date: Mar 2009
Location: Pune
Posts: 123
Rep Power: 17 |
thats a nice modification .... relaxing alpha1 after every loop does converge faster ....
I am working on something similar and facing a ver different problem ... pls do help ... link to thread discussion here .: http://www.cfd-online.com/Forums/ope...tml#post270776 |
|
|
|
|
|
|
August 9, 2010, 06:52
|
|
#5 | |
|
Member
Maolong LIU
Join Date: Apr 2010
Location: USA
Posts: 31
Rep Power: 16 |
If you don't mind, please send me your code, so I can help to check.
Quote:
|
||
|
|
|
||
|
August 9, 2010, 06:56
|
|
#6 |
|
Senior Member
Balkrishna Patankar
Join Date: Mar 2009
Location: Pune
Posts: 123
Rep Power: 17 |
whats ur email ???
|
|
|
|
|
|
|
August 9, 2010, 06:58
|
|
#7 |
|
Member
Maolong LIU
Join Date: Apr 2010
Location: USA
Posts: 31
Rep Power: 16 |
||
|
|
|
|
|
August 11, 2010, 02:23
|
|
#8 |
|
Senior Member
Balkrishna Patankar
Join Date: Mar 2009
Location: Pune
Posts: 123
Rep Power: 17 |
In the twoLiquidMixingFoam what exactly is alpha1 ?? according to the o/p statement ,
it represents volume fraction , but how is the conservation equation formed on the basis of volume fraction ? i mean conservation equation is : d/dt(rho_fluid*mass_K) + divergence(rho_fluid*U*mass_K)= convection + source terms .... where rho_fluid is density of the mixture .... mass_K is the mass fraction of the Kth component .... .... How is this written in terms of the volume fraction alpha ???? |
|
|
|
|
|
|
August 11, 2010, 05:05
|
|
#9 |
|
Senior Member
Balkrishna Patankar
Join Date: Mar 2009
Location: Pune
Posts: 123
Rep Power: 17 |
got it ... the formulation is correct ....
|
|
|
|
|
|
|
October 27, 2010, 22:26
|
|
#10 |
|
New Member
Join Date: Jul 2010
Posts: 11
Rep Power: 16 |
Has anyone seen negative values for alpha1 in twoLiquidMixingFoam? If so, any idea what changes I can make to correct the problem? I'm runnig a simulation of a pipe with RE around 2000 that has water and a higher viscosity fluid (~25cP). Any help is appreciated.
Thanks in advance! Phase 1 volume fraction = 0.000377962 Min(alpha1) = -2.37485 Max(alpha1) = 1 Last edited by phinallydone; October 27, 2010 at 22:48. |
|
|
|
|
|
|
October 28, 2010, 03:39
|
|
#11 |
|
Senior Member
Balkrishna Patankar
Join Date: Mar 2009
Location: Pune
Posts: 123
Rep Power: 17 |
that is not possible .... ur solution will diverge ....
|
|
|
|
|
|
|
October 28, 2010, 12:00
|
|
#12 |
|
New Member
Join Date: Jul 2010
Posts: 11
Rep Power: 16 |
That's just it... it's not diverging. I've tried refining the mesh, relaxing alpha1, and adjustnig tols. Any ideas?
|
|
|
|
|
|
|
November 12, 2012, 18:25
|
|
#13 |
|
New Member
Kia
Join Date: Jun 2011
Posts: 6
Rep Power: 15 |
Hey Maolong,
Can you tell me how you changed the solver to steady state. I am working on a similar problem and my delta T is really small and will like to make it steady state to avoid the Courant number restriction. |
|
|
|
|
|
 |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| SimpleFoam with body force gravity | alexice | OpenFOAM Running, Solving & CFD | 15 | June 16, 2022 15:43 |
| TwoLiquidMixingFoam | shawn | OpenFOAM Running, Solving & CFD | 65 | March 18, 2020 15:56 |
Linear Mode
