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

zero second derivative boundary condition

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 5, 2014, 22:37
Default zero second derivative boundary condition
  #1
Member
 
Lianhua Zhu
Join Date: Aug 2011
Location: Wuhan, China
Posts: 35
Rep Power: 14
zhulianhua is on a distinguished road
Hi! FOAMers

Does anybody know if there's a zero second derivative boundary condition already implemented in OpenFOAM ?

I wan't to let the fvc::grad() correct the gradient at the boundary face cells to the gradient of the next layer cells.


Regards,

Lianhua
zhulianhua is offline   Reply With Quote

Old   August 6, 2014, 13:59
Default
  #2
New Member
 
CZ
Join Date: Jun 2013
Posts: 9
Rep Power: 12
titanchao is on a distinguished road
Hey, Lianhua,

Could you define a scalar field that is the first-order derivative of your desired function and then impose zeroGradient BC of the defined scalar field.

Best of luck!
titanchao is offline   Reply With Quote

Old   August 6, 2014, 23:20
Default
  #3
Member
 
Lianhua Zhu
Join Date: Aug 2011
Location: Wuhan, China
Posts: 35
Rep Power: 14
zhulianhua is on a distinguished road
Hi, titanchao!

Thanks for you kindly reply.

Yes, I use the lsGrad scheme, from its source code, I can see the returned gradient field is initialized with "zeroGradient" by default.

Suppose I want get the gradient g of the scalar volume field s, I do not care about the boundary value of the g, but the it's value in the cells adjacent to the boundary faces.

The problem is: When lsGrad scheme process these cells adjacent to the boundary, it uses the values of s on the boundary faces.
However I want the boundary values of s to be linearly extrapolated form the inner cells (which is zero 2nd order B.C. I think).

Until now, I can't find how to extrapolate linearly the inner cell values of s to the boundary faces.
I found some people asked this question on this forum, but the replies say it's non-sense to do that physically!


Regards,

Lianhua

Quote:
Originally Posted by titanchao View Post
Hey, Lianhua,

Could you define a scalar field that is the first-order derivative of your desired function and then impose zeroGradient BC of the defined scalar field.

Best of luck!
zhulianhua is offline   Reply With Quote

Old   August 7, 2014, 11:29
Default
  #4
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 197
Blog Entries: 1
Rep Power: 18
egp is on a distinguished road
Quote:
Originally Posted by zhulianhua View Post
Hi! FOAMers

Does anybody know if there's a zero second derivative boundary condition already implemented in OpenFOAM ?

I wan't to let the fvc::grad() correct the gradient at the boundary face cells to the gradient of the next layer cells.


Regards,

Lianhua
Search before asking! While the SNR has gotten very low on the forum, there are some nuggets out there for you to find!

http://www.cfd-online.com/Forums/ope...condition.html
egp is offline   Reply With Quote

Old   August 7, 2014, 22:13
Default
  #5
Member
 
Lianhua Zhu
Join Date: Aug 2011
Location: Wuhan, China
Posts: 35
Rep Power: 14
zhulianhua is on a distinguished road
Thanks, Eric!

The post in the link is what I'm looking for and really a nugget!

Sorry, I should have done a good search work.


Regards,

Lianhua
zhulianhua 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
No-slip condition for non-resolved boundary layer in open channel banks Lupocci Main CFD Forum 1 January 17, 2013 03:11
domain imbalance for enrgy equation happy CFX 14 September 6, 2012 01:54
inlet velocity boundary condition murali CFX 5 August 3, 2012 08:56
External Radiation Boundary Condition for Grid Interface CFD XUE FLUENT 0 July 9, 2010 02:53
External Radiation Boundary Condition (Two sided wall), Grid Interface CFD XUE FLUENT 0 July 8, 2010 06:49


All times are GMT -4. The time now is 15:37.