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

Uniquesness of the solution

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 2 Post By ghorrocks
  • 1 Post By Opaque

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 26, 2021, 17:54
Default Uniquesness of the solution
  #1
Senior Member
 
Jiri
Join Date: Mar 2014
Posts: 218
Rep Power: 13
Jiricbeng is on a distinguished road
Hello,

this is a mathematical question but I would like to clarify that with respect to ocassional comutational difficulties.
Is the numerical solution of Reynolds averaged NS equations unique within prescribed tolerance? The issue we all sometimes encounter: running a steady state analysis being initialized from another case, and the solution does not develop as stable, oscillates and you cannot force it to get better, simply speaking. But by initializing from another case or no initialization leads to much different computational progress and different result in fact. So, may that be caused by the innate properties of the equations that the solution may ''fall'' locally into a state from which it cant move to the desired (''correct'') state? Can this be a property of the equations for given boundary conditions not only due to smootheness of discretization and convergence limit? Or is the solution unique and this is caused just by the numerical limits that cut the possibility for the solution to move into different state?
Jiricbeng is offline   Reply With Quote

Old   May 26, 2021, 18:55
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,703
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Quote:
Is the numerical solution of Reynolds averaged NS equations unique within prescribed tolerance?
Some times yes, sometimes no. It is flow dependant.

Your other comments are interesting. Yes, there is a numerical accuracy component to this. Numerical inaccuracy can cause all sorts of non-physical behaviour.

But the Navier Stokes equations allows:
* Steady state, unique solutions (eg lid driven cavity at low Reynolds number - it does not matter what initial condition you specify the final flow is the same and steady)
* Steady state, non-unique solutions (eg detached flow in a diffuser - it can attach to one wall or the other in the diffuser. Which wall it attaches to depends on the initial conditions or asymmetric forcing you use. The final result will be steady state, but it could be attached to the top or bottom wall)
* Transient, non-unique (Von-Karman vortex street, flow around a bluff body etc).

Some background reading for you:
https://en.wikipedia.org/wiki/Fluidics
https://www.google.com/search?q=flui...OZagWGeK_NQDkM
Opaque and aero_head like this.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
ghorrocks is offline   Reply With Quote

Old   May 27, 2021, 08:43
Default
  #3
Senior Member
 
Join Date: Jun 2009
Posts: 1,804
Rep Power: 32
Opaque will become famous soon enough
Glenn's examples are great.

However, I am a bit blunter on the answer. CFD is nothing more than a root finder on a complex non-linear system of equations (PDE+Algebraic ones).

Because it is a non-linear system, it is highly possible that there are multiple roots. In some flow regimes, the non-linearities are weak and the system is fairly linear, and why the linearized equations produce very accurate solutions.

On top of the above, the numerical method used to solve the system could introduce false roots. Also, some methods can approach one of the roots easier than others, or none at all (stiff systems).

My 2cents
aero_head likes this.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
Opaque is offline   Reply With Quote

Reply


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
Filtering DNS solutions Vs projecting filtered DNS solution juliom Main CFD Forum 5 May 19, 2016 16:06
Why 3D solid-pore geometry showing diverged solution? Sargam05 OpenFOAM 0 December 3, 2012 15:45
Analytic solution for 2D steady Euler equations jojo81 Main CFD Forum 0 October 15, 2012 12:05
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Discussion about Mesh independant solution Seb Main CFD Forum 13 May 22, 2001 13:37


All times are GMT -4. The time now is 03:19.