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

vectorField Uwall(Up.patchInternalField() - Up)

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

Like Tree5Likes
  • 1 Post By boeleman
  • 1 Post By Lieven
  • 1 Post By boeleman
  • 2 Post By vigneshTG

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 25, 2013, 13:23
Default vectorField Uwall(Up.patchInternalField() - Up)
  #1
New Member
 
Join Date: May 2012
Posts: 17
Rep Power: 13
boeleman is on a distinguished road
Dear Foamers,

I am trying to understand the calculation of the dynamic contact angle in dynamicAlphaContactAngleFvPatchScalarField.C and I am having trouble to understand the line:

vectorField Uwall(Up.patchInternalField() - Up);

I know Up is U_.boundaryField()[patchi], but why subtract the velocities on patchi from their values next to patchi ("Up.patchInternalField()")?

Thanks,

Arnout
chery1986 likes this.
boeleman is offline   Reply With Quote

Old   January 25, 2013, 17:48
Default
  #2
Senior Member
 
Lieven
Join Date: Dec 2011
Location: Leuven, Belgium
Posts: 299
Rep Power: 22
Lieven will become famous soon enough
It seems to me that Uwall is the relative velocity with respect to the wall.
I expect that Up = 0 for fixed walls and nonzero for moving walls.

You can easily verify this by adding
Info << Up << endl;
to the code

Cheers,

L
vigneshTG likes this.
Lieven is offline   Reply With Quote

Old   January 27, 2013, 18:15
Default
  #3
New Member
 
Join Date: May 2012
Posts: 17
Rep Power: 13
boeleman is on a distinguished road
I had not thought of the possibility of moving walls. Now it makes sense. Thanks for your reply!
boeleman is offline   Reply With Quote

Old   July 10, 2014, 09:42
Default
  #4
Member
 
Vignesh
Join Date: Oct 2012
Location: Darmstadt, Germany
Posts: 66
Rep Power: 13
vigneshTG is on a distinguished road
Boeleman !!

I am also trying to understand the dynamic contact angle model implemented in the openfoam. I want to know whether "patch().nf()" refers to face unit normals or something else? Also, whats the difference between nf and nhat ?

Regards
Vignesh
vigneshTG is offline   Reply With Quote

Old   July 17, 2014, 16:36
Default
  #5
New Member
 
Join Date: May 2012
Posts: 17
Rep Power: 13
boeleman is on a distinguished road
nf(), does indeed refer to the face unit normal, so Uwall -= (nf & Uwall)*nf;calculates the velocity parallel to the wall, by subtracting the normal component from the complete velocity vector. nhat() is the normal to the interface of the two phases.
vigneshTG likes this.
boeleman is offline   Reply With Quote

Old   August 15, 2014, 16:11
Default
  #6
New Member
 
Anjishnu Choudhury
Join Date: Jul 2014
Posts: 4
Rep Power: 11
Anjishnu Choudhury is on a distinguished road
Hey,
I am really not getting the physical meaning of the term vectorField Uwall(Up.patchInternalField() - Up);

Uwall seems to be of type vectorField, but what does the term in brackets mean ?

Thanks in advance!
Anjishnu Choudhury is offline   Reply With Quote

Old   August 15, 2014, 19:02
Default
  #7
Member
 
Vignesh
Join Date: Oct 2012
Location: Darmstadt, Germany
Posts: 66
Rep Power: 13
vigneshTG is on a distinguished road
Quote:
Originally Posted by Anjishnu Choudhury View Post
Hey,
I am really not getting the physical meaning of the term vectorField Uwall(Up.patchInternalField() - Up);

Uwall seems to be of type vectorField, but what does the term in brackets mean ?

Thanks in advance!
Up contains the values of U in the boundary and Up.patchInternalField is the values of U in the cells next to the boundary. Their difference gives the relative velocity !!

Hope this helps
__________________
Thanks and Regards

Vignesh
vigneshTG is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[OpenFOAM] Visualisation of components of a vectorField nico ParaView 2 September 19, 2021 00:42
scalar points to vectorField deji OpenFOAM 2 January 26, 2012 13:28
FvPatchScalarField access operator maka OpenFOAM 3 July 1, 2008 09:32
Bug in vectorField hjasak OpenFOAM Running, Solving & CFD 1 September 7, 2007 09:09
Changes in OF 14 and compilation issues alberto OpenFOAM Running, Solving & CFD 13 May 13, 2007 16:21


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