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

Defining Interface as wall

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 29, 2009, 05:31
Default Defining Interface as wall
  #1
New Member
 
Join Date: Mar 2009
Posts: 28
Rep Power: 7
Hectux is on a distinguished road
Hi,

I am using interFoam with two fluids in a channel.

is it possible to define an interface beetween the two fluids as a wall BC?
I want to avoid convection beetween the two phases.

Or can anyone tell me wehere to find the sourcecode regarding the interface calculation and/or calculcation of the BCs?


Best regards
Hectux is offline   Reply With Quote

Old   June 3, 2009, 03:39
Default
  #2
New Member
 
Join Date: Mar 2009
Posts: 28
Rep Power: 7
Hectux is on a distinguished road
I ran several cases with segmented flow types.

My problem is. Actually I am expecting internal circulations in the segments because the interface should be - as described in my literature - "a sharp boundary with no flow through it."

But my flow looks like the pictures attached.

I have chosen a realistic surface tension of 7.46 mN/m.
The wall is no-slip.
Solver: interFoam.

I do not know what's wrong. The convection looks like it goes through the interface.
Attached Images
File Type: jpg slug1.jpg (8.0 KB, 20 views)
File Type: jpg slug2.jpg (15.7 KB, 26 views)
File Type: jpg slug3.jpg (13.0 KB, 26 views)
Hectux is offline   Reply With Quote

Old   June 4, 2009, 12:09
Default
  #3
Senior Member
 
Join Date: Mar 2009
Posts: 248
Rep Power: 8
jaswi is on a distinguished road
Dear Adam

I have read your message as well. Here is all I can say...

Yes, I have worked with the interFoam but i do not understand what exactly you are trying to setup. As far as my understanding goes the interface between two fluids (in case of interfoam) is not something you can apply boundary conditions to. One can apply BC to an interface when it is explicitly defined and tracked. Also please excuse me of my stupid question but I do not see what do you mean by wall boundary condition on the interface.

When you define a two phases in interfoam , it is represented by the scalar field gamma. And if you want to apply some boundary condition on the interface then you have to explicitly look for the cells where your interface lies and then modify accordingly. That will require some implementation effort.

Once again please specify what is your physical setup ?
What do you mean by saying that there should be no convection between the fluids ?
Do you mean the convection of concentration of a species ?

Frankly speaking my friend, I am a bit confused by your problem setup and cannot say much in this case.

Best Regards
Jaswi
jaswi is offline   Reply With Quote

Old   June 5, 2009, 03:41
Question Need the relative velocity
  #4
New Member
 
Join Date: Mar 2009
Posts: 28
Rep Power: 7
Hectux is on a distinguished road
I think I know where the problem is.

The simulation is correct. There is no flux through the interface. The problem is:

The vectors show me the absolute velocity. As the slugs (segments) move in x-direction (from left to right) all vector arrows are rightist.
What I need to see is the relative velocity field.

But how the get the relative velocity of one phase (gamma=1) to e.g. the other face (gamma=0)?
Then I should see the velocity field in the phase gamma=1 without the overlaid translational movement to the right.
Hectux is offline   Reply With Quote

Old   June 5, 2009, 06:50
Default
  #5
Senior Member
 
Join Date: Mar 2009
Posts: 248
Rep Power: 8
jaswi is on a distinguished road
Dear Adam

Greetings

The interFoam uses VOF method i.e. the velocity has a single field representation. If you would like to have relative velocity at the interface then multiplying the velocity of cells at the interface with the volume fraction will give you the velocity of the heavier phase. For the start just mutliplying the complete velocity field volVectorField U with the complete volScalarField gamma and that will give you the velocity field for the heavier phase. May be this is what you need.

Hope that was what you were looking for
BR
Jaswi
jaswi is offline   Reply With Quote

Old   June 5, 2009, 09:26
Thumbs up Works good but.....
  #6
New Member
 
Join Date: Mar 2009
Posts: 28
Rep Power: 7
Hectux is on a distinguished road
Ok. Thats a start.

I am really thankful.

As I am a newbie in OpenFOAM I copied the some stuff from the vorticity utility and added the rest as I thought it was right.

My code snipplet of the important calculation is now - as you sad - like this:
---------------------------------------------

volVectorField relU
(
IOobject
(
"relU",
runTime.timeName(),
mesh,
IOobject::NO_READ
),
U * gamma
);

----------------------------------------------

This works. I now have only the vectors of my phase gamma=1 and the other vectors have the length 0 cos I multiplied it with gamma=0,

BUT! I still have one question.

U is a vector field and accourdingly U * gamma. I would now subtract from every vector in the new field relU (=U * gamma) a constant vector (e.g. (-1 0 0)).

How can I do this?

I cannot write something like this:
U * gamma - (1 0 0)
cos on the left side I got a field and on the right one a vector.
Hectux is offline   Reply With Quote

Old   June 8, 2009, 04:45
Lightbulb
  #7
New Member
 
Join Date: Mar 2009
Posts: 28
Rep Power: 7
Hectux is on a distinguished road
Got it.

I just created a Vector Field for the Mesh with setFields which expresses the translative movement of the segments.

Then I could substract the vector fields from each other and I only got the relative velocity field within the segments. Which is a recirculation
Hectux 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
Fluent3DMeshToFoam simvun OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 48 May 14, 2012 05:20
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug unoder OpenFOAM Installation 11 January 30, 2008 20:30
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 15:55
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12
Replace periodic by inlet-outlet pair lego CFX 3 November 5, 2002 20:09


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