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

twoLiquidMixingFoam: alpha becomes unbounded

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 23, 2018, 07:45
Default twoLiquidMixingFoam: alpha becomes unbounded
  #1
Member
 
Join Date: Apr 2017
Posts: 68
Rep Power: 9
kaaja is on a distinguished road
Hi.
I get an alpha at 1.6.
I think alpha is supposed to be bounded between 0 and 1.
I have the following settings in fvSolutions:

Code:
"alpha.water.*"
    {
        nAlphaCorr      1;
        nAlphaSubCycles 2;
        cAlpha          1;
        icAlpha         0; 

        MULESCorr       yes;
        nLimiterIter    10;
        alphaApplyPrevCorr  yes;

        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-8;
        relTol          0;
        minIter         1;
    }
Does anybody have an idea about how to solve this?
kaaja is offline   Reply With Quote

Old   March 23, 2018, 09:09
Default
  #2
Senior Member
 
piu58's Avatar
 
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15
piu58 is on a distinguished road
This is far to less information. We need your case for telling something valuable.
__________________
Uwe Pilz
--
Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950)
piu58 is offline   Reply With Quote

Old   March 23, 2018, 09:22
Default
  #3
Member
 
Join Date: Apr 2017
Posts: 68
Rep Power: 9
kaaja is on a distinguished road
Quote:
Originally Posted by piu58 View Post
This is far to less information. We need your case for telling something valuable.
Hi.
Maybe a picture can explain the case somewhat. The following shows my alpha, scaled from 1 to its max value (so the problem area shows easily)

alphaProblem.jpg

It is a flow of two fluids in a tank, and there is an obstacle at the top.
The flow is going from right to left.

Edit: I might add that the phase equal to 1 is at the bottom, so the part that is larger than 1 comes from nowhere.

Maybe not relevant, but in case:
The top BC for U infront of the obstacle is "inletOutlet" with internalField. In practice only the zeroGradient part is activiated in this area, since the flow goes up. The rest of the top as slip BC. This is a result of lots of trial and error. Slip over the whole does not work (Get blowup at a single cell at the corner of the obstacle infront cloeses to the screenm which is a symmetryPlane side). Inlet outlet over the full top gave very strange results.

Last edited by kaaja; March 24, 2018 at 05:32.
kaaja is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Regarding unbounded alpha values -interFoam phani45 OpenFOAM Running, Solving & CFD 8 September 25, 2021 06:05
alphaEqn.H in twoPhaseEulerFoam cheng1988sjtu OpenFOAM Bugs 15 May 1, 2016 16:12
setfields for alpha belkadi OpenFOAM Programming & Development 7 November 25, 2014 11:53
Bounded scheme gives unbounded solution su_junwei OpenFOAM Running, Solving & CFD 3 November 18, 2011 02:35
Determining alpha and beta for porous baffle Liaqat Khan Siemens 1 October 27, 2000 04:44


All times are GMT -4. The time now is 04:29.