CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

fvm::div(phi, Ug)...how to insert properly?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 17, 2012, 05:11
Default (U_j+Ug_j)*d(U_i+Ug_i)/dx_j...how to insert properly?
  #1
Member
 
Tom
Join Date: Jan 2012
Location: France
Posts: 43
Rep Power: 4
sawyer86 is on a distinguished road
Hi there,

My question can be seen even from this prospective. I have to insert in the equation this term:

(U_j+Ug_j)*d(U_i+Ug_i)/dx_j

where Ug_i are the components of a given velocity field and U is the velocity I want to solve. This above is the convective term in the expression "u*grad(u)"-like

how can I do? Is it possible?

Thank you in advance!

Last edited by sawyer86; October 17, 2012 at 13:08.
sawyer86 is offline   Reply With Quote

Old   October 17, 2012, 05:23
Default
  #2
Member
 
Join Date: Mar 2012
Location: Munich, Germany
Posts: 67
Rep Power: 4
treima is on a distinguished road
Hi,

how did you initialise Ug, is it equal to the initialisation of U? Or perhaps there are some problems with the implicit and explicit treatment. You can have a look at this thread: Incompatible fields for operation. There is a description of the same problem.


regards

treima
treima is offline   Reply With Quote

Old   October 17, 2012, 05:24
Default
  #3
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 912
Rep Power: 16
akidess will become famous soon enough
You do not solve for Ug, hence that term is explicit. phi is not the variable you solve for, it's the data from the previous iteration.
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Check out the scientific computing exchange http://scicomp.stackexchange.com
akidess is offline   Reply With Quote

Old   October 17, 2012, 05:32
Default (U_j+Ug_j)*d(U_i+Ug_i)/dx_j
  #4
Member
 
Tom
Join Date: Jan 2012
Location: France
Posts: 43
Rep Power: 4
sawyer86 is on a distinguished road
Quote:
Originally Posted by akidess View Post
You do not solve for Ug, hence that term is explicit. phi is not the variable you solve for, it's the data from the previous iteration.
Hi,

one question. Do you say that the volume velocity flux is given from the previous iteration...that's so the way how icoFoam treats the non linear-convective term in the way to linearize it, am I right?

This linearization influence the calculation of UEqn.H()...but if I do insert that term explict how does UEqn.H() must change?

My question can be seen even from this prospective. I have to insert in the equation this term:

(U_j+Ug_j)*d(U_i+Ug_i)/dx_j

how can I do? Is it possible?

Thank you very much
sawyer86 is offline   Reply With Quote

Old   October 17, 2012, 05:36
Default
  #5
Member
 
Tom
Join Date: Jan 2012
Location: France
Posts: 43
Rep Power: 4
sawyer86 is on a distinguished road
Quote:
Originally Posted by akidess View Post
You do not solve for Ug, hence that term is explicit. phi is not the variable you solve for, it's the data from the previous iteration.
(U_j+Ug_j)*d(U_i+Ug_i)/dx_j

Sorry I did not explain this term. It is essentially the convective term when you do a kind of linearization starting from a base flow...like in the RANS, but different...

Then the term i exactly want to add is the part of the above expression:

U_j*dUg_i/dx_j

...?

Last edited by sawyer86; October 17, 2012 at 06:24.
sawyer86 is offline   Reply With Quote

Old   October 17, 2012, 05:45
Default
  #6
Member
 
Tom
Join Date: Jan 2012
Location: France
Posts: 43
Rep Power: 4
sawyer86 is on a distinguished road
Quote:
Originally Posted by treima View Post
Hi,

how did you initialise Ug, is it equal to the initialisation of U? Or perhaps there are some problems with the implicit and explicit treatment. You can have a look at this thread: Incompatible fields for operation. There is a description of the same problem.


regards

treima
Thank you for the answer but i had a look of that thread and this is not just what I am looking for. By the way thank you again!

Regards
sawyer86 is offline   Reply With Quote

Reply

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
insert picture carno Main CFD Forum 4 May 11, 2011 07:33
Meshing a 2D airfoil properly jms Main CFD Forum 2 March 28, 2011 03:30
GAMBIT Error Msg database file not properly closed Janomano ANSYS Meshing & Geometry 8 November 23, 2010 04:14
Insert a velocity waveform Hila FLUENT 0 April 21, 2006 07:04
3D Insert modelling in Mold Flow Dan Rößler Main CFD Forum 0 December 20, 2005 06:37


All times are GMT -4. The time now is 23:22.