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

VOF method

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

Like Tree9Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   October 12, 2010, 07:57
Default paper request
  #61
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 9
ata is on a distinguished road
Hi Dear FOAMers
May any one tell me how can I get the paper
"Drop impact onto a liquid layer of finite thickness: Dynamics of the cavity evolution"?
Best regards
ata is offline   Reply With Quote

Old   October 12, 2010, 08:05
Default
  #62
Senior Member
 
Suresh kumar Kannan
Join Date: Mar 2009
Location: Luxembourg, Luxembourg, Luxembourg
Posts: 129
Rep Power: 9
kumar is on a distinguished road
Hi ata give me your email I.D . I can send it to you.

regards
K.Suresh kumar
kumar is offline   Reply With Quote

Old   October 12, 2010, 08:15
Default
  #63
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 9
ata is on a distinguished road
Hi Kumar
I sent you via a private message in the forum.
Thank you very much
ata is offline   Reply With Quote

Old   February 4, 2011, 13:56
Default
  #64
New Member
 
Alexey Ginevsky
Join Date: Oct 2010
Location: Moscow
Posts: 6
Rep Power: 7
algin is on a distinguished road
Hi ata,
Could you send me paper "Drop impact onto a liquid layer of finite thickness: Dynamics of the cavity evolution"?

My email: alexeyginevsky@gmail.com
Thank you in advance.
algin is offline   Reply With Quote

Old   February 5, 2011, 00:05
Default Drop impact onto a liquid layer of finite thickness Dynamics of the cavity evolution
  #65
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 9
ata is on a distinguished road
Hi
I sent it.
Good luck
ata is offline   Reply With Quote

Old   April 8, 2011, 04:26
Default
  #66
New Member
 
Klaus Schnitzlein
Join Date: Mar 2009
Posts: 7
Rep Power: 9
schnitzlein is on a distinguished road
Hi ata,
please could you send me your paper "Drop impact onto a liquid layer of finite thickness: Dynamics of the cavity evolution", too?

My email: klaus@rt.tu-cottbus.de
Thanks in advance.

linoth
schnitzlein is offline   Reply With Quote

Old   April 8, 2011, 10:39
Default
  #67
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 9
ata is on a distinguished road
Hi
I sent it.
Good luck
ata is offline   Reply With Quote

Old   September 8, 2011, 02:35
Default
  #68
New Member
 
NieYongguang
Join Date: Sep 2010
Posts: 27
Rep Power: 7
nygbook is on a distinguished road
Dear foamer
I have review the code of interFoam and interPhaseChangeFoam. However, I can't understand the function usage of MULES. It write in interPhaseChangeFoam: MULES::explicitSolve(oneField(), alpha1, phi, phiAlpha, Sp, Su, 1, 0). What mean does 'Su' in the code. Anybody can me a reference about MULES method?
nygbook is offline   Reply With Quote

Old   September 8, 2011, 03:39
Default
  #69
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 1,112
Rep Power: 20
akidess will become famous soon enough
Sp and Su are implicit and explicit source terms (e.g. for cavitation simulations).
akidess is offline   Reply With Quote

Old   September 8, 2011, 22:49
Default
  #70
New Member
 
NieYongguang
Join Date: Sep 2010
Posts: 27
Rep Power: 7
nygbook is on a distinguished road
Quote:
Originally Posted by akidess View Post
Sp and Su are implicit and explicit source terms (e.g. for cavitation simulations).
Thank you for your reply! However, I still have a question. The code of pEqn.H in interPhaseChangeFoam is below
Code:
fvScalarMatrix pEqn
        (
            fvc::div(phi) - fvm::laplacian(rUAf, p)
          - (vDotvP - vDotcP)*pSat + fvm::Sp(vDotvP - vDotcP, p)
        );
why isn't fvm::Su((vDotvP-vDotcP)*pSat, p)?
nygbook is offline   Reply With Quote

Old   September 9, 2011, 02:11
Default
  #71
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 1,112
Rep Power: 20
akidess will become famous soon enough
Because there is no fvm::Su() function. Explicit source terms are used without any additional notation.
akidess is offline   Reply With Quote

Old   September 9, 2011, 02:38
Default
  #72
New Member
 
NieYongguang
Join Date: Sep 2010
Posts: 27
Rep Power: 7
nygbook is on a distinguished road
Quote:
Originally Posted by akidess View Post
Because there is no fvm::Su() function. Explicit source terms are used without any additional notation.
no! no! I think it's wrong. The code of fvm::Su() can be found in fvmSup.c
Code:
00032 template<class Type>
00033 Foam::tmp<Foam::fvMatrix<Type> >
00034 Foam::fvm::Su
00035 (
00036     const DimensionedField<Type, volMesh>& su,
00037     const GeometricField<Type, fvPatchField, volMesh>& vf
00038 )
00039 {
00040     const fvMesh& mesh = vf.mesh();
00041 
00042     tmp<fvMatrix<Type> > tfvm
00043     (
00044         new fvMatrix<Type>
00045         (
00046             vf,
00047             dimVol*su.dimensions()
00048         )
00049     );
00050     fvMatrix<Type>& fvm = tfvm();
00051 
00052     fvm.source() -= mesh.V()*su.field();
00053 
00054     return tfvm;
00055 }
It is evident that fvm::Su mean explicit source term multiply by the volume of mesh. The code of fvm::Sp can also be found in fvmSup.C that
Code:
00098 template<class Type>
00099 Foam::tmp<Foam::fvMatrix<Type> >
00100 Foam::fvm::Sp
00101 (
00102     const DimensionedField<scalar, volMesh>& sp,
00103     const GeometricField<Type, fvPatchField, volMesh>& vf
00104 )
00105 {
00106     const fvMesh& mesh = vf.mesh();
00107 
00108     tmp<fvMatrix<Type> > tfvm
00109     (
00110         new fvMatrix<Type>
00111         (
00112             vf,
00113             dimVol*sp.dimensions()*vf.dimensions()
00114         )
00115     );
00116     fvMatrix<Type>& fvm = tfvm();
00117 
00118     fvm.diag() += mesh.V()*sp.field();
00119 
00120     return tfvm;
00121 }
The relative theory can be found in p.87 of Jasak' thesis.
nygbook is offline   Reply With Quote

Old   September 9, 2011, 03:52
Default
  #73
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 1,112
Rep Power: 20
akidess will become famous soon enough
Mea culpa! I do still believe that writing "fvm::Su((vDotvP-vDotcP)*pSat, p)" is equivalent to "(vDotvP-vDotcP)*pSat*p". Anyway, I think the gist of your question is why the source term is treated implicitly? My guess is that (vDotvP-vDotcP) is usually larger zero, so when including the expression in the matrix instead of moving it to the right hand side of the equation system you will increase diagonal dominance.
akidess is offline   Reply With Quote

Old   September 9, 2011, 05:51
Default
  #74
New Member
 
NieYongguang
Join Date: Sep 2010
Posts: 27
Rep Power: 7
nygbook is on a distinguished road
Quote:
I do still believe that writing "fvm::Su((vDotvP-vDotcP)*pSat, p)" is equivalent to "(vDotvP-vDotcP)*pSat*p".
Yes, It should equal to (vDotvP-vDotcP)*pSat*p.mesh().v(). I mean if it doesn't use fvm::Su, the explicit source item isn't multiplied by mesh().v(). Is it right?
nygbook is offline   Reply With Quote

Old   March 25, 2012, 11:18
Default
  #75
Member
 
ehsan
Join Date: Mar 2009
Posts: 98
Rep Power: 9
ehsan is on a distinguished road
Dear all

I need a reference which describes by more details the new local-time stepping (LTS) approach applied in VOF in openfoam.

Thanks
ehsan is offline   Reply With Quote

Old   May 17, 2012, 23:20
Talking Summary + some questions
  #76
New Member
 
Pharg Mandadapu
Join Date: Jul 2011
Posts: 16
Rep Power: 6
pharg_yrartibra is on a distinguished road
I found this thread when I was looking for reference documentation for InterFoam, and decided to make like a short summary of it and at the same time ask some interesting questions that I think is still left unanswered.

Quote:
Originally Posted by Henry Weller (Henry) View Post
> i.e. is it CICSAM or a different one

It's a new technique I developed a few years ago to resolve some of the fundamental problems of CICSAM and other traditional VOF interface compression methods. The differences have already been debated at length, have a look through previous threads on the subject.
Can you please link to any of these threads? I don't know the name of your method so I don't really know what to search for. I have already tried searching for CICSAM and InterFoam, but the hits I get are no good.

Quote:
Originally Posted by Ali (Ali) View Post
I was wondering what's the actual and preferred expression for the compressive term in gamma equation and a very brief explanation of how it is derived. The formula for 'phir' in the code differs from what Henrik has mentioned in (Eq. (4.15)). In the code:

phir=cGamma()*mag(phi/mesh.magSf())*interface.nHatf()

while Henrik's suggests (if I have formulated it correctly):

phir=cGamma*nhatf*max(nhat * phi / mesh.magSf()**2 )
Hi Ali, where can I find that expression in the interFoam code? Could you refer me to the file and line number, please?

Quote:
Originally Posted by henry View Post
Sorry, I haven't written many papers; it's too time consuming and stops me doing the interesting work.
Well, you should! Your method seems very interesting and I'm sure there are a lot of people (including me of course) who would be very interested in reading how it works. Unless there already is a lot of information about it on the web, but as I said, I haven't been able to find any...

Quote:
Originally Posted by hjasak View Post
There is a decent description of the implemented Interface-Capturing Methodology (probably a bit out of date now) in a PhD Thesis by dr. Henrik Rusche:

Computational Fluid Dynamics of Dispersed Two-Phase Flows at High Phase Fractions, Henrik Rusche, Imperial College of Science, Technology & Medicine, December 2002.

Hrv
Quote:
Originally Posted by henry View Post
That describes an old version of my method, not the one currently in interFoam.
Okay.

Quote:
Originally Posted by sega View Post
I'm currently having a look into the gammaEqn.H from OF 1.4.1

Well, some parts look like the equation (4.15) from Rusches PhD thesis.

But still I'm not sure about some issues:

- Where is the time derivative of gamma?
- What is interface.nHatf()? Is it the normal vector to the cell face?
- What doese MULES::explicitSolve01(gamma, phi, phiGamma) mean?

Thanks so far!
Quote:
Originally Posted by caw View Post
Dear Sebastian,

the time derivative is not there because of the usage of MULES. Look at OpenFOAM-1.5.x\src\finiteVolume\fvMatrices\solvers\MULES\MU LES.H for details:

" MULES: Multidimensional universal limiter with explicit solution
So I guess MULES is the name of the advection scheme used in InterFoam. Good to know!

Quote:
Originally Posted by sega View Post
Dear Edin.

Thanks for your answere.
So, the equation will take this form?
Hi Sebastian, did you get any luck in finding the correct expression for the equation used in InterFoam that describes the advection of gamma?

Quote:
Originally Posted by kumar View Post
Hello Edin Berberovic and santiago,
I have been following this post closely to understand the MULES and VOF methodology implemented in OpenFOAM.
I just got the paper "Drop impact onto a liquid layer of finite thickness: Dynamics of the cavity evolution". Thanks santiago for mentioning it, since i have been looking for a paper which explains the MULES.
I want to confirm that the MULES concept explained in that paper is the one implemented in OF-1.5, in the interFoam solver
What about OF-1.6.x, in which the pressure difference is handled in a different way. So could anybody tell me what is the major difference between the MULES explained in the paper and the one currently in OF-1.6.x, interFoam solver.

regards
K.Suresh kumar
Quote:
Originally Posted by santiagomarquezd View Post
Kumar, Edin's paper is one of the most comprehensive papers about interFoam that I have found. It's true, it explains interFoam version 1.5 but the part about MULES is valid in 1.6.x also.
Hi Kumar and Santiago, Edin's paper looks very interesting, thank you for sharing this information with us. So I guess this is like the documentation for MULES?

Quote:
Originally Posted by eberberovic View Post
Hi Isabel,

I've just sent you the paper PRE 2009. There is one typo in it, namely Eq. (24) represents the face flux rather than velocity at the cell-face. Sorry for this, it went also unnoticed by the reviewers.

Edin.
Okay, good to know. Once again, your paper seems very interesting (I'm about to read it soon!). I have a question about MULES, though, does it conserve mass? Can mass get 'clipped' away, or is it just the flux between cells that can get clipped?

Quote:
Originally Posted by ehsan View Post
Dear all

I need a reference which describes by more details the new local-time stepping (LTS) approach applied in VOF in openfoam.

Thanks
I didn't even know that local-time stepping existed, but I it seems like a good idea if a very small fraction of all cells need to be solved with a much smaller time step than the rest of the cells, to prevent the whole simulation to be significantly slowed down by a few "troublesome cells". All I could find about the method was what was written on this page. I would also really like to see some reference for how this method works... Maybe the LTS approach isn't a part of InterFoam, but it would be interesting to know more about it now that it has been brought up.

--Pharg
Pirlu likes this.

Last edited by pharg_yrartibra; May 17, 2012 at 23:48.
pharg_yrartibra is offline   Reply With Quote

Old   June 12, 2013, 23:24
Default
  #77
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 742
Rep Power: 9
sharonyue is on a distinguished road
Hi All,

I found all of you guys are working around VOF method, I am a newbie here, but I am facing a problem regarding the changeable viscosity in VOF. For example: if the liquid is kind of non-newtonian fluid. More details I wrote it here, interFoam, wierd nu.

Can you guys help me check this out? Is this a problem only regarding with the viscosity model or with VOF? Thanks in advance.
sharonyue is offline   Reply With Quote

Old   November 5, 2014, 15:32
Default
  #78
Member
 
james wilson
Join Date: Aug 2014
Location: Orlando, Fl
Posts: 36
Rep Power: 3
jameswilson620 is on a distinguished road
@isabel

"
In the interFoam solver, what these lines in the file gamma.Eqn.H mean?

surfaceScalarField phiGamma =
fvc::flux
(
phi,
gamma,
gammaScheme
)
+ fvc::flux
(
-fvc::flux(-phir, scalar(1) - gamma, gammarScheme),
gamma,
gammarScheme
);
"

this refers to the compression term used in MULES. note (gamma) and (1-gamma). These filter the velocity field and gives us a way to evaluate Ur = U1 - U2 which is used in the interface compression term. Also look up fvcFlux.* in your install directory. this will shed some light on what phiGamma becomes.

see: http://infoscience.epfl.ch/record/130534

and

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

to see the interface compression term.
jameswilson620 is offline   Reply With Quote

Old   June 8, 2015, 08:05
Default
  #79
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 319
Blog Entries: 1
Rep Power: 9
Tushar@cfd is on a distinguished road
Quote:
Originally Posted by santiagomarquezd View Post
Kumar, Edin's paper is one of the most comprehensive papers about interFoam that I have found. It's true, it explains interFoam version 1.5 but the part about MULES is valid in 1.6.x also.
Changes in pressure treatment from pd to p are not affecting the MULES part, I've been studying MULES and it's important only in solving the non-linear advection equation (Eq. 15), in fact it is not mentioned in Edin's paper. Changes in pd affects only in mommentum equation formulation (Eq. 20).
Calling MULES is a strategy to solve Eq. 15, taking a look of alphaEqn.H we have,

Code:
00009     for (int aCorr=0; aCorr<nAlphaCorr; aCorr++)
00010     {
00011         surfaceScalarField phiAlpha =
00012             fvc::flux
00013             (
00014                 phi,
00015                 alpha1,
00016                 alphaScheme
00017             )
00018           + fvc::flux
00019             (
00020                 -fvc::flux(-phir, scalar(1) - alpha1, alpharScheme),
00021                 alpha1,
00022                 alpharScheme
00023             );
00024 
00025         MULES::explicitSolve(alpha1, phi, phiAlpha, 1, 0);
00026 
00027         rhoPhi = phiAlpha*(rho1 - rho2) + phi*rho2;
00028     }
by means of this solving and Eq. 31 a new flux is assembled in order
to solve the mommentum equation.
Hope these insights can be helpful, and will be valuable if other
foamers correct me and add more info. Please share your findings.

Regards
Dear S M Damian and others,

Since this post is very old and some of you might have dig into the code. Hence, I am posting my doubt in this thread.

If we look at the MULES solver as pointed by SM Damian, new flux is assembled for the momentum equation.

Doubt: In case of interface, the "rhoPhi" term can be easily related to the equation. Now, How to relate the "rhoPhi" term for the phase 1 or phase 2 cells (where no interface is present)?

Please correct me, if I misinterpret anything.

-
Thanks in Advance
Tushar@cfd is offline   Reply With Quote

Old   June 13, 2015, 01:50
Default
  #80
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 319
Blog Entries: 1
Rep Power: 9
Tushar@cfd is on a distinguished road
Dear All,

Thanks again for your time. I have understood the term "rhoPhi". Hence I am sharing this info, may be it will be beneficial to others.

For phase 1:
rhoPhi = phi*rho1;

For phase 2:
rhoPhi = phi*rho2;

For interface:
rhoPhi = phiAlpha*(rho1 - rho2) + phi*rho2;

-
Tushar@cfd 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
Comparing between the Fractional step method and the SIMPLE method ghlee Main CFD Forum 1 April 10, 2012 16:59
Finite volume method vs finite difference method? superfool Main CFD Forum 4 October 21, 2006 14:37
Pressure Correction method, Finite Volume Method Seeker01 Main CFD Forum 2 January 13, 2003 03:49
hess-smith method and fvm method yangqing FLUENT 0 March 20, 2002 20:25
Projection method and Block-off method Leo Main CFD Forum 0 June 14, 2001 07:22


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