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

IcoFoam query

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 23, 2008, 22:19
Default Hello I find the following
  #1
Member
 
srinath
Join Date: Mar 2009
Location: Champaign, USA
Posts: 90
Rep Power: 8
srinath is on a distinguished road
Hello

I find the following quite puzzling.
In icoFoam.C, i find the foll. lines

solve(UEqn == -fvc::grad(p));

// --- PISO loop

for (int corr=0; corr<nCorr; corr++)
{
volScalarField rUA = 1.0/UEqn.A();

U = rUA*UEqn.H();
......
}

Solving UEqn presumably changes the value of U. But what is the point if couple of lines later, we store rUA*Ueqn.H() in U

If solving UEqn is a predictor for a U value, what is the point if we never use the value and overwrite it?

I would appreciate it if someone could clarify this as i find it quite puzzling.

Regards
Srinath
p.s. commenting out the solve UEqn line, seems to give visually acceptable results for the cavity case, but diff says the results are different. So clearly it is not a bug, but is doing something quite subtle
srinath is offline   Reply With Quote

Old   December 5, 2008, 05:00
Default Hello, Thank you for your qu
  #2
Member
 
Fabian Peng Karrholm
Join Date: Mar 2009
Posts: 61
Rep Power: 8
fabianpk is on a distinguished road
Hello,
Thank you for your question. This gives me an excellent opportunity to advertise my thesis. This exact question is covered in Appendix A.
The thesis main topic is diesel spray simulations using OpenFOAM, as well as simulation of cavitating flow in a model injector. However, since I have pondered on questions similar to the one you put out during my PhD studies, I put it as an appendix for myself to read, should I ever forget what happens in icoFoam.
A free version is available at:

http://powerlab.fsb.hr/ped/kturbo/Op...olmPhD2008.pdf

If you want a paper copy, contact me by email.

/Fabian
fabianpk is offline   Reply With Quote

Old   December 6, 2008, 10:01
Default Thanks Fabian Appendix A in
  #3
Member
 
srinath
Join Date: Mar 2009
Location: Champaign, USA
Posts: 90
Rep Power: 8
srinath is on a distinguished road
Thanks Fabian

Appendix A in your thesis was very useful.

Am also reading the rest of it as it is very interesting.

Cheers
Srinath
srinath is offline   Reply With Quote

Old   December 6, 2008, 11:49
Default Hello, Fabian My name is Sur
  #4
Senior Member
 
Suresh kumar Kannan
Join Date: Mar 2009
Location: Luxembourg, Luxembourg, Luxembourg
Posts: 129
Rep Power: 8
kumar is on a distinguished road
Hello, Fabian
My name is Suresh. I have just stārted using OpenFOAM a few weeks back. I saw your thesis, your work seems to be quiet interesting. I would like to let you know that My P.hD research work is also similar to your topic. I amalso using OpenFOAM to study the break up of a liquid sheet from a nozzle. But a major part of my research work will include Nonlinear instability analysis of the liquid sheet.
I am using lesInterFOAM of OpenFOAM v1.5 to do these simulations and I have already performed one simulation on a simple geometry.
Since I have just started using OpenFOAM, i have a few basic questions and I hope you wouldnt mind answering them:
1) I found that the icoFOAM, interFOAM, lesInterfOAM are segregated solvers based on the PISO algorithm, just correct me if i am wrong in understanding .(with VOF for surface tracking and LES for turbulence modeling)

2) But I found in your thesis that you are including the compressibility effects, isnt it advicable to solve the equations in a coupled manner like in a coupled solver to solve a compressible fluid problem.
3) I also have a very basic question, is there any coupled solver in OpenFOAM?
I have some experience in developing a coupled solver, so i understand the structure of a coupled solver. I also have some idea about the basic methodology of a segeregated solver as explained in basic textbooks by S.V.Patankar and Malalkasera.
4) For my simulations in OpenFOAM I have a setup like injecting 200bar pressure fluid into the chamber with 10bar pressure. This is just a test case given to me by Delphi.

bye
thanks for your reply inadvance
with regards
K.Suresh kumar
kumar is offline   Reply With Quote

Old   December 9, 2008, 03:49
Default Hello, 1) Yes, to find an e
  #5
Member
 
Fabian Peng Karrholm
Join Date: Mar 2009
Posts: 61
Rep Power: 8
fabianpk is on a distinguished road
Hello,

1) Yes, to find an example of a solver that's not PISO, look at simpleFoam.
2) Perhaps it might be suitable, and there are quite a few coupled solver codes out there now, but PISO works fine for compressible flow (in my experience).
3) Not as far as I know, no. I think to develop a coupled solver in OpenFOAM you would have to dig into the source code much further than the top-level code.

/Fabian
fabianpk is offline   Reply With Quote

Old   December 9, 2008, 04:46
Default Hello Fabian, Th
  #6
Senior Member
 
Suresh kumar Kannan
Join Date: Mar 2009
Location: Luxembourg, Luxembourg, Luxembourg
Posts: 129
Rep Power: 8
kumar is on a distinguished road
Hello Fabian,
Thankyou very much for your reply. I will have a look at simpleFoam.
But for my problem I dont have to worry about the coupled solver as my problem corresponds to an incompressible fluid.
I was just curious to know more about the solver.
I just have one more question can you provide me some reference on the algorithm and some details about the lesInterFOAM solver, if you have any?
thankyou very much for reply
bye
K.Suresh kumar
kumar 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
IcoFoam aap OpenFOAM Running, Solving & CFD 15 May 28, 2012 08:30
Density in icoFoam Densidad en icoFoam manuel OpenFOAM Running, Solving & CFD 8 September 22, 2010 04:10
About phi in icoFoam kar OpenFOAM Running, Solving & CFD 3 February 20, 2008 06:20
Possible bug in icoFoam msrinath80 OpenFOAM Bugs 6 November 19, 2007 18:35
IcoFoam on AIX 53 ds2taieb OpenFOAM Installation 1 March 24, 2006 04:22


All times are GMT -4. The time now is 13:15.