|
[Sponsors] |
November 29, 2017, 17:15 |
Forces.C code explaination
|
#1 |
Member
Join Date: Jun 2017
Posts: 73
Rep Power: 9 |
Hi,
I have two questiosn belonging the following code: Code:
vectorField fN ( rho(p)*Sfb[patchi]*(p.boundaryField()[patchi] - pRef) ); vectorField fT(Sfb[patchi] & devRhoReffb[patchi]); |
|
November 29, 2017, 18:04 |
|
#2 |
Member
Join Date: Jun 2017
Posts: 73
Rep Power: 9 |
Ok I found out that R() represents the reynolds stress tensor.
Reff() is the effective stress tensor, does anyone know what is meant by effective here? Any why there is a "b" at the end of Sf "b" and RhoReff "b"()? Boundary maybe? |
|
November 30, 2017, 08:09 |
|
#4 |
Member
Join Date: Jun 2017
Posts: 73
Rep Power: 9 |
Thanks, that helped a lot.
Just two more question: Why is there a "b" at the end? In "devRhoReffb" Rho is stands for a density based equation. I am calculationg an incompressible case and still getting values of it. Does this equation represent a general Approach to calculate fT? Maybe something changed, I am using OpenFOAM 4.1. |
|
November 30, 2017, 08:40 |
|
#5 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52 |
As we are only interested in the force acting on a patch, the b postfix just shows one that it is related to boundary values, which can be already seen in the snipped you provided. You can name the variables as you want.
__________________
Keep foaming, Tobias Holzmann |
|
November 30, 2017, 12:35 |
|
#6 |
Member
Join Date: Jun 2017
Posts: 73
Rep Power: 9 |
Ok thanks I thought that.
Btw I have found out that for the incompressible case devRhoReff equals devReff. Ok so the last part that I am missing is the missing normal vector. There should be a normal vector for calculating the pressure force and another one for the friction force. Are they in Sfb[patchi]? |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
The FOAM Documentation Project - SHUT-DOWN | holger_marschall | OpenFOAM | 242 | March 7, 2013 13:30 |
How to make code run in parallel? | cwang5 | OpenFOAM Programming & Development | 1 | May 30, 2011 05:47 |
Open Source Vs Commercial Software | MechE | OpenFOAM | 28 | May 16, 2011 12:02 |
Small 3-D code | Zdravko Stojanovic | Main CFD Forum | 2 | July 19, 2010 11:11 |