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

Convective outlet boundary condition for Unsteady flows

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

Like Tree10Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   November 10, 2008, 21:12
Default I would not change many things
  #61
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 698
Rep Power: 12
msrinath80 is on a distinguished road
I would not change many things together. My recommendation is to take a case which gives you a fairly smooth profile for Cd/Cl and *only* change the outlet B/C to convectiveOutlet. If you see spikes now, it is very likely that this is due to the outlet B/C. If you don't see any spikes, go ahead and change the time discretization to "crankNicholson 1" and see what happens.

All said and done, my experience with OpenFOAM tells me that this is because of using "crankNicholson 1" in fvSchemes. For a more detailed exposition, read this[1] thread where Dr. Weller clearly explains why using "full crankNicholson" isn't really a good idea.


References:
[1] http://www.cfd-online.com/OpenFOAM_D...tml?1195511709
immortality likes this.
msrinath80 is offline   Reply With Quote

Old   November 11, 2008, 08:36
Default Hai pU|, I have run the sim
  #62
New Member
 
Join Date: Mar 2009
Location: adelaide, SA, Australia
Posts: 27
Rep Power: 8
mali is on a distinguished road
Hai pU|,

I have run the simulation again, and this time i was using backward second order time discretization scheme. No problems arise this time and i got smooth force oscillation.

I think you are right. This is not because of the convectiveOutlet, but more probably to 'full' crankNicholson scheme.

Thank you again.
__________________
mali
mali is offline   Reply With Quote

Old   November 21, 2008, 04:00
Default Hi all, In case if any one
  #63
Super Moderator
 
Takuya OSHIMA
Join Date: Mar 2009
Location: Niigata City, Japan
Posts: 518
Blog Entries: 1
Rep Power: 11
7islands is on a distinguished road
Hi all,

In case if any one of you needs a clarification for the discretizations and testing results of the convectiveOutlet bc code, I have written a quick summary [1]. The description starts from section 2.4 at page 5.

[1] http://oshima.eng.niigata-u.ac.jp/Op..._0287-open.pdf
7islands is offline   Reply With Quote

Old   November 21, 2008, 05:35
Default Wow, proceedings of Inter-nois
  #64
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: South Bend, IN, USA
Posts: 688
Blog Entries: 9
Rep Power: 12
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Wow, proceedings of Inter-noise 2008 (Shanghai), Should have pay you a visit, haha.
__________________
~
Daniel WEI
-------------
NatHaz Modeling Laboratory
Department of Civil & Environmental Engineering & Earth Sciences
University of Notre Dame, USA
Email || My Personal CFD Blog
lakeat is offline   Reply With Quote

Old   December 3, 2008, 07:44
Default Greetings all, As a final t
  #65
Member
 
Johan Lorentzon
Join Date: Mar 2009
Location: Lunds University, Sweden
Posts: 71
Rep Power: 14
pi06jl6 will become famous soon enough
Greetings all,

As a final touch in my forthcoming graduation, i am about to use convective outlet, this in order to minimize the corridor in my FSI program, but i cannot use wmake, i am guessing i am missing something but neverthenless the error message indicates no proper definitions, all of the error messages refers to "no matching function" followed by "candidates are", then i get like a few hundreds of lines, all related to fixedGradientPatchField. My educated guess is that wrong ld path is used, so here comes the settings in the option file: -lfiniteVolume. Anyone please help me?

Best Regards

Johan
pi06jl6 is offline   Reply With Quote

Old   December 3, 2008, 08:06
Default Hi Johan, It sounds like a co
  #66
Super Moderator
 
Takuya OSHIMA
Join Date: Mar 2009
Location: Niigata City, Japan
Posts: 518
Blog Entries: 1
Rep Power: 11
7islands is on a distinguished road
Hi Johan,
It sounds like a compilation error, not a linker error. Can you post a full build log?

Takuya
7islands is offline   Reply With Quote

Old   December 3, 2008, 20:10
Default Hi Johan, I saw the log via e
  #67
Super Moderator
 
Takuya OSHIMA
Join Date: Mar 2009
Location: Niigata City, Japan
Posts: 518
Blog Entries: 1
Rep Power: 11
7islands is on a distinguished road
Hi Johan,
I saw the log via email. I guess you are trying to compile one of the older versions. Can you try the newest one, 20070905 (see my post on Sep. 05, 2007) if my guess is correct?

Takuya
7islands is offline   Reply With Quote

Old   December 4, 2008, 10:11
Default Hi Johan, Thank you for the
  #68
Super Moderator
 
Takuya OSHIMA
Join Date: Mar 2009
Location: Niigata City, Japan
Posts: 518
Blog Entries: 1
Rep Power: 11
7islands is on a distinguished road
Hi Johan,

Quote:
Thank you for the answer, i am quite positive that i am missing something
like a flag or enviromental settings, since other people have used it. To
your question, of course, the latest as you suggested on the thread. Take
a look on the include -I, is that correct?

Best Regards

> SOURCE=convectiveOutletFvPatchFields.C ; g++ -m32 -Dlinux -DDP -Wall
> -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -O3
> -DNoRepository
> -ftemplate-depth-40
> -I/home/pi06jl6/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude
> -IlnInclude -I. -I/home/pi06jl6/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/
> lnInclude -fPIC -pthread -c $SOURCE -o
> Make/linuxGccDPOpt/convectiveOutletFvPatchFields.o
In fact, no. Since Make/options in the 20070905 distribution contains the following lines
<pre>EXE_INC = \
-I$(LIB_SRC)/finiteVolume/lnInclude \
-DWM_PROJECT_VERSION_MAJOR=`echo $(WM_PROJECT_VERSION) | awk 'BEGIN{FS="."}{print $$1;}'` \
-DWM_PROJECT_VERSION_MINOR=`echo $(WM_PROJECT_VERSION) | awk 'BEGIN{FS="."}{print $$2;}'` \
-DWM_PROJECT_VERSION_PATCH=`echo $(WM_PROJECT_VERSION) | awk 'BEGIN{FS="."}{if (length($$3) > 0) print $$3; else print "0";}'`</pre>their corresponding lines should appear in the log, and that's why I suspected if you were trying an old version. Can you check if your Make/options contains the lines? If the file does contain the lines, can you try compiling some other user libraries floating around in this forum or on the Wiki in order to see if your environment is properly set up?

Takuya

p.s. please keep the discussion on the forum whenever possible.
7islands is offline   Reply With Quote

Old   December 4, 2008, 15:52
Default Thank you, this solved my prob
  #69
Member
 
Johan Lorentzon
Join Date: Mar 2009
Location: Lunds University, Sweden
Posts: 71
Rep Power: 14
pi06jl6 will become famous soon enough
Thank you, this solved my problem!
pi06jl6 is offline   Reply With Quote

Old   December 13, 2008, 04:53
Default Hi all; I was running a 2-d
  #70
New Member
 
Join Date: Mar 2009
Location: adelaide, SA, Australia
Posts: 27
Rep Power: 8
mali is on a distinguished road
Hi all;

I was running a 2-d square cylinder at low reynolds number and defined the outlet as convectiveOutlet. However the results of the vorticity contour was not as i expected. Instead of Karman vortex shedding distrubution, i got weird vorticity contour near the outlet.

I have tried many type of spatial and temporal discretization schemes. But the results were still the same.


1) anybody know what was actually happened.
2) What type / how should i defined on the initial condition of the pressure at the outlet.

Thank you

-sukri-
__________________
mali
mali is offline   Reply With Quote

Old   December 15, 2008, 20:12
Default Hi all; I recently noticed
  #71
New Member
 
Join Date: Mar 2009
Location: adelaide, SA, Australia
Posts: 27
Rep Power: 8
mali is on a distinguished road
Hi all;

I recently noticed the mistake i have made when doing the convectiveOutlet simulation. I defined the startFrom in the controlDict as latestTime. Now, i have fixed it, and the results seem good for me.

Thank you Takuya for sharing the convective outlet code.

Thank you again

-sukri-
__________________
mali
mali is offline   Reply With Quote

Old   March 6, 2009, 06:33
Default Hi all I am a new user. I w
  #72
New Member
 
Brugiere Olivier
Join Date: Mar 2009
Posts: 27
Rep Power: 8
brugiere_olivier is on a distinguished road
Hi all

I am a new user. I was running a diffuser and defined the outlet as convectiveOutlet. But I've the same problem than Cedric. When I use the convectiveOutlet condition in parallel, I've this message :

[1] --> FOAM FATAL ERROR :
gradientInternalCoeffs cannot be called for a defaultFvPatchField (actual type convectiveOutlet)
on patch sortie of field U in file "/craya/data/brugiere/poursuitesortie/processor1/765.52/U"
You are probably trying to solve for a field with a default boundary condition.
[1]
[1] From function defaultFvPatchField<type>::gradientInternalCoeffs( ) const
[1] in file fields/fvPatchFields/basic/default/defaultFvPatchField.C at line 694.
[1]
FOAM parallel run exiting
[1]
[compute-0-0.local:08265] MPI_ABORT invoked on rank 1 in communicator MPI_COMM_WORLD with errorcode 1

[0]
[0] --> FOAM FATAL ERROR :
gradientInternalCoeffs cannot be called for a defaultFvPatchField (actual type convectiveOutlet)
on patch sortie of field U in file "/craya/data/brugiere/poursuitesortie/processor0/765.52/U"
You are probably trying to solve for a field with a default boundary condition.
[0]
[0] From function defaultFvPatchField<type>::gradientInternalCoeffs( ) const
[0] in file fields/fvPatchFields/basic/default/defaultFvPatchField.C at line 694.
[0]
FOAM parallel run exiting
[0]
[compute-0-8.local:15020] MPI_ABORT invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 1

Somebody have any ideas ?

Thank you

Olivier
brugiere_olivier is offline   Reply With Quote

Old   March 6, 2009, 06:39
Default Hi all I am a new user. I w
  #73
New Member
 
Brugiere Olivier
Join Date: Mar 2009
Posts: 27
Rep Power: 8
brugiere_olivier is on a distinguished road
Hi all

I am a new user. I was running a diffuser and defined the outlet as convectiveOutlet. But I've the same problem than Cedric. When I use the convectiveOutlet condition in parallel, I've this message :

[1] --> FOAM FATAL ERROR :
gradientInternalCoeffs cannot be called for a defaultFvPatchField (actual type convectiveOutlet)
on patch sortie of field U in file "/craya/data/brugiere/poursuitesortie/processor1/765.52/U"
You are probably trying to solve for a field with a default boundary condition.
[1]
[1] From function defaultFvPatchField<type>::gradientInternalCoeffs( ) const
[1] in file fields/fvPatchFields/basic/default/defaultFvPatchField.C at line 694.
[1]
FOAM parallel run exiting
[1]
[compute-0-0.local:08265] MPI_ABORT invoked on rank 1 in communicator MPI_COMM_WORLD with errorcode 1

[0]
[0] --> FOAM FATAL ERROR :
gradientInternalCoeffs cannot be called for a defaultFvPatchField (actual type convectiveOutlet)
on patch sortie of field U in file "/craya/data/brugiere/poursuitesortie/processor0/765.52/U"
You are probably trying to solve for a field with a default boundary condition.
[0]
[0] From function defaultFvPatchField<type>::gradientInternalCoeffs( ) const
[0] in file fields/fvPatchFields/basic/default/defaultFvPatchField.C at line 694.
[0]
FOAM parallel run exiting
[0]
[compute-0-8.local:15020] MPI_ABORT invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 1

Somebody have any ideas ?

Thank you

Olivier
brugiere_olivier is offline   Reply With Quote

Old   April 13, 2009, 23:26
Default
  #74
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: South Bend, IN, USA
Posts: 688
Blog Entries: 9
Rep Power: 12
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Quote:
Originally Posted by msrinath80 View Post
@Eugene:

I just pulled out an excerpt from a very recent paper in Journal of Computational Physics (Vol: 206, page 661), and I can see now what you mean by 'slightly confused'.

"...The Convective Boundary Condition (CBC) has been used by Najjar and Vanka [4] and Najjar and Balachandran [5] in computing uniform flow past a normal flat plate. Cheng and Armfield [6] have employed it in computing uniform, two-dimensional flow past a circular cylinder. It has also been employed by Pauley et al. [7], Arnal et al. [8], Sohankar et al. [9] and Biswas and co-workers [10] for the computations of uniform, two-dimensional flow past a square cylinder. However, it is felt that the CBC lacks in a proper physical basis for elliptic and parabolic problems and it is also somewhat awkward to implement. The quantity U in the Equation is loosely defined in the literature as different authors have defined it in a different manner [4-*7]. The value of U, which yields minimum distortion in the vorticity structure at the exit, has to be determined by trial..."

Sorry, please bear with me, I still can't get the point, Why convectiveOulet is prefered by many? Can anybody give me some good explanations?
Is it bcoz CBC can shorten the length towards outlet? Do you have a complete review of this bc that I may refer to?

Another question, what are the requirements for the other turbulence quantities at the outlet, such as nuSgs, if CBC is used, is there any caution I should note?

Thanks for your patience!
__________________
~
Daniel WEI
-------------
NatHaz Modeling Laboratory
Department of Civil & Environmental Engineering & Earth Sciences
University of Notre Dame, USA
Email || My Personal CFD Blog
lakeat is offline   Reply With Quote

Old   April 29, 2009, 03:39
Default
  #75
New Member
 
Join Date: Apr 2009
Posts: 26
Rep Power: 8
askjak is on a distinguished road
Hi everybody,

I am looking into getting a convective outlet boundary condition for use in LES.

I can see that Takuya have done a lot of great work creating the convectiveOutlet BC which is supported for older version OpenFOAM and some has managed to use it under OF 1.5.

In the release documentation for OF 1.5 it is stated that it includes a generalized advective outflow boundary condition which is the basis for for example the waveTransmissive BC.

My question is now the following: Has anyone tried to make a convective BC based on the generalized advective code included in OF 1.5?

/Ask
askjak is offline   Reply With Quote

Old   May 15, 2009, 03:54
Default
  #76
New Member
 
Join Date: Apr 2009
Posts: 26
Rep Power: 8
askjak is on a distinguished road
Hi again,

I am trying to create a convectiveOutlet BC which uses advectiveFvPatchField (see post above). The idea is to set a constant speed at the outlet by overriding advectionSpeed() from advectiveFvPatchField.

Everything is in place but I can't manage to return a constant scalarField (see below) due to a lacking overview of the OF code and not optimal c++ skills.

template<class Type> tmp<scalarField> convectiveOutletFvPatchField<Type>::advectionSpeed () const

Please help with this simple problem. When I have tested that everything works I will post the BC to the forum for everybody benefit.

Best regards,Ask
askjak is offline   Reply With Quote

Old   June 12, 2009, 08:17
Default Convective outlet based on generalized advective outlet
  #77
New Member
 
Join Date: Apr 2009
Posts: 26
Rep Power: 8
askjak is on a distinguished road
I figured out how to solve the above problem.

I have made a convective outlet BC based on the generalized advective bc included in OF 1.5. I have not done thorough testing and have assumed that the generalized advective BC has been implemented correctly. I have only tested with the LES solver oodles.

Attached is the code for the boundary condition and a simple example.

convectiveOutlet20090612.tar.gz

pitzDaily_convectiveOutlet.tar.gz

Let me know if there is any problems.

/Ask
cfdonline2mohsen likes this.
askjak is offline   Reply With Quote

Old   September 10, 2009, 21:28
Default
  #78
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: South Bend, IN, USA
Posts: 688
Blog Entries: 9
Rep Power: 12
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Thanks, Ask!
__________________
~
Daniel WEI
-------------
NatHaz Modeling Laboratory
Department of Civil & Environmental Engineering & Earth Sciences
University of Notre Dame, USA
Email || My Personal CFD Blog
lakeat is offline   Reply With Quote

Old   September 22, 2009, 21:24
Default
  #79
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: South Bend, IN, USA
Posts: 688
Blog Entries: 9
Rep Power: 12
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Hi Takuya, and hi Foamers,

A good day to you!
When I use O-type grids to simulate my flow past a circular cylinder, I cannot get good results using convectiveOutlet. The streamwise velocity (UMean) near the outlet drops greatly. So in the end I enlarge my domain, and using zeroGradient and now the results looks pretty good, but whenever I use convectiveoutlet, the outlet flow field was greatly disturbed.

I don't know who invented this boundary condition, and I don't quite understand why there'are many who speak good of this boundary condition. Does it make sense by saying?
Quote:
"This condition ensures that vortices can approach and pass
the outflow boundary without significant disturbances or reflections into the inner domain"
And I don't know whether there's something wrong in the code itself.

Is there something wrong on my side? What's your experience, can you justify your results using convectiveOutlet?
__________________
~
Daniel WEI
-------------
NatHaz Modeling Laboratory
Department of Civil & Environmental Engineering & Earth Sciences
University of Notre Dame, USA
Email || My Personal CFD Blog
lakeat is offline   Reply With Quote

Old   September 22, 2009, 22:14
Default
  #80
Super Moderator
 
Takuya OSHIMA
Join Date: Mar 2009
Location: Niigata City, Japan
Posts: 518
Blog Entries: 1
Rep Power: 11
7islands is on a distinguished road
Hi Daniel,

Is the outlet of your O-mesh perpendicular to the mean stream? convectiveOutlet tries to convect physical quantity in the perpendicular direction to the boundary. Or alternatively you might be able to try the advective BC in the standard OF release.

As written in README I don't use this particular convectiveOutlet code for any of CFD but this type of BC worked very well when I did unsteady flow-past-a-cylinder simulations with my in-house LES code in the past.

Takuya
7islands 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
Convective boundary condition andrea_barbera OpenFOAM Running, Solving & CFD 4 March 4, 2010 05:36
convective boundary condition Mani CFX 7 February 2, 2008 17:25
Convective Boundary Condition garni FLUENT 0 September 25, 2005 13:00
Convective boundary condition STN Main CFD Forum 5 May 29, 2002 08:47
convective boundary condition frederic felten Main CFD Forum 1 April 25, 2001 01:20


All times are GMT -4. The time now is 23:35.