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

LES Outlet BC: Convection Boundary

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 31, 2013, 12:49
Default LES Outlet BC: Convection Boundary
  #1
Member
 
Anonymous
Join Date: Jan 2012
Location: Canada
Posts: 49
Rep Power: 5
Industrial_CFD is on a distinguished road
Hi Everyone,

I was wondering is the convection boundary condition for LES suggested by Peric and Ferzinger (d(phi)/dt +U*d(phi)/dx=0) has been implemented in OF 2.1

Previous threads have indicated that modifications have been made to OF 1.6; but since it would be a very useful BC, I thought that maybe it has been implemented as part of the standard release.

Currently I have a simulation blowing up with a constant pressure outlet, or an inletOutlet BC. I was would like to try a more physically realistic BC that taken the transient nature of LES into account.

Cheers
Industrial_CFD is offline   Reply With Quote

Old   January 31, 2013, 18:59
Default
  #2
Senior Member
 
Fumiya Nozaki
Join Date: Jun 2010
Location: Yokohama, Japan
Posts: 170
Rep Power: 7
fumiya is on a distinguished road
Hi,

I think the advective BC is what you are looking for.

Fumiya
fumiya is offline   Reply With Quote

Old   January 31, 2013, 23:01
Default reply
  #3
Member
 
Anonymous
Join Date: Jan 2012
Location: Canada
Posts: 49
Rep Power: 5
Industrial_CFD is on a distinguished road
Hi Fumiya,

From the OpenFOAM documentation advection BC is d(phi)/dt =0. So this isn't really true, since mass conservation isn't achieved.
Industrial_CFD is offline   Reply With Quote

Old   February 1, 2013, 04:45
Default
  #4
Senior Member
 
Lieven
Join Date: Dec 2011
Location: Mol, Belgium
Posts: 295
Rep Power: 13
Lieven will become famous soon enough
You could also consider using fixedMeanValue to set a surface averaged pressure. This causes less influence on the velocity field and might solve the 'blowing up' problem your having. Of course you should verify whether this BC is physically applicable to your case...

Cheers,

L
Lieven is offline   Reply With Quote

Old   February 1, 2013, 12:08
Default Reply
  #5
Member
 
Anonymous
Join Date: Jan 2012
Location: Canada
Posts: 49
Rep Power: 5
Industrial_CFD is on a distinguished road
Thanks Lieven, I will give that a try.
I am really confused about why this simulation is being such a bugger....LES I guess!
Industrial_CFD is offline   Reply With Quote

Old   February 1, 2013, 19:55
Default Reply
  #6
Member
 
Anonymous
Join Date: Jan 2012
Location: Canada
Posts: 49
Rep Power: 5
Industrial_CFD is on a distinguished road
Hi Lieven,

Is there anything special I have to do (e.g.download and compile) that BC?

Cheers: Adam

p.s. trying pimpleFoam to add stability...still no luck, solution is diverging.
Industrial_CFD is offline   Reply With Quote

Old   February 2, 2013, 13:32
Default
  #7
Senior Member
 
Lieven
Join Date: Dec 2011
Location: Mol, Belgium
Posts: 295
Rep Power: 13
Lieven will become famous soon enough
You can find a link to the sources in the following discussion:
Pressure outlet boundary condition

Cheers,

L
Lieven is offline   Reply With Quote

Old   February 5, 2013, 00:02
Default Update
  #8
Member
 
Anonymous
Join Date: Jan 2012
Location: Canada
Posts: 49
Rep Power: 5
Industrial_CFD is on a distinguished road
Hi Guys,

I implemented the meanFixedValue pressure outlet BC. I am getting what I think is strange behaviour, check my output. Pressure residual blows up but Co number is stable...anyone seen this before?




Cheers
Attached Files
File Type: txt log.txt (33.9 KB, 13 views)
Industrial_CFD is offline   Reply With Quote

Old   February 5, 2013, 23:26
Default Question
  #9
Member
 
Anonymous
Join Date: Jan 2012
Location: Canada
Posts: 49
Rep Power: 5
Industrial_CFD is on a distinguished road
Still having trouble figuring out why this is happening.
The same solution settings work perfectly in Fluent, yet in OpenFoam the pressure residual stabilizes around 0.85, even when the Co number is stable (around 0.035).

Any ideas why? Has anyone else seen anything similar?

Cheers
Industrial_CFD is offline   Reply With Quote

Old   February 6, 2013, 04:01
Default
  #10
Senior Member
 
Lieven
Join Date: Dec 2011
Location: Mol, Belgium
Posts: 295
Rep Power: 13
Lieven will become famous soon enough
Very strange indeed.
Unfortunately, I've never seen this type of behaviour before and I don't really have a clue what could be wrong... What are the other BC for the pressure?

Cheers,

L
Lieven is offline   Reply With Quote

Old   February 6, 2013, 07:11
Default Reply
  #11
Member
 
Anonymous
Join Date: Jan 2012
Location: Canada
Posts: 49
Rep Power: 5
Industrial_CFD is on a distinguished road
The outlet BC for pressure is fixedMean Value; meanValue 0;

It's really strange. I am going to post all of my dictionaries. When this type of thing happens it is usually some fundamental mistake...but I can't see it.

I really appreciate all of the help. Its hard to make a direct comparison between FLUENT and OpenFOAM, but I find it strange that OpenFOAM would be giving such odd behavior.
Attached Files
File Type: txt nuTilda.txt (1.4 KB, 7 views)
File Type: txt nuSgs.txt (1.2 KB, 9 views)
File Type: txt k.txt (1.5 KB, 9 views)
File Type: txt B.txt (1.4 KB, 4 views)
File Type: txt U.txt (1.6 KB, 18 views)
Industrial_CFD is offline   Reply With Quote

Old   February 6, 2013, 07:12
Default Reply
  #12
Member
 
Anonymous
Join Date: Jan 2012
Location: Canada
Posts: 49
Rep Power: 5
Industrial_CFD is on a distinguished road
and also the 'p' dictionary.

Cheers
Attached Files
File Type: txt p.txt (1.2 KB, 29 views)
Industrial_CFD is offline   Reply With Quote

Old   February 6, 2013, 16:05
Default Update
  #13
Member
 
Anonymous
Join Date: Jan 2012
Location: Canada
Posts: 49
Rep Power: 5
Industrial_CFD is on a distinguished road
This is really confusing me. I have (once again) run the simulation in fluent, using piso, second order discretization in space and time; same initial conditions, same BC for U and p (constant V in, zero gradient out, constant average p out). Fluent is stable, converges nicely, doesn't blow up.

Openfoam won't converge the pressure residual; I have tried many different BCs. Could someone please please look at my dictionaries I posted above? There must be something fundamentally wrong. How can two solvers who have been compared many times be acting so differently on the same mesh?

I have been trying to solve this for so long that I need a second set of eyes...unless there is a bug in OpenFoam i'm not aware of.

Also, just to double check, when monitoring convergence, it's the "initial residual" that matters right? The "final residual" is related to the linear equation solver if I am correct?

Cheers
Industrial_CFD 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
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
outlet boundary conditions Xhoven OpenFOAM Running, Solving & CFD 3 June 11, 2011 19:33
Disturbed flow field at outlet boundary (Multiphase flow through pipe) Michiel CFX 17 April 21, 2010 10:14
VOF Outlet boundary condition in cfd - ace JM Main CFD Forum 0 December 15, 2006 09:07
Solver error message!!! IoSa CFX 1 September 14, 2006 04:48


All times are GMT -4. The time now is 05:16.