CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Suspicious results

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

Like Tree4Likes
  • 1 Post By wyldckat
  • 1 Post By fportela
  • 1 Post By fportela
  • 1 Post By fportela

Reply
 
LinkBack Thread Tools Display Modes
Old   September 1, 2013, 11:42
Default Suspicious results
  #1
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
why this suspicious result in p field has obtained?
what may be the cause?what experiences are there?
Attached Images
File Type: jpg suspicious2.jpg (24.8 KB, 22 views)
File Type: jpg suspicious3.jpg (14.2 KB, 20 views)
File Type: jpg suspicious.jpg (21.0 KB, 19 views)
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.

Last edited by immortality; September 1, 2013 at 14:30.
immortality is offline   Reply With Quote

Old   September 1, 2013, 13:29
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,312
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Uhm... I already answered to you here: How to visualize the difference of fields btween two simulation post #28.
immortality likes this.
wyldckat is offline   Reply With Quote

Old   September 2, 2013, 11:38
Default
  #3
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
Hi Bruno
above images was for kOmegaSST case,but for the case with realizableKE the results are smooth as attached at the same snapshot time.
maybe this model has damped the variations more strongly?
this site should be renamed to:www.Bruno-online.com!
Attached Images
File Type: jpg p.jpg (10.7 KB, 10 views)
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   September 2, 2013, 12:02
Default
  #4
Member
 
Felipe Portela
Join Date: Dec 2012
Location: London
Posts: 58
Rep Power: 4
fportela is on a distinguished road
Hi Eshan,

I think I had some problem similar to this in rhoCentralFoam running on a supersonic channel flow: rhoCentralFoam for channel flow

By applying near wall refinement I was able to reduce those strange waves. However, I was never able to get a nice turbulent flow...
immortality likes this.
fportela is offline   Reply With Quote

Old   September 2, 2013, 16:17
Default
  #5
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
Hi Felipe
thanks for consideration,
then you used LES and still this problem existed?
in the other run I did with realizableKE as you observe above,its smoother in high pressure parts of the tube.
so what may be the cause in your opinion?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   September 2, 2013, 16:26
Default
  #6
Member
 
Felipe Portela
Join Date: Dec 2012
Location: London
Posts: 58
Rep Power: 4
fportela is on a distinguished road
Honestly I do not know. I might have done some mistake setting up the case, but I highly doubt it since I was able to get reasonable results with rhoPimpleFoam.

My whole point of using rhoCentralFoam was to try and get the implicit dissipation to work as my SGS model, in the end, as I mentioned, the flow was rather laminar so adding an explicit model did not change anything.

My case was supposed to have almost uniform pressure due to the presence of the forcing term, but still I got the wavy pattern as you see in my other post.

What exactly is your case setup? Maybe we are computing completely different things, rendering this discussion pointless :P

Quote:
Originally Posted by immortality View Post
Hi Felipe
thanks for consideration,
then you used LES and still this problem existed?
in the other run I did with realizableKE as you observe above,its smoother in high pressure parts of the tube.
so what may be the cause in your opinion?
fportela is offline   Reply With Quote

Old   September 2, 2013, 17:21
Default
  #7
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
Hi
my case is a tube like the tutorial of shockTube,but with unsteady boundaries in both sides,a shock originate like shock tube from one side and moves to the other(front) side.
but in times has shown in images above,shouldn't remain strong shocks.
I only changed the turbulent model,do you think its because realizableKE always use wall functions while in kOmegaSST case may it use values on the wall directly sometimes during the run?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   September 2, 2013, 17:28
Default
  #8
Member
 
Felipe Portela
Join Date: Dec 2012
Location: London
Posts: 58
Rep Power: 4
fportela is on a distinguished road
Hi Eshan,

The shockTube tutorial is 1D. The plots in the first post do not seem to represent a 1D case.

Quote:
Originally Posted by immortality View Post
Hi
my case is a tube like the tutorial of shockTube,but with unsteady boundaries in both sides,a shock originate like shock tube from one side and moves to the other(front) side.
but in times has shown in images above,shouldn't remain strong shocks.
I only changed the turbulent model,do you think its because realizableKE always use wall functions while in kOmegaSST case may it use values on the wall directly sometimes during the run?
fportela is offline   Reply With Quote

Old   September 3, 2013, 04:53
Default
  #9
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
Hi Felipe
no I said its general shape,I converted it to my problem needs,then my case is 2D as you told.has it any effect if I use cellMDLimited in grad schemes?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   September 3, 2013, 04:56
Default
  #10
Member
 
Felipe Portela
Join Date: Dec 2012
Location: London
Posts: 58
Rep Power: 4
fportela is on a distinguished road
what boundary conditions do you use for the top and bottom boundaries?

have you tried an clean case (without any limiters and whatsoever)? this should shed some light on what is going on.

Quote:
Originally Posted by immortality View Post
Hi Felipe
no I said its general shape,I converted it to my problem needs,then my case is 2D as you told.has it any effect if I use cellMDLimited in grad schemes?
fportela is offline   Reply With Quote

Old   September 3, 2013, 05:04
Default
  #11
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
upper and lower are walls but right and left are patches,yes I used schemes without any limiters till now,but maybe a limiter can result a smoother solution I think,whats your opinion?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   September 3, 2013, 05:11
Default
  #12
Member
 
Felipe Portela
Join Date: Dec 2012
Location: London
Posts: 58
Rep Power: 4
fportela is on a distinguished road
Hi Eshan,

By walls you mean you impose no-slip conditions right?

Have you tried refining the mesh (still without limiters) to see what happens? In my problem with rhoCentralFoam (which I posted above), I was able to remove the waves by applying near wall refinement and using a rather fine mesh.

I would see to try refining and identify if the pattern you see is related to the mesh size, or if it stays constant with refinement. In case of the former, then the issue is related with your numerical model (in which case limiters might help), in case of the later, the issue should be related to you physical model.

Keep me updated.

Cheers,
Felipe

Quote:
Originally Posted by immortality View Post
upper and lower are walls but right and left are patches,yes I used schemes without any limiters till now,but maybe a limiter can result a smoother solution I think,whats your opinion?
immortality likes this.
fportela is offline   Reply With Quote

Old   September 3, 2013, 11:30
Default
  #13
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
but if I make mesh finer near wall its y+ will be too low,I want it be more than 30 till I have a high Re mesh.

edit: yplus turns below 30 in many snapshot times like the attached one,so I use mutUSpaldingWallFunction that someone said to me that its fine for these situations,is it OK in your opinion?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.

Last edited by wyldckat; September 4, 2013 at 16:19.
immortality is offline   Reply With Quote

Old   September 3, 2013, 12:27
Default
  #14
Member
 
Felipe Portela
Join Date: Dec 2012
Location: London
Posts: 58
Rep Power: 4
fportela is on a distinguished road
My point in the previous post is that if you refine (to near DNS resolution) then you shouldn't need limiter or explicit models (someone correct me if I'm wrong). This is a nice way to isolate the problem, I would say.

Quote:
Originally Posted by immortality View Post
yplus turns below 30 in many snapshot times like the attached one,so I use mutUSpaldingWallFunction that someone said to me that its fine for these situations,is it OK in your opinion?
wyldckat likes this.

Last edited by wyldckat; September 4, 2013 at 16:20. Reason: fixed broken link in quote
fportela is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
lid driven cavity varying results yasmil OpenFOAM Running, Solving & CFD 1 July 12, 2013 09:57
CFD results not close to experimental results cider STAR-CCM+ 0 July 8, 2013 07:53
OpenFOAM - Validation of Results Ahmed OpenFOAM Running, Solving & CFD 9 June 22, 2011 18:59
Different Results from Fluent 5.5 and Fluent 6.0 Rajeev Kumar Singh FLUENT 6 December 19, 2010 12:33
Airfoil y+ Results Query asd Main CFD Forum 1 June 22, 2007 08:02


All times are GMT -4. The time now is 00:17.