|
[Sponsors] |
August 5, 2014, 22:37 |
zero second derivative boundary condition
|
#1 |
Member
Lianhua Zhu
Join Date: Aug 2011
Location: Wuhan, China
Posts: 35
Rep Power: 14 |
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 |
|
August 6, 2014, 13:59 |
|
#2 |
New Member
CZ
Join Date: Jun 2013
Posts: 9
Rep Power: 12 |
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! |
|
August 6, 2014, 23:20 |
|
#3 |
Member
Lianhua Zhu
Join Date: Aug 2011
Location: Wuhan, China
Posts: 35
Rep Power: 14 |
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 |
|
August 7, 2014, 11:29 |
|
#4 | |
Senior Member
|
Quote:
http://www.cfd-online.com/Forums/ope...condition.html |
||
August 7, 2014, 22:13 |
|
#5 |
Member
Lianhua Zhu
Join Date: Aug 2011
Location: Wuhan, China
Posts: 35
Rep Power: 14 |
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 |
|
|
|
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 |