
[Sponsors] 
July 16, 2011, 02:10 

#41 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,907
Rep Power: 27 
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 Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

July 16, 2011, 03:03 

#42  
Senior Member
Gijsbert Wierink
Join Date: Mar 2009
Posts: 383
Rep Power: 10 
Hi Tanay,
It would be helpful to post the compile error. But just guessing from Quote:
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); EDIT: Whoops, missed Alberto's post. Well, hope it is helpful anyway
__________________
Regards, Gijs Last edited by gwierink; July 16, 2011 at 03:10. Reason: Missed post 

June 19, 2012, 08:50 

#43  
New Member
Join Date: Oct 2011
Posts: 19
Rep Power: 6 
Quote:
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())))? 

June 19, 2012, 09:01 

#44 
Senior Member
Balkrishna Patankar
Join Date: Mar 2009
Location: Pune
Posts: 122
Rep Power: 9 
According to the documentation , OpenFOAM doc
Code:
A limitedLinear scheme based on a flux phi with is specified as default by phi = 1. Probably , thats why Gauss limited linear doesnt work. 

September 26, 2013, 00:22 
two liquid mixing in a pipe flow

#45 
New Member
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 8 
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. 

September 26, 2013, 11:44 

#46 
Senior Member
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 136
Rep Power: 9 
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 

September 26, 2013, 19:00 

#47 
New Member
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 8 
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 17:35. 

November 4, 2013, 13:59 

#48 
New Member
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 5 
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! 

November 4, 2013, 14:14 

#49 
Senior Member
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 136
Rep Power: 9 
Hi,
In the folder "system" in the file: turbulenceProperties you can set desired turbulence model. If you choose: laminar, then standard NS 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 

November 4, 2013, 14:46 

#50 
New Member
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 5 
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! 

November 4, 2013, 15:53 

#51 
Senior Member
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 136
Rep Power: 9 
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 

November 5, 2013, 11:55 

#52 
New Member
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 5 
Thanks for the answer, that will be a nice help for me to work on.


January 14, 2014, 18:43 
schemes for alpha1

#53 
New Member
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 8 
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!! 

January 22, 2014, 15:39 
twoLiquidMixingDyMFoam

#54 
New Member
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 8 
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?


January 23, 2014, 22:20 
twoLiquidMixingDyMFoam

#55 
New Member
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 8 
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 

February 19, 2014, 14:38 
LES model

#56 
New Member
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 5 
Hello everyone,
How can i implement the LES model in the solver twoLiquidMixingFoam? Thank you very much! 

March 20, 2014, 13:41 
fvSchemes

#57 
New Member
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 5 
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. 

June 25, 2014, 13:26 
twoLiquidMixingDyMFoam

#58 
New Member
Carlos Alberto Flórez Arias
Join Date: Aug 2010
Posts: 8
Rep Power: 8 
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?


September 23, 2014, 13:50 
Front Velocity

#59 
New Member
Christiano Molossi
Join Date: Sep 2013
Posts: 10
Rep Power: 5 
Hello again people,
I was wondering... How do you guys do to calculate the front velocity from a gravity current simulation? 

December 8, 2014, 07:35 
Regarding Kinematic Viscosity

#60  
Member
Himanshu Sharma
Join Date: Jul 2012
Posts: 84
Rep Power: 6 
Hi ziemowitzima,
your reply was useful but since there are two fluids and in the transport property file we specify kinematic viscosity of both the fluid,hence my question is which fluid kinematic viscosity will be used to solve the UEqn. will be weighted same as density > rho = rho_1*alpha1 +rho_2*(1alpha_1) as there are nothing being mentioned on this part. Thank you Quote:


Thread Tools  
Display Modes  

