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

TwoLiquidMixingFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 16, 2011, 03:10
Default
  #41
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,888
Rep Power: 25
alberto will become famous soon enoughalberto will become famous soon enough
The version of twoLiquidMixingFoam in OpenFOAM 2.0 already uses PIMPLE. A tutorial was added today to 2.0.x (in git repository). Maybe take a look.

Best,
__________________
Alberto

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
GeekoCFD 32bit - The 32bit edition of GeekoCFD.
GeekoCFD text mode - A smaller version of GeekoCFD, text-mode only, with only OpenFOAM. Available in a variety of virtual formats.
alberto is offline   Reply With Quote

Old   July 16, 2011, 04:03
Default
  #42
Senior Member
 
Gijsbert Wierink
Join Date: Mar 2009
Posts: 382
Rep Power: 8
gwierink is on a distinguished road
Hi Tanay,

It would be helpful to post the compile error. But just guessing from

Quote:
Originally Posted by Tanay View Post
surfaceScalarField muEff ("muEff", (alpha1*mu1 + (scalar(1)-alpha1)*mu2));
volScalarField rho ("rho", (alpha1*rho1 + (scalar(1)-alpha1)*rho2));
surfaceScalarField rhoIntoPhi ("rhoIntoPhi", rho*phi);
it seems that at least one of the problems is that you're trying to assign a volScalarField (rho) to something that is a surfaceScalarField (rhoIntoPhi). You can make rho into a surfaceScalarField for example by:

Code:
surfaceScalarField muEff ("muEff", (alpha1*mu1 + (scalar(1)-alpha1)*mu2));
volScalarField rho ("rho", (alpha1*rho1 + (scalar(1)-alpha1)*rho2));
surfaceScalarField rhof = fvc::interpolate(rho);
surfaceScalarField rhoIntoPhi ("rhoIntoPhi", rhof*phi);
Hope that is of any help

EDIT: Whoops, missed Alberto's post. Well, hope it is helpful anyway
__________________
Regards, Gijs

Last edited by gwierink; July 16, 2011 at 04:10. Reason: Missed post
gwierink is offline   Reply With Quote

Old   June 19, 2012, 09:50
Default
  #43
New Member
 
Join Date: Oct 2011
Posts: 19
Rep Power: 5
okroud is on a distinguished road
Quote:
Originally Posted by Tanay View Post

Here's the fvSchemes file if anyone wishes to view it-
ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
grad(U) Gauss linear;
grad(muEff) Gauss linear;
grad(p_rgh) Gauss linear;
}

snGradSchemes
{
default corrected;
}

divSchemes
{
default none;
div(phi,U) Gauss limitedLinearV 1;
div(phi,k) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
div(phi,R) Gauss limitedLinear 1;
div(R) Gauss linear;
div(phi,nuTilda) Gauss limitedLinear 1;
div((nuEff*dev(grad(U).T()))) Gauss linear;
div(phi,alpha1) Gauss limitedLinear 1;
div(phi) Gauss limitedLinear 1;
div(rhoPhi,U) Gauss limitedLinear 1;
div(rho*phi,U) Gauss limitedLinear 1;
}
It's long ago, but do you still know why you have set
div((nuEff*dev(grad(U).T()))) Gauss linear;
And not limitedLinear 1 like all the other divSchemes. I know Gauss limitedLinear 1 exactly doesn't work (which is actually my problem), but isn't there a similar Scheme like limitedLinearV1 oder limitedLinear phi that works for div((nuEff*dev(grad(U).T())))?
okroud is offline   Reply With Quote

Old   June 19, 2012, 10:01
Default
  #44
Senior Member
 
Balkrishna Patankar
Join Date: Mar 2009
Location: Pune
Posts: 119
Rep Power: 7
balkrishna is on a distinguished road
According to the documentation , OpenFOAM doc
Code:
A limitedLinear scheme based on a flux phi with  is specified as default by phi = 1.
In ur definition, "div((nuEff*dev(grad(U).T()))) " there is no phi involved.

Probably , thats why Gauss limited linear doesnt work.
balkrishna is offline   Reply With Quote

Old   September 26, 2013, 01:22
Default two liquid mixing in a pipe flow
  #45
New Member
 
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 6
carlosfloar is on a distinguished road
Hi, I'm a master student who works about the mixing of two liquids in batching transport in a long pipes. I need to use a moving mesh because the high computational cost involved in the solutión of the entire domain. I want to use openFOAM for the solution of that problem specificaly employing the solver twoliquidmixingFoam. ¿It is possible to do?

thanks in advance for your help.
carlosfloar is offline   Reply With Quote

Old   September 26, 2013, 12:44
Default
  #46
Senior Member
 
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 118
Rep Power: 7
ziemowitzima is on a distinguished road
What kind of flow, are u planning to model ? Turbulent one, what is your Re ?
Do you mean ADAPTIVE MESH or MOVING mesh (changing geometry of the domain).
Both are possible. Did you look in this forum ?

Best
ziemowitzima is offline   Reply With Quote

Old   September 26, 2013, 20:00
Default
  #47
New Member
 
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 6
carlosfloar is on a distinguished road
I try to model a turbulent flow with Reynolds about 40000, I need to use moving mesh because the extension of the pipe, too many km and my interest is centered in the mixed fluid. At the moment I use a libraries written in c++ to attempt to solve this problem but I want to use OpenFOAM to solve it. ¿Anybody can help me with this? please, documentation, examples, tutorials...

Last edited by carlosfloar; September 27, 2013 at 18:35.
carlosfloar is offline   Reply With Quote

Old   November 4, 2013, 13:59
Default
  #48
New Member
 
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 3
m.Chrsitiano is on a distinguished road
Hi everyone,

How do I know what equations that twoLiquidMixingFoam uses?
For exemple, if I set no turbulence model (laminar), where is the laplacian term in the UEqn? And where the solver uses the kinematic viscosity (nu)?

Any replay would be very helpful, since I've just started using OpenFOAM.

Best regards!
m.Chrsitiano is offline   Reply With Quote

Old   November 4, 2013, 14:14
Default
  #49
Senior Member
 
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 118
Rep Power: 7
ziemowitzima is on a distinguished road
Hi,
In the folder "system" in the file:
turbulenceProperties
you can set desired turbulence model. If you choose: laminar,
then standard N-S equations will be solved without any turbulence model, namely in the U eqn:
fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
+ turbulence->divDevReff(U)
);
the diffusion term
turbulence->divDevReff(U)
will be "replaced" by :
laplacian(nu,U)
where nu is your viscosity

if you choose e.g.:
simulationType RASModel;
then you have to have one more file in "system":
RASproperties
then in
turbulence->divDevReff(U)
corresponding turbulence model will be solved

see some of the examples in the tutorial/run/.../cases
ziemowitzima is offline   Reply With Quote

Old   November 4, 2013, 14:46
Default
  #50
New Member
 
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 3
m.Chrsitiano is on a distinguished road
Thank you very much! It was very enlightening!

Now make sense to me.
One more thing, about transport equation, where am I able to find it? Where is the file that it is writing?

Thank you very much again for your attention!
m.Chrsitiano is offline   Reply With Quote

Old   November 4, 2013, 15:53
Default
  #51
Senior Member
 
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 118
Rep Power: 7
ziemowitzima is on a distinguished road
every solver is located at:
/opt/openfoam210/applications/solvers

or sth similar depending what version of OF and what OS you have.

see
http://openfoamwiki.net/index.php/Main_Page
and
http://openfoamwiki.net/index.php/IcoFoam

best
ziemowitzima is offline   Reply With Quote

Old   November 5, 2013, 11:55
Default
  #52
New Member
 
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 3
m.Chrsitiano is on a distinguished road
Thanks for the answer, that will be a nice help for me to work on.
m.Chrsitiano is offline   Reply With Quote

Old   January 14, 2014, 18:43
Default schemes for alpha1
  #53
New Member
 
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 6
carlosfloar is on a distinguished road
hi
I want to ask about the alphaDiffusionEqn.H file, in that file what's the meaning of the second line in:

fvm::ddt(alpha1)
- fvc::ddt(alpha1)
- fvm::laplacian
(
volScalarField("Dab", Dab + alphatab*turbulence->nut()),
alpha1
)
I think the meaning of the first line, fvm::ddt(alpha1) is (((rho_P*alpha1_P)^n)-((rho_P*alpha1_P)^o))/(Delta t)
it's this right? and, the second line?
please help me!!
carlosfloar is offline   Reply With Quote

Old   January 22, 2014, 15:39
Default twoLiquidMixingDyMFoam
  #54
New Member
 
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 6
carlosfloar is on a distinguished road
hi, I'm working in the implementation of Dynamic mesh in twoLiquidMixingFoam solver, I'm stuck at the moment, I try to use InterDyMFoam for a guide but don't obtain any results....any suggestions?
carlosfloar is offline   Reply With Quote

Old   January 23, 2014, 22:20
Default twoLiquidMixingDyMFoam
  #55
New Member
 
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 6
carlosfloar is on a distinguished road
I want to simulate a parcel composed of two fluids moving with the mean flow velocity, I try to implement DyM twoLiquidMixingFoam, but when I try to simulate a case the delta of time begins to decrease until it finally stops the solution ... I use a Courant=0.1 ,I initially tested on lockForExchange with boundary conditions fixedValue=0 on the walls ..... any suggestions?

Last edited by carlosfloar; January 23, 2014 at 22:20. Reason: new title
carlosfloar is offline   Reply With Quote

Old   February 19, 2014, 14:38
Default LES model
  #56
New Member
 
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 3
m.Chrsitiano is on a distinguished road
Hello everyone,

How can i implement the LES model in the solver twoLiquidMixingFoam?

Thank you very much!
m.Chrsitiano is offline   Reply With Quote

Old   March 20, 2014, 14:41
Default fvSchemes
  #57
New Member
 
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 3
m.Chrsitiano is on a distinguished road
Hey guys,

Nobody know how to setup the solver to LES?

Another thing is, how do i change the numerical schemes in the fvSchems file? For exemple, the InterpolationSchemes?

Please, i would be very greatful for any replay.

Thank you in advance,

Christiano.
m.Chrsitiano is offline   Reply With Quote

Old   June 25, 2014, 14:26
Default twoLiquidMixingDyMFoam
  #58
New Member
 
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 6
carlosfloar is on a distinguished road
Hi, I implemented the funcionality of dynamic mesh in the solver twoLiquidMixingFoam to solve the mixing of two fluids in a pipeline tracking a control volume of the domain large enough so that the interface between the two liquids is always contained therein. I ran tests with both sovers, twoLiquidMixingFoam and "twoLiquidMixingDyMFoam" but when comparing the amount of liquid that is mixed, it is rather more with the solver twoLiquidMixingDyMFoam. I use the same configuration in both solvers, the same fvSchemes, fvSolution, etc. anyone can help me with the Dynamic Mesh error?
carlosfloar is offline   Reply With Quote

Old   September 23, 2014, 14:50
Default Front Velocity
  #59
New Member
 
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 3
m.Chrsitiano is on a distinguished road
Hello again people,

I was wondering...

How do you guys do to calculate the front velocity from a gravity current simulation?
m.Chrsitiano 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 16:28.