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

Arglist

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 28, 2007, 12:15
Default Hi i am studing about the pr
  #1
New Member
 
morteza
Join Date: Mar 2009
Posts: 18
Rep Power: 8
morteza is on a distinguished road
Hi
i am studing about the programming in openfoam and i don't know the effect of this line at the potentialFoam.C solver:

argList::validOptions.insert("writep","")

Thanks

morteza
morteza is offline   Reply With Quote

Old   January 28, 2007, 18:58
Default Start potentialFoam with no ca
  #2
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,912
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Start potentialFoam with no case and no directory and you will be enlightend: in addition to complainig about the missing case/directory it will tell you that -writep is a perfectly valid option for that command. And why is it valid? Because the solver has been told so in the program line you mentioned.

What the option switches on/off, I don't know.
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   February 6, 2007, 09:27
Default Thank you Bernhard. I started
  #3
New Member
 
morteza
Join Date: Mar 2009
Posts: 18
Rep Power: 8
morteza is on a distinguished road
Thank you Bernhard.
I started potentialFoam with no case and directory and found what you say.

Thanks
Morteza
morteza is offline   Reply With Quote

Old   February 6, 2007, 09:40
Default Hi bernard Do you know how po
  #4
New Member
 
morteza
Join Date: Mar 2009
Posts: 18
Rep Power: 8
morteza is on a distinguished road
Hi bernard
Do you know how potentialFoam solves these two equations:
div(U)=0
laplacian(p)=0
as mentioned at programmers guide.
at programmers guide we have this main line for solving these two equations:

laplacian([0 0 1 0 0 0 ],p)=div(phi)
where phi is U.A.

Thanks

morteza
morteza is offline   Reply With Quote

Old   March 30, 2007, 06:10
Default Hello morteza! This might b
  #5
Member
 
Erik Arlemark
Join Date: Mar 2009
Location: Eindhoven, Netherlands
Posts: 47
Rep Power: 8
Erik is on a distinguished road
Hello morteza!

This might be a bit late but its a discusion on this on the following thread:

http://www.cfd-online.com/cgi-bin/Op...how.cgi?1/4014

/Erik
Erik is offline   Reply With Quote

Old   March 30, 2007, 06:13
Default Sorry the link didnt lead wher
  #6
Member
 
Erik Arlemark
Join Date: Mar 2009
Location: Eindhoven, Netherlands
Posts: 47
Rep Power: 8
Erik is on a distinguished road
Sorry the link didnt lead where I thought it would.

Here is a thread I saved from the forum:

Continuity Eq: div(U) = 0
and with source term: div(U) = S

Discretised momentum eq. A*U = H - grad(p)
->
U = H/A -(1/A)*grad(p) = Ustar - (1/A)*grad(p)
H/A is the momentum-predictor value of U, ie Ustar.
if we insert this in the continuity eq. it yields

div(U) = div(Ustar - (1/A)*grad(p)) = S
->
div(Ustar) - laplacian(1/A, p) = S
or
laplacian(1/A, p) = -S + div(phi)
where phi is Ustar evaluated on faces.
have a look at icoFoam and you will easily see where to insert S.

/Erik
Erik 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



All times are GMT -4. The time now is 02:53.