CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

gauss grad no correction for coupled boundary??!!

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

LinkBack Thread Tools Display Modes
Old   December 19, 2011, 19:19
Default gauss grad no correction for coupled boundary??!!
Senior Member
Join Date: Nov 2009
Location: Michigan
Posts: 135
Rep Power: 10
doubtsincfd is on a distinguished road
In the code for gauss gradient (gaussGrad.c), the grad calcuation at the boundary is done as:

forAll(vsf.boundaryField(), patchi)
if (!vsf.boundaryField()[patchi].coupled())
const vectorField n
/ vsf.mesh().magSf().boundaryField()[patchi]

gGrad.boundaryField()[patchi] += n *
- (n & gGrad.boundaryField()[patchi])

So is there no correction for coupled patches?? I am getting embarrassingly wrong values of gradient at cyclic patch.

Second question: In the same code, calculation of grad has been done for three sets of mesh entities:
1. On face centers on patch
2. On cell centers adjacent to patch
3. On all inside cell centers

Is my understanding correct?
doubtsincfd is offline   Reply With Quote


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
Problem in fvschemes divSchemes cannot use Gauss linearUpwind nico765 OpenFOAM Bugs 27 June 3, 2014 05:29
please help me which solver is better for my application Ger_US OpenFOAM 8 January 8, 2013 13:29
natural convection mehrdadeng CFX 10 February 25, 2011 06:25
solution diverges when linear upwind interpolation scheme is used subash OpenFOAM 0 May 29, 2010 01:23
Stability problems with kepsilon in external aero edreed OpenFOAM Running, Solving & CFD 21 July 16, 2008 15:00

All times are GMT -4. The time now is 09:54.