|
[Sponsors] |
October 24, 2008, 11:35 |
Hi,
I want to use interPhas
|
#1 |
New Member
Matthias Hofmann
Join Date: Mar 2009
Posts: 2
Rep Power: 0 |
Hi,
I want to use interPhaseChangeFoam on my case, but I can't find any information about the coefficients needed for the Kunze- Merkle- or SchnerrSauer phase change models, nor can I find any information about the models themselves. Can anyone please tell me the dimensions of the coefficients? Some typical values would be great too. Or does someone know a link which could help me? Best regards Matthias |
|
June 30, 2009, 13:05 |
|
#2 |
Senior Member
isabel
Join Date: Apr 2009
Location: Spain
Posts: 171
Rep Power: 17 |
Here you have a tutorial in wich perhaps you find interesting information:
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2008/NaixianLu/REPORT_interPhaseC hangeFoam.pdf |
|
July 15, 2009, 05:34 |
dimension of coefficients
|
#3 |
New Member
Martina Friedrich
Join Date: Jul 2009
Posts: 5
Rep Power: 17 |
The dimensions of the coefficients are:
KunzCoeffs { Cc Cc [0 0 0 0 0 0 0] 1000; Cv Cv [0 0 0 0 0 0 0] 10000; UInf UInf [0 1 -1 0 0 0 0] 6; tInf tInf [0 0 1 0 0 0 0] 1; } MerkleCoeffs { Cc Cc [0 0 0 0 0 0 0] 1000; Cv Cv [0 0 0 0 0 0 0] 10000; UInf UInf [0 1 -1 0 0 0 0] 6; tInf tInf [0 0 1 0 0 0 0] 1; } SchnerrSauerCoeffs { Cc Cc [0 0 0 0 0 0 0] 1000; Cv Cv [0 0 0 0 0 0 0] 10000; n n [0 -3 0 0 0 0 0] 10000; dNuc dNuc [0 1 0 0 0 0 0] 1e-06; } The values of the coefficients depends on the simulation case. |
|
July 15, 2009, 05:37 |
dimension of coefficients
|
#4 |
New Member
Martina Friedrich
Join Date: Jul 2009
Posts: 5
Rep Power: 17 |
cavitation
{ pSat pSat [1 -1 -2 0 0 0 0] -18000; restart no; rampN 200; startN 10000; } |
|
March 1, 2012, 20:42 |
Singhal model
|
#5 |
Member
vahid
Join Date: Feb 2012
Location: Mashhad-Iran
Posts: 80
Rep Power: 13 |
Hi !!!
I transform Saur model to Zwart model successfully but transforming Saur model to Singhal model is faild!!! I see this error : phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C:48: error: no matching function for call to ‘Foam::dimensioned<double>::dimensioned()’ Can you help me??? |
|
March 2, 2012, 16:58 |
AW: Singhal model
|
#6 |
New Member
Martina Friedrich
Join Date: Jul 2009
Posts: 5
Rep Power: 17 |
Hi,
to help you, having your implementation would be great. Within the data you attached, I did not find any "Singhal model" implementation. Or have you used the phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.c, where also your error is located? By the way: a documentation of the "Singhal model" could also help. Best regards, MartinaF |
|
March 7, 2013, 14:38 |
|
#7 |
Senior Member
|
Hi,
I would like also to implement th Singhal model in OpenFoam. Which is the best solver to use with this model? Moreover, I have to account for liquid, vapour and some non condensable gas (dissolved and undissolved). I am quite new in implementing a new model in OpenFoam. What are the steps that I have to follow? Regards Marco |
|
March 8, 2013, 01:11 |
Re.
|
#8 |
Member
vahid
Join Date: Feb 2012
Location: Mashhad-Iran
Posts: 80
Rep Power: 13 |
Hi.
I do this completly. Please Ask your question more clearly? |
|
March 8, 2013, 03:23 |
|
#9 |
Senior Member
|
Hi Vahid,
Have you already done it? I try to explain myself better...I would like to implement a cavitation model in OpenFOAM! The cavitation model is the FULL CAVITATION MODEL developed by Singhal at al. The reference are : INDUSTRIAL TWO PHASE FLOW (Full cavitation model) Von karman Institute `Mathematical basis and validaion of the full cavitation model` Singhal, A. I have never implemented a new model in OpenFOAM! The interPhaseChangeFoam solver is the right solver to do it? My model will account for two types of cavitation: the vapuor cavitation (due to a decrease of the pressure under the equilibrium vapor pressure) and the gas-cavitation (due to a release of the dissolved non condensable gas from the liquid) ! So my mixture would be liquid+vapour+dissoved nonCondensableGas+ undissoved nonCondensableGas. Does the interPhaseChangeFoam account just for two phase? It is difficult extend it for more than two components? Moreover, I would like to account for mass transfer between liquid and vapour as well as dissolved and undissolved Gas. Anyway, as you understood I `m quite new with such staff! So please can you give some hints or outlines to start my project? Thanks in advance Regards Marco |
|
March 10, 2013, 04:46 |
Re.
|
#10 |
Member
vahid
Join Date: Feb 2012
Location: Mashhad-Iran
Posts: 80
Rep Power: 13 |
Hi dear Marco.
yes,I made this solver in five or six month ago. and its working very very good. Now I'm working on a paper,Next week i can get more help for you!!! Please send me your email? |
|
March 11, 2013, 02:42 |
|
#11 | |
Senior Member
|
Goodmorning Vahid,
it sounds great!!! I am very happy that you can share some knowledge with me!! I am very in trouble with this model! My mail is: Quote:
Thank you very much! Regards Marco |
||
April 3, 2013, 13:10 |
|
#12 |
New Member
Join Date: Apr 2011
Posts: 3
Rep Power: 15 |
Hi all,
I already invested a lot of time, figuring out the source code of interPhaseChangeFoam and still there are some essential details I can't find any information to. In the UEqn.H, there is the follwing equation i try to understand: Code:
fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U) + turbulence->divDevRhoReff(rho, U) );
Diego |
|
April 4, 2013, 04:18 |
|
#13 |
Senior Member
|
Hi Diego,
"rho" is the mixture density! The Sp term is used to create a larger diagonal term (to aid the solver) so it only can be used if the linearization of the source term has a negative dependency on the variable being solved for! The turbulence term accounts for the turbulence correction of the viscosity! I hope that it helps you Regards Marco |
|
April 5, 2013, 08:12 |
|
#14 |
New Member
Join Date: Apr 2011
Posts: 3
Rep Power: 15 |
Hi Marco,
thank you for the answer, it helped me indeed. Just one more thing: Is it correct to say, that the Sp-term has no physical significance? In comparison with the momentum equation, I can not find any analogy to this part. Would it therefore be valid to comment it out, as long as convergence is assured? Diego |
|
April 8, 2013, 09:29 |
|
#15 |
Member
ABE
Join Date: Jul 2012
Posts: 46
Rep Power: 14 |
Hi Diego,
You are right. I think Sp had been added to improve convergence and coupling between equations. PS: in OF22, it has been removed! I have a question from you, which OF version do you use (2.1 or oledr)? Thank you in advance ABE |
|
April 8, 2013, 09:48 |
|
#16 |
Senior Member
|
Hi Diego,
Abe and yo are right! Sorry For the late response! Regards Marco |
|
April 9, 2013, 04:33 |
|
#17 |
New Member
Join Date: Apr 2011
Posts: 3
Rep Power: 15 |
I am using 2.1. Thats interesting, I will try simulating without Sp in 2.1. Thanks for the good advice!
Diego |
|
June 6, 2013, 11:49 |
|
#18 |
Senior Member
|
Hallo Faomers,
I am using interPhaseChangeFoam (OpenFoam 2.2.x) to simulate cavitatation in an orifice (I have implemented the singhal mass transfer model)! Unfortunatly, my simulation blows uo due to too high turbulent variable values (kEpsilon model) Have abyone experice to set boundary condition for tutrbulent on interPhaseChangeFoam Regards Marco |
|
June 6, 2013, 11:53 |
|
#19 |
Senior Member
|
Hi Foamers,
Sorry I was wrong and I posted the same post two times! Here are my turbulent boundary conditions dimensions [0 2 -3 0 0 0 0]; internalField uniform $turbulentEpsilon; boundaryField { OUTLET { type zeroGradient; } WALLS { type epsilonWallFunction; value $internalField; } SYM1 { type symmetryPlane; } SYM2 { type symmetryPlane; } INLET { type fixedValue; value $internalField; } #include "include/initialConditions" dimensions [0 2 -2 0 0 0 0]; internalField uniform $turbulentKE; boundaryField { OUTLET { type zeroGradient; } WALLS { type kqRWallFunction; value $internalField; } SYM1 { type symmetryPlane; } SYM2 { type symmetryPlane; } INLET { type fixedValue; value $internalField; Thanks Regards Marco |
|
April 28, 2015, 05:16 |
|
#20 | |
New Member
zhouhoucun
Join Date: Dec 2014
Posts: 12
Rep Power: 11 |
Quote:
I did the same job with you, I have add the Zwart and Singhal model into this solver without any error.But the simulation results of hydrofoil NACA0015 are quiet different from that the Kunz and Schnerr model. Could you give me some advices? Or could you send me your models? My email address is 470861844@qq.com Thank you in advance! |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Looking for interPhaseChangeFoam tutorial of 15 | zjucfd | OpenFOAM Running, Solving & CFD | 12 | April 11, 2013 04:01 |
CFD for aerodynamics coefficients | amalahama | Main CFD Forum | 0 | February 4, 2008 07:45 |
LRR coefficients | amelia | OpenFOAM Pre-Processing | 1 | October 18, 2005 12:06 |
a coefficients | cmv | Siemens | 1 | March 13, 2005 14:49 |
Coefficients | Henrique Argentieri | Phoenics | 0 | August 5, 2003 10:53 |