CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Monitoring residuals , is it safe ?

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

Like Tree13Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 13, 2017, 15:18
Default
  #21
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by BlnPhoenix View Post
So did you find some average temperature in your experiment? Have you measured exactly after one hour? Try to describe your experimental findings. Is the temperature constant over time, if yes what do you expect if you would measure after five minutes? Same temp as after one hour?

If you cannot afford the simulation to run to one hour, i would concider redoing your measurements and run it till you have experimental data to compare.

Mass imbalance means that the mass flow in the systems equals the mass flow out of the system. It's very important to ensure mass balance in every time step. So please check for that and as mentioned set your target residual values accordingly. Also for temperature and velocities as i described above. Your residuals sound good but try to look for the actual change in values to see what they actually mean.
My case is a dynamic and thermal study of air jets .Temperature was not constant over time because i was dealing with hot air jets in a room, and with time of corse the temperature increases , for each station ( axialy and radially) i mesured the temperature at that station to calculate the reduced temperature ( ti- ta)/(tmax- ta) where ti is the temperature at each station , ta is the ambiant temperature and tmax is the maximum température at the inlet. Until now i didnt find a numerical way (in postprocessing) to get the température the same way like i got it in experiance. Because experimantaly the mesurements of temperature was not with a constant time step, it depended to when the velocity on the device's screen became constant, so i waited until i get a velocity mesurements then i get a temperature mesurement (ti and ta) .
If there is no solution, i think i will just do a dynamic study, because it is impossible for me experimentaly to get température mesurements in a constant time step, because i dont know if in that moment the velocity will be stable on the device's screen.

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid is offline   Reply With Quote

Old   July 13, 2017, 16:48
Default
  #22
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Your case will be stiff using a compressible flow model. Are you preconditioning the equations? Could you switch to a pure incompressible solver? Try to evaluate the normalized residuals
I did a small search about your questions, well first, fluent hide the different solvers from the users, when you select different density specifications, it will activate different solvers. I am setting density constant where there is other choices like ideal gas, incompressible-ideal-gas..... so the i guess that i am using an incompressible solver.
For the preconditioning, what i understand is that it has un effect of reducing the acoustic speed of the system to accelerate the convergence in incompressible flows as the mach number goes to zero wright? But i dont know how to put that in the software

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid is offline   Reply With Quote

Old   July 13, 2017, 16:55
Default
  #23
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by medaouarwalid View Post
I did a small search about your questions, well first, fluent hide the different solvers from the users, when you select different density specifications, it will activate different solvers. I am setting density constant where there is other choices like ideal gas, incompressible-ideal-gas..... so the i guess that i am using an incompressible solver.
For the preconditioning, what i understand is that it has un effect of reducing the acoustic speed of the system to accelerate the convergence in incompressible flows as the mach number goes to zero wright? But i dont know how to put that in the software

Sent from my F1f using CFD Online Forum mobile app

http://www.afs.enea.it/project/neptu...th/node378.htm

but you can use also the NITA
medaouarwalid likes this.
FMDenaro is offline   Reply With Quote

Old   July 13, 2017, 17:09
Default
  #24
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I am using a bounded second order implicit for transient formulation in a pressure based solver

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid is offline   Reply With Quote

Old   July 13, 2017, 17:29
Default
  #25
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Okay, now i set non iterative time advancement. Any other quggestions befor i restart calculation sir ?

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid is offline   Reply With Quote

Old   July 13, 2017, 17:36
Default
  #26
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by medaouarwalid View Post
Okay, now i set non iterative time advancement. Any other quggestions befor i restart calculation sir ?

Sent from my F1f using CFD Online Forum mobile app

Have you infos about a physical initial condition? Otherwise, if you set an arbitrary initial condition you have to wait enough time to let the solution correlating in a physical way from the arbitrary numerical one.
medaouarwalid likes this.
FMDenaro is offline   Reply With Quote

Old   July 13, 2017, 17:50
Default
  #27
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Have you infos about a physical initial condition? Otherwise, if you set an arbitrary initial condition you have to wait enough time to let the solution correlating in a physical way from the arbitrary numerical one.
No, only boundary conditions. About residaul normalize, i let the default iteration number (5) and normalization factor (0) for all ???

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid is offline   Reply With Quote

Old   July 13, 2017, 18:17
Default
  #28
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Divergence detected in AMG solver: x-momentum

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid is offline   Reply With Quote

Old   July 13, 2017, 22:16
Default
  #29
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by medaouarwalid View Post
Divergence detected in AMG solver: x-momentum

Sent from my F1f using CFD Online Forum mobile app

This must be NITA of fluent. It is its main problem.
medaouarwalid likes this.
arjun is offline   Reply With Quote

Old   July 14, 2017, 03:56
Default
  #30
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by medaouarwalid View Post
Divergence detected in AMG solver: x-momentum

Sent from my F1f using CFD Online Forum mobile app
This is a very general advice and can be due to several issues...

Does it appear after the first time step or during the transient?
It could be to an error in the setting of the BC.s or a numerical instability ...
medaouarwalid likes this.
FMDenaro is offline   Reply With Quote

Old   July 14, 2017, 04:00
Default
  #31
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Are you setting the fractional step or the PISO? Did you set a second order time-advancing?
medaouarwalid likes this.
FMDenaro is offline   Reply With Quote

Old   July 14, 2017, 11:44
Default
  #32
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
This is a very general advice and can be due to several issues...

Does it appear after the first time step or during the transient?
It could be to an error in the setting of the BC.s or a numerical instability ...
It appeared after the 4th time step

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid is offline   Reply With Quote

Old   July 14, 2017, 11:51
Default
  #33
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Are you setting the fractional step or the PISO? Did you set a second order time-advancing?
I was using The piso one, i will try FSM then i will tell you + yes bounded second order implicit

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid is offline   Reply With Quote

Old   July 14, 2017, 12:24
Default
  #34
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
- i used hybrid initialization , a warning displays
Warning: convergence tolerance of 1.000000e-06 not reached
during Hybrid Initialization.

-i started calculation
nb : there is also warnings from the begining of the calculation:
reversed flow in 4884 faces on pressure-outlet 6. reversed flow in 4884 faces on pressure-outlet 7.
- after incrementing to the 3rd time step :
# Divergence detected in AMG for x-momentum: protective actions enabled!
# Divergence detected in AMG for x-momentum, temporarily solve with BCGSTAB!
Error: Divergence detected in AMG solver: x-momentum
Error Object: #f
medaouarwalid is offline   Reply With Quote

Old   July 14, 2017, 13:03
Default
  #35
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
my advise is that run first few time steps with implicit method without worry about cost.

This is reason, FVUS have a setting that says fraactional solver after this many time steps.

In the start it is not converged and NITA has problems.
medaouarwalid likes this.
arjun is offline   Reply With Quote

Old   July 14, 2017, 13:35
Default
  #36
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by arjun View Post
my advise is that run first few time steps with implicit method without worry about cost.

This is reason, FVUS have a setting that says fraactional solver after this many time steps.

In the start it is not converged and NITA has problems.
What implicit methode ? I am already using it i guess

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid is offline   Reply With Quote

Old   July 14, 2017, 13:49
Default
  #37
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Check your BC.s, maybe they do not satisfy the incompressibility constraint
FMDenaro is offline   Reply With Quote

Old   July 14, 2017, 13:51
Default
  #38
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Check your BC.s, maybe they do not satisfy the incompressibility constraint
How is so ? 7 Hot air jets blowing of 7 nozzles in atmospher in a room with velocity 7 m/s for each jet. I dont think it is cmplicated BCs .
The roof is wall, nozzles wall also, and the rest i set it pressure outlet because i wanted to limite the domaine, should i put the real dimensions of the room and set every thing as wall ?

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid is offline   Reply With Quote

Old   July 14, 2017, 13:56
Default
  #39
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by medaouarwalid View Post
How is so ? 7 Hot air jets blowing of 7 nozzles in atmospher in a room with velocity 7 m/s for each jet. I dont think it is cmplicated BCs .
The roof is wall, nozzles wall also, and the rest i set it pressure outlet because i wanted to limite the domaine, should i put the real dimensions of the room and set every thing as wall ?

Sent from my F1f using CFD Online Forum mobile app

Pressure outlet must allow the mass flow rate entering from the nozzles to exit. If that does not happen, the continuity equation is not satisfied, the mass increases and the error in the divergence-free constraint acts as a source driving to divergence.
medaouarwalid likes this.
FMDenaro is offline   Reply With Quote

Old   July 14, 2017, 14:07
Default
  #40
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Pressure outlet must allow the mass flow rate entering from the nozzles to exit. If that does not happen, the continuity equation is not satisfied, the mass increases and the error in the divergence-free constraint acts as a source driving to divergence.
i noticed experimentaly that velocitie at let saying 2 meters in axial direction is close to zero , So to reduced the dimensions of the room i set pressure outlet at 2 meter instead of putting the real dimension of the room and set it as wall. I will change the geometri by setting the real dimensions and real boundary conditions . After that I will share with you the result.
A question, could the temperature of the jet in my case have an influence in compressibility ?

Sent from my F1f using CFD Online Forum mobile app
medaouarwalid 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
ANSYS Workbench on "Uncertified" Linux Distros hsr CFX 289 April 20, 2023 10:23
Monitoring Residuals - Initial or Final? BigPapi34 OpenFOAM Running, Solving & CFD 1 November 3, 2014 06:31
Monitoring Residuals - Initial or Final? BigPapi34 Main CFD Forum 0 November 3, 2014 06:23
Monitoring residuals in Workbench-Fluent using Remote Solver fadiga Main CFD Forum 2 October 16, 2013 16:43
judging convergence through residuals MachZero Main CFD Forum 7 December 25, 2012 13:18


All times are GMT -4. The time now is 18:14.