|
[Sponsors] |
September 22, 2006, 06:09 |
Conservation target in transient simulations
|
#1 |
Guest
Posts: n/a
|
Hello, does anbody know, where can I see, which equation does not fulfill the conservation target in my simulation???Is it possible to see it from the out file where linear solution is shown? (Without conservation taget checked 2 or 3 coefficient loops are used, but in case when I choose conservation target (even it does not matter is it 0.01 or 0.9 ) then the maximum number of coeff. loops always remain with the best wishes Target
|
|
September 22, 2006, 08:16 |
Re: Conservation target in transient simulations
|
#2 |
Guest
Posts: n/a
|
||
September 22, 2006, 08:37 |
Re: Conservation target in transient simulations
|
#3 |
Guest
Posts: n/a
|
I will try once again, may be somebody could help me to solve the problem...
The problem is that my transient heat transfer/CFD simulation does not satisfy the "conservation target". That I know because the solver goes through the maximal number of coeff. loops (set 10 or 20)even then when I set conservation target to maximal possible 0.99. In this simulation timestep is small enough, in case of not asking conservation target to be satisfied,- convergence of Max and rms residuals is very fast. It would be important for me to know, which equation exactly does not satisfy the conservation target. It could give me an answer why good convergent (when looking on rms or max residuals solution is not fully reliable). |
|
September 22, 2006, 08:43 |
Re: Conservation target in transient simulations
|
#4 |
Guest
Posts: n/a
|
Post a copy of your command file.
|
|
September 22, 2006, 08:55 |
Re: Conservation target in transient simulations
|
#5 |
Guest
Posts: n/a
|
Hi, it is various pages... I am speaking about that part in definition exactly.
SOLVER CONTROL: ADVECTION SCHEME: Option = High Resolution END CONVERGENCE CONTROL: Maximum Number of Coefficient Loops = 8 Timescale Control = Coefficient Loops END CONVERGENCE CRITERIA: Conservation Target = 0.99 Residual Target = 1e-06 Residual Type = MAX END TRANSIENT SCHEME: Option = Second Order Backward Euler Timestep Initialisation Option = Previous Timestep |
|
September 22, 2006, 09:26 |
Re: Conservation target in transient simulations
|
#6 |
Guest
Posts: n/a
|
Make the following changes:
CONVERGENCE CONTROL: Maximum Number of Coefficient Loops = 8 -> Change to 3 to 5 Timescale Control = Coefficient Loops END CONVERGENCE CRITERIA: Conservation Target = 0.99 <- Remove this entirely. This was the cause of your problem. Understand what it does before trying to use it again. Residual Target = 1e-06 <- Set too low, change to (1e-5 - 1e-4) Residual Type = MAX <- RMS may be more appropriate for your application depending on your mesh quality. END |
|
September 22, 2006, 09:55 |
Re: Conservation target in transient simulations
|
#7 |
Guest
Posts: n/a
|
You can monitor equation imbalances in the solver manager... make a new monitor and look at the 'Imbalance' fields.
You can also add the expert parameter 'monitor totals=t' to get the data dumped to your out file every timestep. |
|
September 22, 2006, 09:59 |
Re: Conservation target in transient simulations
|
#8 |
Guest
Posts: n/a
|
Thank you Joe for suggestion. I know that it is working without conservation target and with rms is converging very well. Unfortunately solution becomes not physical - not stable and that is why I wanted to keep balances in each inner loop and started to experiment with changing the value of conservation target - 0.01 did not work so I tried larger numbers and none of it worked (solver reached the number of maximal coeff. loop and continued the next timestep).
Yes, that is exactly what I am interested in - to understand the role of conservation target to find, which equation does not allow to satisfy the balances inside each inner loop of transient simulation. |
|
September 22, 2006, 10:07 |
Re: Conservation target in transient simulations
|
#9 |
Guest
Posts: n/a
|
Thank you Phil!!! I plotted equation imbalances! I can see that problem is with P-mas imbalance - after some iterations it is just equal to 1.
I have a fixed amount of fluid without inlet or outlet but I calculate transport equation for virtual additional variable inside it to simulate some other processes. Could P-mas imbalance =1 because of that or I have to look for some other problem? |
|
September 22, 2006, 10:11 |
Re: Conservation target in transient simulations
|
#10 |
Guest
Posts: n/a
|
If your solution is becoming unphysical it is much more likely that something else is causing the problem, rather than just poor equation imbalances.
|
|
September 22, 2006, 10:55 |
Re: Conservation target in transient simulations
|
#11 |
Guest
Posts: n/a
|
I call unphysical the situation when the solution, which is expected to be axisymmetrical becomes non-axisymmetrical. In some cases it remains symmetric, in some not - that is why I wanted to look more careful on balances...
Is there a possibility to apply conservation target to each equation separately (not to apply for one equation)? |
|
September 24, 2006, 17:14 |
Re: Conservation target in transient simulations
|
#12 |
Guest
Posts: n/a
|
Yes, I think that is the problem. I am quite sure you can define conservation targets of each equation. Either in Pre, otherwise in the command editor.
Another suggestion is to use double precision. This might help when you want to achieve deep convergence with large and small numbers.......... Astird |
|
September 25, 2006, 08:20 |
Re: Conservation target in transient simulations
|
#13 |
Guest
Posts: n/a
|
My guess is that P-mass equation imbalance =1 is because of density temperature dependence, included in the model.
My present question is - does anybody knows how is it possible to set conservation target for each equation seperately? |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Building OpenFOAM1.7.0 from source | ata | OpenFOAM Installation | 46 | March 6, 2022 14:21 |
Best practice for transient simulations? | siw | CFX | 5 | October 30, 2010 06:45 |
OpenFOAM14 for Mac OSX Darwin 104 | gschaider | OpenFOAM Installation | 118 | July 20, 2008 06:19 |
Momentum conservation in transient simulations | Latslosh | CFX | 0 | June 18, 2008 10:43 |
How are a transient simulations performed?? | Luk | Main CFD Forum | 0 | October 19, 2007 11:09 |