interPhaseChangeFoam
hello every1!
any1 familiar with this solver? interPhaseChangeFoam i was wondering whether cavitation takes place on both phases ? (perhaps this question is not applicable at all, sry but im not that familiar with the cavitation phenonemon :( ) *C file states: Code:
******************************* say we have water and oil.. which fits the description of 2 immiscible fluids.. in this case, what would VOF be tracking? the interface water-vapor or water-oil? when postprocessing the respective tutorial, /home/userx/OpenFOAM/userx-2.1.0/run/multiphase/interPhaseChangeFoam , i only get to c the alpha parameter regarding phases http://img252.imageshack.us/img252/5328/caviex.jpg transport dict states: Code:
/*--------------------------------*- C++ -*----------------------------------*\ so can any1 give me some insights? this case is rather complex and slow to carry out many tests.. from what i could glean it seems that what the solver takes for immiscible 2 fluids is the water and its counterpart which appears with cavitation and henceforward the solver uses VOF to track the interface between water and its vapor.... is that right? at least this is what the interpretation of results tells me!.. if so i guess the description of the solver should change a little .. idk.. maybe it's just me :o thanks a lot! |
Calim,
the solver only captures two phases. Therefore you would have to severly extend the solver to capture three or even four phases (oil, water, water- vapor and "oil-vapor"). So at the moment one phase plus its vapor-phase is captured. Best Kathrin |
Quote:
|
kinetic energy turbulence:
please help me…
I want to add the kinetic energy Turbulence(k) in model <<sauer>> in solver <<interPhaseChangeFoam>>. when I Type wmake in terminal. I seen this message : phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C:79: error: ‘k_’ was not declared in this scope make: *** [Make/linux64GccDPOpt/SchnerrSauer.o] Error 1 Although I added the following line in the<<options>>: -I$(LIB_SRC)/turbulenceModels/incompressible/RAS/RASModel \ Then with enter the <<k_>>in model <<sauer>>as following : // * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * // Foam::tmp<Foam::volScalarField> Foam: PhaseChangeTwoPhaseMixtures::SchnerrSauer::rRb ( const volScalarField& limitedAlpha1 ) const { return pow ( ((*4*constant::mathematical: pi*n_)/3)*k_ *limitedAlpha1/(1.0 + alphaNuc() - limitedAlpha1), 1.0/3.0 ); } |
I want to use OF to simulate cavitation. I can get relatively good results with older versions (i.e. 1.5). But by using new versions (2.1 or 2.2), results are not good (I can say awful).
has anybody got good cavitating results with OF2.1 or 2.2? Any idea or comment is highly appreciated, and thank you in advance. ABE |
Same problem here. I created a bug report:
http://www.openfoam.org/mantisbt/view.php?id=817 Best regards David |
Quote:
Thanks for reply, and creating the bug report. ABE |
Thank you david,
I have the same problem and I posted a new thread a week ago. http://www.cfd-online.com/Forums/ope...m-2-2-0-a.html My results (a simple cavitating 2D wing profile) are awful with OF 2.2.x. Instead, with OF 2.1.1 they correlate well with available experiments. Stefano |
I saw your post. Your results with OF2.1.1, generally speaking, seems to be good. However, I have used OF2.1.x and results were not so accurate.
Moreover, I have encountered another issue. Right after restarting the simulation in OF22x, if the amount of vapor be considerable, the pressure field will be disrupted. ABE |
Hi , who know what is the meanings of "cAlpha"? And how to give its numbers ?
|
Hi Sandy,
You can find the description of cAlpha here: http://www.geocities.jp/penguinitis2...interFoam.html cAlpha is a constant representing the degree of compression used to calculate the relative flux, as long as I know. In my cases, I set it equal to zero because I did not consider relative velocity and compressibility of phases interface. ABE |
Quote:
|
Quote:
|
Hi Sandy,
I meant that in the case that you stop the simulation, and then want to continue it, the pressure filed will be completely corrupted in the first iterations of new run. In the small amount of vapor, you may not get in trouble by this issue but when the amount of vapor is considerable, getting convergence would be very difficult. ABE |
Quote:
|
cAlpha1
Quote:
|
From the type of error, it could be anything. I do not have any specific idea now, but I am curious about the differences between values which are stored for next iterations and values which are written out.
When you restart the simulation, I think all of the fluxes, and face values are recalculated based on the cell center values which you had written out, which for sure is different from the continuous simulation. |
The reason for the bad results is a bug. An interim solution is presented in the bug report. The results seem to agree with version 2.0 but are still a bit different from 1.6.
If you have problems directly after the restart, it could help to remove correctPhi from the solver: pimpleControl pimple(mesh); // #include "../interFoam/correctPhi.H" #include "CourantNo.H" #include "setInitialDeltaT.H" Best regards David |
True, I have checked, and I should say that by adding the interim correction (http://www.openfoam.org/mantisbt/view.php?id=817), the problems will be solved (in my case, flow behind a 2d flat plate).
Now the shape of the cavity and pressure fields are more accurate, and besides, the restart issue that I mentioned before has been solved. I just have added a small modification. In my case, in some time steps the MULES solver was not able to bound the alpha value between zero and one (although this values were used as the inputs of the solver) which could lead to divergence. So, I added the following limiter right after the MULES solver as interim solution: alpha1 = min(max(alpha1, scalar(0.0)),scalar(1.0)); About emitting the correctPhi, thanks David I will try that. |
All times are GMT -4. The time now is 21:04. |