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

About interFoam solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 9, 2005, 05:39
Default Hi! I am trying to put up a c
  #61
bjornchr
Guest
 
Posts: n/a
Hi!
I am trying to put up a complete second order setup for a dambreak problem in interFoam.

Does anybody know a good choice of schemes?

I also wonder what backward means. It says that it is second order, is it better than the Crank Nicholson?

I hope for answer.

Bjorn
  Reply With Quote

Old   December 9, 2005, 06:40
Default Backward is slightly more diss
  #62
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 724
Rep Power: 11
eugene is on a distinguished road
Backward is slightly more dissipative than Crank Nicholson, but on the up side it also tends to be slightly more stable.
eugene is offline   Reply With Quote

Old   December 10, 2005, 08:30
Default Thanks a lot Eugene! I am s
  #63
bjornchr
Guest
 
Posts: n/a
Thanks a lot Eugene!

I am setting up a dambreak problem and I want to disregard viscousity. Is it as simple as setting mu to zero in the transportproperties file?

Best regards
Bjørn
  Reply With Quote

Old   April 24, 2007, 02:35
Default Hallo, I want to couple the i
  #64
New Member
 
Karl-Heinz Leitz
Join Date: Mar 2009
Posts: 16
Rep Power: 7
khleitz is on a distinguished road
Hallo,
I want to couple the interFoam-Solver with the laplacianFoam-Solver in order to calculate the heat conduction in a liquid.
I have changed the damBreak example so that I have a liquid in a pot and the heat source is at the top.
However I can't manage to give the air phase and the liquid phase different transport coefficients DT.
Can anybody help me.
Best regards,
Karl-Heinz
khleitz is offline   Reply With Quote

Old   April 24, 2007, 02:42
Default Hallo, I want to simulate a l
  #65
New Member
 
Karl-Heinz Leitz
Join Date: Mar 2009
Posts: 16
Rep Power: 7
khleitz is on a distinguished road
Hallo,
I want to simulate a local pressure on a liquid surface. Therefore I have modified the damBreak example so that I have a liquid in a pot.
I want to simulate a local pressure/force on the liquid surface in order to get a little dent.
Can anybody help me?
Best regards,
Karl-Heinz
khleitz is offline   Reply With Quote

Old   April 24, 2007, 04:24
Default Hi KH! About the varying DT
  #66
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,834
Rep Power: 38
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi KH!

About the varying DT: this is relatively easy. Just replace the scalar DT with a volScalarField. The value of that field is calculated every time before the calculation of the temperature equation:

DT=DTfluid*gamma+DTgas*(1-gamma);

(Or something similar)
__________________
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   November 7, 2007, 03:07
Default Hallo, can anybody explain m
  #67
New Member
 
Karl-Heinz Leitz
Join Date: Mar 2009
Posts: 16
Rep Power: 7
khleitz is on a distinguished road
Hallo,
can anybody explain me the pEqn in the interFoam-Solver?
What is rUAf?
Why is there a phi and a phiU? What's the difference?
What does adjustPhi(phi, U, pd) do?
What is the sense of U += rUA*fvc::reconstruct((phi - phiU)/rUAf)?
Best regards,
KH

{
volScalarField rUA = 1.0/UEqn.A();
surfaceScalarField rUAf = fvc::interpolate(rUA);

U = rUA*UEqn.H();

surfaceScalarField phiU
(
"phiU",
(fvc::interpolate(U) & mesh.Sf()) + fvc::ddtPhiCorr(rUA, rho, U, phi)
);

phi = phiU +
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(g amma)
- ghf*fvc::snGrad(rho)
)*rUAf*mesh.magSf();

adjustPhi(phi, U, pd);

for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
{
fvScalarMatrix pdEqn
(
fvm::laplacian(rUAf, pd) == fvc::div(phi)
);

pdEqn.setReference(pdRefCell, pdRefValue);

if (corr == nCorr-1 && nonOrth == nNonOrthCorr)
{
pdEqn.solve(mesh.solver(pd.name() + "Final"));
}
else
{
pdEqn.solve(mesh.solver(pd.name()));
}

if (nonOrth == nNonOrthCorr)
{
phi -= pdEqn.flux();
}
}

U += rUA*fvc::reconstruct((phi - phiU)/rUAf);
U.correctBoundaryConditions();
}
khleitz is offline   Reply With Quote

Old   November 7, 2007, 06:38
Default I wonder what you will make fr
  #68
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,741
Rep Power: 20
hjasak will become famous soon enough
I wonder what you will make from the answers...

Q: What is rUAf?

It is the inverse of momentum matrix diagonal interpolated to the face

Q: Why is there a phi and a phiU? What's the difference?

Because phi is a volume flux and phiU is only its convection-diffusion component. The rest of the correction, which appears on the rhs of the pressure equation is added after the pressure solution to allow face-based treatment of the buoyancy and surface tension terms.

Q: What does adjustPhi(phi, U, pd) do?

A: In cases where the pressure boundary condition is zero gradient all the way around, it adjust the outlet flow to satisfy global continuity.

Q: What is the sense of U += rUA*fvc::reconstruct((phi - phiU)/rUAf)?

That's the rest of the momentum terms, which were on the rhs of the pressure equation. At that point, U only held the convection-diffusion part - see comments above.

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   March 8, 2008, 07:17
Default Hallo, can anybody explain me
  #69
New Member
 
Karl-Heinz Leitz
Join Date: Mar 2009
Posts: 16
Rep Power: 7
khleitz is on a distinguished road
Hallo,
can anybody explain me how the continuity equation is implemented in the interfoam solver?
I guess it is included in the gammaEqn.
Can anybody explain this equation to me, I mean the different terms and their connection to the mass conservation?
Best regards,
Karl-Heinz
khleitz is offline   Reply With Quote

Old   March 8, 2008, 07:55
Default Try a PhD Thesis form Henrik R
  #70
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,741
Rep Power: 20
hjasak will become famous soon enough
Try a PhD Thesis form Henrik Rusche - all explained there.

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   August 18, 2008, 10:52
Default hi all, i have a transport(co
  #71
Member
 
davey david
Join Date: Mar 2009
Posts: 54
Rep Power: 7
suredross is on a distinguished road
hi all,
i have a transport(concentration) equation in my solver and would like to set initial conditions for half my mesh.i have cyclic,fixedwalls and electrodes as my boundaries.there is no inlet or outlet.anyone with ideas??
thanks in advance.
davey
suredross is offline   Reply With Quote

Old   August 18, 2008, 16:06
Default Hi Davey! "Initial conditio
  #72
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,834
Rep Power: 38
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi Davey!

"Initial conditions for half of my mesh": I assume you want to set Initial conditions in the WHOLE mesh. Just one value in one half and another value in the other half. There is a utility (setFields) that comes with OpenFOAM that helps you do that. If you're going for more complicated stuff, there is another utility floating around (funkySetFields is the name, I think). But if setFields is sufficient for your purposes, then I recommend that.

Bernhard
__________________
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   August 19, 2008, 09:04
Default Hi Bernhard, yes,its about se
  #73
Member
 
davey david
Join Date: Mar 2009
Posts: 54
Rep Power: 7
suredross is on a distinguished road
Hi Bernhard,
yes,its about setting a different value in each half(concentration)and i have had a look at setfields.unfortunately its not working for me.since it uses boxToCell and i want the value for half the mesh,i am having some problems.any idea of the input method for celltocell??
thanks in advance.

davey
suredross is offline   Reply With Quote

Old   August 19, 2008, 11:19
Default Just make the box big enough t
  #74
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,834
Rep Power: 38
gschaider will become famous soon enoughgschaider will become famous soon enough
Just make the box big enough to cover half the mesh and everything should be OK

Bernhard
__________________
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   August 19, 2008, 11:29
Default another thing with the dimensi
  #75
Member
 
davey david
Join Date: Mar 2009
Posts: 54
Rep Power: 7
suredross is on a distinguished road
another thing with the dimensions of the box,its labelled(x,y,z),right??with the coordinates representing a point in the mesh??i am quite hazy about how its done in the tutorial.also,my solver is a modified icofoam,do i get to assign the other physical properties(rho)and nu,again??i need to realise ,with time, the mixing as well.
sorry for my rash thougths-guess am having a long day...
cheers
davey
suredross is offline   Reply With Quote

Old   August 20, 2008, 06:09
Default Hi, figured it out now....run
  #76
Member
 
davey david
Join Date: Mar 2009
Posts: 54
Rep Power: 7
suredross is on a distinguished road
Hi,
figured it out now....running on my case.
thanks

davey
suredross is offline   Reply With Quote

Old   September 17, 2008, 07:01
Default Is there a way to plot the hor
  #77
Senior Member
 
musaddeque hossein
Join Date: Mar 2009
Posts: 251
Rep Power: 8
musahossein is on a distinguished road
Is there a way to plot the horizontal displacement in the dam break problem using Paraview? Would it be the same as plotting Ux using the Plot over line in the Filters Menu?

Thanks

Musa
musahossein is offline   Reply With Quote

Old   September 19, 2008, 04:28
Default (1) Pick up the patch in the R
  #78
Member
 
a a saha
Join Date: Mar 2009
Posts: 67
Rep Power: 7
asaha is on a distinguished road
(1) Pick up the patch in the Region Status and select gamma in Vol Field Status.
(2) Define a probeline from Filters menu.
(3) Select Contours gamma=0.5 from the Filters menu.
(4) Integrate Variables from Filters menu. Keep the contour flag visible by making the data visible.
(5) Then Plot Selection over time should give you the desired displacement with time.

I hope this helps. If you come to know a much simpler method please let me also know.
asaha is offline   Reply With Quote

Old   September 19, 2008, 13:59
Default Can the interFoam solver handl
  #79
Senior Member
 
chegdan's Avatar
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 508
Rep Power: 17
chegdan will become famous soon enough
Can the interFoam solver handle cyclic boundaries? I have borrowed some of the code from channelOodles to allow channel flow and extract gradP, but it seems that gamma is not passed through the cyclic boundaries. Does anyone have any insight on what else needs to be changed to pass gamma through the cyclic boundary? Any help is much appreciated.

Dan
__________________
Dan

Find me on twitter @dancombest
chegdan is offline   Reply With Quote

Old   September 19, 2008, 16:24
Default Saha: Thankyou very much for
  #80
Senior Member
 
musaddeque hossein
Join Date: Mar 2009
Posts: 251
Rep Power: 8
musahossein is on a distinguished road
Saha:
Thankyou very much for takning the time to write these instructions. I was tried it, but did not get any results. As a result I have a few questions. They are appended to your instruction as appropriate:

(1) Pick up the patch in the Region Status and select gamma in Vol Field Status.

Musa: I did that. That will leave the "internal mesh" only, correct?

(2) Define a probeline from Filters menu.

Musa: I did not see a probeline option. I saw a probe option in the Filters menu. I decided to put the probe at the intersection of the right side of the wall of water and the bottom of the tank. In other words, I specified the point as (1,0,0) assuming the format is x,y,z.

(3) Select Contours gamma=0.5 from the Filters menu.

Musa: I selected Contours from the Filters menu. However I did not see a option to specify gamma=0.5. Did you mean set the contour line interval at 0.5? Please let me know.

(4) Integrate Variables from Filters menu. Keep the contour flag visible by making the data visible.

Musa: Upon selecting Integrate Variables, I also selected "show data" option in the display tab. However, there is a pull down menu that permits selection of "point data", "cell data" and "field data". Do I need to select any of these or just stick with Point data? I hope you will be able to clarify.

(5) Then Plot Selection over time should give you the desired displacement with time.

Musa: I did, but nothing happened.

Any further suggestion/instructions will be appreciated.

Thanks!
musahossein 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
InterFoam MULES solver jaswi OpenFOAM Running, Solving & CFD 4 November 21, 2012 08:56
Wmake problem interFoam solver feijooos OpenFOAM Running, Solving & CFD 4 December 8, 2008 11:01
DICPCG solver in interFoam m9819348 OpenFOAM Running, Solving & CFD 1 September 20, 2007 13:10
About interfoam solver qiu OpenFOAM Running, Solving & CFD 0 May 6, 2007 22:48
Need documentation for interFOAM solver mer OpenFOAM Running, Solving & CFD 5 May 31, 2006 12:22


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