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

Forces.C code explaination

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 29, 2017, 16:15
Default Forces.C code explaination
  #1
Member
 
Join Date: Jun 2017
Posts: 73
Rep Power: 8
Friendly is on a distinguished road
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]);
Does anyone know what is meant by Sfb and RhoReffb or Reffb? I think Sfb is the area of the face i. "dev" subtracts the trace from the tensor. RhoReffb should be something like mu*grad(u).
Friendly is offline   Reply With Quote

Old   November 29, 2017, 17:04
Default
  #2
Member
 
Join Date: Jun 2017
Posts: 73
Rep Power: 8
Friendly is on a distinguished road
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?
Friendly is offline   Reply With Quote

Old   November 30, 2017, 02:19
Default
  #3
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Pool Gor mathematics, numerics, derivations and openfoam. There you find the definition of DevRhoREff
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   November 30, 2017, 07:09
Default
  #4
Member
 
Join Date: Jun 2017
Posts: 73
Rep Power: 8
Friendly is on a distinguished road
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.
Friendly is offline   Reply With Quote

Old   November 30, 2017, 07:40
Default
  #5
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
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
Tobi is offline   Reply With Quote

Old   November 30, 2017, 11:35
Default
  #6
Member
 
Join Date: Jun 2017
Posts: 73
Rep Power: 8
Friendly is on a distinguished road
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]?
Friendly 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
The FOAM Documentation Project - SHUT-DOWN holger_marschall OpenFOAM 242 March 7, 2013 12:30
How to make code run in parallel? cwang5 OpenFOAM Programming & Development 1 May 30, 2011 04:47
Open Source Vs Commercial Software MechE OpenFOAM 28 May 16, 2011 11:02
Small 3-D code Zdravko Stojanovic Main CFD Forum 2 July 19, 2010 10:11


All times are GMT -4. The time now is 21:17.