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

Divergence problem related to interMixingFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 5, 2016, 05:38
Default Divergence problem related to interMixingFoam
  #1
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
Hi Foamers,

I am doing a simulation related to 3 incompressible fluids (ethanol, water and air), two of them are misible, from the Openfoam User Guide, I used interMixingFoam, my concentration is on the mixing process of water and ethanol. The physical process is firstly put the ethanol and water in a cup that have a inlet hole in the bottom, and then pump air into the cup from the bottom to help the mixing process (the geometry of the cup is shown in attatched pictures, the second picture shows that the ethanol is initially at the center of the water), the air is pumped into the cup within a specified time interval (which means the pressure at inlet is a time-dependent variable, and the inlet velocity should be calculated via this inlet pressure).

At the first time, I create the cup geometry by blockMeshDict and setup the varying pressure by uniformTotalPressure (I'm not familiar about this boundary condition, before this I have attempt timeVaryingTotalPressure, timeVaryingUniformFixedValue, timeVaryingMappedFixedValue and uniformFixedValue but all failed, the Openfoam returns an error indicate that there is no such a BC for the former two, I used Openfoam 2.3.1, for the latter two, I don't know how to use them since I didn't found a clear tutorial), after I start the simulation, it diverges at the begining and I don't know how to fix this problem.

After that, I tried to change the time-dependent boundary condition since I dont't familiar about it, I modified the inlet pressure to a fixed value and restart the simulation, but the diverge problem still exist. I have tried refine the mesh, decrease (or increase) the Courant number, decrease the deltaT value but all failed.

Then I have changed the geometry of the cup to a simple cylinder, and the simulation converges in the first 0.06s but diverges after that (under the fixed inlet pressure BC).

I have attatched the zip files of my simulation code (include varying pressure and fixed pressure) , could anyone help me to check my case file? I don't where is the problem, and if possible, could you help me to check if my time-dependent pressure BC is correct? Any suggestion will be appreaciated!

Best,

weicent
Attached Images
File Type: jpg forum 1.jpg (20.1 KB, 59 views)
File Type: png forum 2.png (23.1 KB, 49 views)
Attached Files
File Type: zip ethanolWater_fixedValue.zip (11.7 KB, 12 views)
File Type: zip ethanolWater_timeDependent.zip (10.7 KB, 7 views)
weicent is offline   Reply With Quote

Old   September 5, 2016, 22:29
Default
  #2
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
the below is the message come from terminal window, as you can see, the courant number is extremly large, which is not normal.

Courant Number mean: 3250.96 max: 2.49028e+08
Interface Courant Number mean: 0 max: 0
deltaT = 2.2309e-12
--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1055
Increased the timePrecision from 6 to 7 to distinguish between timeNames at time 0.00111111
Time = 0.001111111

--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1055
Increased the timePrecision from 7 to 8 to distinguish between timeNames at time 0.00111111
diagonal: Solving for alpha.air, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for alpha.other, Initial residual = 1.05882e-06, Final residual = 8.58729e-18, No Iterations 1
Air phase volume fraction = 0.60692 Min(alpha1) = 0 Max(alpha1) = 1.1154
Liquid phase volume fraction = 0.605912 Min(alpha2) = 0 Max(alpha2) = 1.1154
weicent is offline   Reply With Quote

Old   September 5, 2016, 22:30
Default
  #3
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
the below is the message come from terminal window, as you can see, the courant number is extremly large, which is not normal.

Courant Number mean: 3250.96 max: 2.49028e+08
Interface Courant Number mean: 0 max: 0
deltaT = 2.2309e-12
--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1055
Increased the timePrecision from 6 to 7 to distinguish between timeNames at time 0.00111111
Time = 0.001111111

--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1055
Increased the timePrecision from 7 to 8 to distinguish between timeNames at time 0.00111111
diagonal: Solving for alpha.air, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for alpha.other, Initial residual = 1.05882e-06, Final residual = 8.58729e-18, No Iterations 1
Air phase volume fraction = 0.60692 Min(alpha1) = 0 Max(alpha1) = 1.1154
Liquid phase volume fraction = 0.605912 Min(alpha2) = 0 Max(alpha2) = 1.1154
weicent is offline   Reply With Quote

Old   September 7, 2016, 05:28
Default
  #4
Member
 
Join Date: May 2016
Posts: 39
Rep Power: 8
dzordz is on a distinguished road
Try this.

For /U BC:
walls
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}
inlet
{
type fixedValue;
value uniform (0 0 1);
}

For /p BC:
walls
{
type fixedFluxPressure;
}
inlet
{
type zeroGradient;
}
outlet
{
type totalPressure;
p0 uniform 0;
U U;
phi phi;
rho rho;
psi none;
gamma 1;
value uniform 0;
}

Other is good. Hope it works with this combination
dzordz is offline   Reply With Quote

Old   September 7, 2016, 21:54
Default
  #5
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
Hi dzordz;

Thank you very much for you reply!

I have tried your code, but the problem still exist, I put it below, its actually still diverges:

Courant Number mean: 0.000236288 max: 7.53633
Interface Courant Number mean: 0 max: 0
deltaT = 1.68055e-12
--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1055
Increased the timePrecision from 7 to 8 to distinguish between timeNames at time 4.56634e-05
Time = 4.5663439e-05

diagonal: Solving for alpha.air, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for alpha.other, Initial residual = 4.12149e-10, Final residual = 4.12149e-10, No Iterations 0
Air phase volume fraction = 0.392911 Min(alpha1) = 0 Max(alpha1) = 1
Liquid phase volume fraction = 0.60608 Min(alpha2) = -5.37996e-73 Max(alpha2) = 1.00064
--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1055
Increased the timePrecision from 8 to 9 to distinguish between timeNames at time 4.56634e-05
diagonal: Solving for alpha.air, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for alpha.other, Initial residual = 4.12149e-10, Final residual = 4.12149e-10, No Iterations 0
Air phase volume fraction = 0.392911 Min(alpha1) = 0 Max(alpha1) = 1
Liquid phase volume fraction = 0.60608 Min(alpha2) = -5.37996e-73 Max(alpha2) = 1.00064
DICPCG: Solving for p_rgh, Initial residual = 0.00332591, Final residual = 0.00015166, No Iterations 9
DICPCG: Solving for p_rgh, Initial residual = 0.00450806, Final residual = 0.000213074, No Iterations 7
DICPCG: Solving for p_rgh, Initial residual = 0.00645848, Final residual = 0.000316934, No Iterations 7
time step continuity errors : sum local = 4.90275e-07, global = 4.19724e-14, cumulative = -4.96021e-05
DICPCG: Solving for p_rgh, Initial residual = 0.00844886, Final residual = 0.000409655, No Iterations 7
DICPCG: Solving for p_rgh, Initial residual = 0.0127241, Final residual = 0.000561724, No Iterations 8
DICPCG: Solving for p_rgh, Initial residual = 0.0162193, Final residual = 0.000758191, No Iterations 9
time step continuity errors : sum local = 1.17456e-06, global = -5.09239e-15, cumulative = -4.96021e-05
DICPCG: Solving for p_rgh, Initial residual = 0.0242594, Final residual = 0.00102831, No Iterations 11
DICPCG: Solving for p_rgh, Initial residual = 0.0363386, Final residual = 0.00165169, No Iterations 11
DICPCG: Solving for p_rgh, Initial residual = 0.0423467, Final residual = 9.83501e-08, No Iterations 166
time step continuity errors : sum local = 1.48924e-10, global = 1.54726e-11, cumulative = -4.96021e-05
time step continuity errors : sum local = 1.48924e-10, global = 1.54726e-11, cumulative = -4.96021e-05
ExecutionTime = 17.87 s ClockTime = 17 s

thanks for your kind help again;

weicent
weicent is offline   Reply With Quote

Old   September 7, 2016, 22:06
Default
  #6
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
I actually found some problems in setFieldsDict, before using the interMixingFoam, I have been working with interFoam, which is a two phase solver, so this is the first time I tried to manage three fluids in setFieldsDict by VOF method, I misunderstand something when modify this dictionary, now I have changed it, but the diverge problem still exist, I put the new setFieldsDict here so that someone can check my case to see where is the diverge problem.
Attached Files
File Type: zip setFieldsDict.zip (848 Bytes, 14 views)
weicent is offline   Reply With Quote

Old   September 8, 2016, 03:17
Default
  #7
Member
 
Join Date: May 2016
Posts: 39
Rep Power: 8
dzordz is on a distinguished road
Quote:
Originally Posted by weicent View Post

Increased the timePrecision from 8 to 9 to distinguish between timeNames

weicent
Since this happens I would try raising time and write precision in ControlDict to 8 for example.
dzordz is offline   Reply With Quote

Old   September 8, 2016, 05:50
Default
  #8
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
Quote:
Originally Posted by dzordz View Post
Since this happens I would try raising time and write precision in ControlDict to 8 for example.
Hi dzordz;

thanks for your help, I have made some modify in alpha.air, alpha.water and alpha.other, and also in geometry. I have delete the arc edges in blockMeshDict now, since the curved geometry make the mesh quality is not good, I'm planning to establish this curved geometry later on, but before that, I need to check if the modified case can running well, actually, it seems converges untill 0.43, and the simulation still running I'm waiting to see if it can finally yield a stable result.

below is the new alpha.xx documents I used:

PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volScalarField;
    
location    "0";
    
object      alpha.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    
inlet
    
{
        
type            inletOutlet;
        
inletValue      uniform 1;
        
value           uniform 1;
    }
    
outlet
    
{
        
type            zeroGradient;
    }
    
walls
    
{
        
type            zeroGradient;
    }
}

// ************************************************************************* // 

PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volScalarField;
    
location    "0";
    
object      alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    
inlet
    
{
        
type            zeroGradient;
    }
    
outlet
    
{
        
type            zeroGradient;
    }
    
walls
    
{
        
type            zeroGradient;
    }
}

// ************************************************************************* // 
PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volScalarField;
    
location    "0";
    
object      alpha.other;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;


boundaryField
{
    
inlet
    
{
        
type            inletOutlet;
        
inletValue      uniform 0;
        
value           uniform 0;
    }
    
outlet
    
{
        
type            zeroGradient;
    }
    
walls
    
{
        
type            zeroGradient;
    }
}

// ************************************************************************* // 
weicent is offline   Reply With Quote

Old   September 8, 2016, 06:14
Default
  #9
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
the cylinder geometry and mesh should be created like the first and second picture, while my geometry is shown in the third picture, which i think probably a mistake which lead to simulation divergence.

I will try to modify the new blockMeshDict later and update the outcomes.
Attached Images
File Type: jpg cylinder2.jpg (22.3 KB, 20 views)
File Type: jpg cylinderMesh.jpg (127.4 KB, 25 views)
File Type: jpg mesh in top surface.jpg (94.6 KB, 19 views)
weicent is offline   Reply With Quote

Old   September 8, 2016, 06:16
Default
  #10
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
Quote:
Originally Posted by weicent View Post
the cylinder geometry and mesh should be created like the first and second picture, while my geometry is shown in the third picture, which i think probably a mistake which lead to simulation divergence.

I will try to modify the new blockMeshDict later and update the outcomes.
sorry, i just forget put my mistake mesh screenshot, here is it:
Attached Images
File Type: jpg mesh in top surface.jpg (94.6 KB, 19 views)
weicent is offline   Reply With Quote

Old   September 8, 2016, 06:39
Default
  #11
Member
 
Join Date: May 2016
Posts: 39
Rep Power: 8
dzordz is on a distinguished road
What output do you get with checkMesh? There might be some skew faces or big non-orthogonal cells that are messing up your simulation.
dzordz is offline   Reply With Quote

Old   September 8, 2016, 22:17
Default
  #12
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
Quote:
Originally Posted by dzordz View Post
What output do you get with checkMesh? There might be some skew faces or big non-orthogonal cells that are messing up your simulation.
Hi dzorzd:

below is my checkMesh output, thanks for help:


PHP Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh 
for time 0

Time 
0

Mesh stats
    points
:           19481
    faces
:            51300
    internal faces
:   44700
    cells
:            16000
    faces per cell
:   6
    boundary patches
3
    point zones
:      0
    face zones
:       0
    cell zones
:       0

Overall number of cells of each type
:
    
hexahedra:     16000
    prisms
:        0
    wedges
:        0
    pyramids
:      0
    tet wedges
:    0
    tetrahedra
:    0
    polyhedra
:     0

Checking topology
...
    
Boundary definition OK.
    
Cell to face addressing OK.
    
Point usage OK.
    
Upper triangular ordering OK.
    
Face vertices OK.
    
Number of regions(OK).

Checking patch topology for multiply connected surfaces...
    
Patch               Faces    Points   Surface topology                  
    inlet               100      121      ok 
(non-closed singly connected)  
    
outlet              100      121      ok (non-closed singly connected)  
    
walls               6400     6440     ok (non-closed singly connected)  

Checking geometry...
    
Overall domain bounding box (-0.00636396 0 -0.00636396) (0.00636396 0.077 0.00636396)
    
Mesh (non-empty, non-wedgedirections (1 1 1)
    
Mesh (non-empty) directions (1 1 1)
    
Boundary openness (-4.83852e-18 3.40064e-16 8.22549e-18OK.
    
Max cell openness 2.50633e-16 OK.
    
Max aspect ratio 4.71405 OK.
    
Minimum face area 1.125e-08Maximum face area 1.62e-06.  Face area magnitudes OK.
    
Min volume 4.27603e-12Max volume 8.91e-10.  Total volume 6.91938e-06.  Cell volumes OK.
    
Mesh non-orthogonality Max43.7185 average7.73539
    Non
-orthogonality check OK.
    
Face pyramids OK.
    
Max skewness 1.21537 OK.
    
Coupled point location match (average 0OK.

Mesh OK.

End 
weicent is offline   Reply With Quote

Old   September 9, 2016, 01:26
Default
  #13
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
Quote:
Originally Posted by dzordz View Post
Try this.

For /U BC:
walls
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}
inlet
{
type fixedValue;
value uniform (0 0 1);
}

For /p BC:
walls
{
type fixedFluxPressure;
}
inlet
{
type zeroGradient;
}
outlet
{
type totalPressure;
p0 uniform 0;
U U;
phi phi;
rho rho;
psi none;
gamma 1;
value uniform 0;
}

Other is good. Hope it works with this combination
Hi dzordz:

After last time I modify the alpha.xxx and blockMesh, the simulation seems converges, while after 0.82s, it still diverge.

Then I made some modify again and combined your P/U code and finally it worked, the simulation now looks very good and paraview shows a satisfied result, thank you my friend!

For summary, below is what I modified from last time:

1. I have changed the inletValue and value for alpha.air and alhpa.other

for alpha.air, the updated code is:

inlet
{
type inletOutlet;
inletValue uniform 0.5;
value uniform 0.5;
}


for alpha.other, the updated code is:

inlet
{
type inletOutlet;
inletValue uniform 0.5;
value uniform 0;
}


2. I have replace the P/U boundary condition with your code.

3. I have disable the arc edges

4. I set the inlet velocity as 'uniform (0 0.5 0)', which is 0.5m/s



But I actually don't understand what does 'value' specified? In my opinion, in alpha.xxx file, the 'inletValue' command specify the volume phase fraction at inlet, 1 represents inlet is full of that phase, 0 means there is no that phase at inlet (I'm not sure if I am right), so what on earth the 'value'command specify?


Best;

weicent
weicent is offline   Reply With Quote

Old   September 9, 2016, 03:19
Default
  #14
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
Quote:
Originally Posted by weicent View Post
Hi dzordz:

After last time I modify the alpha.xxx and blockMesh, the simulation seems converges, while after 0.82s, it still diverge.

Then I made some modify again and combined your P/U code and finally it worked, the simulation now looks very good and paraview shows a satisfied result, thank you my friend!

For summary, below is what I modified from last time:

1. I have changed the inletValue and value for alpha.air and alhpa.other

for alpha.air, the updated code is:

inlet
{
type inletOutlet;
inletValue uniform 0.5;
value uniform 0.5;
}


for alpha.other, the updated code is:

inlet
{
type inletOutlet;
inletValue uniform 0.5;
value uniform 0;
}


2. I have replace the P/U boundary condition with your code.

3. I have disable the arc edges

4. I set the inlet velocity as 'uniform (0 0.5 0)', which is 0.5m/s



But I actually don't understand what does 'value' specified? In my opinion, in alpha.xxx file, the 'inletValue' command specify the volume phase fraction at inlet, 1 represents inlet is full of that phase, 0 means there is no that phase at inlet (I'm not sure if I am right), so what on earth the 'value'command specify?


Best;

weicent
The updated code works good, but still diverges after 1.75s, I open the paraview and it shows a good result before 1.75s.

After the simulation diverges, I open the controlDict and change the start time from 0s to 1.75s (also change the 'startFrom' from 'latestTime' to 'startTime') and then restart the simulation, it seems the simulation running again.

I will update if there is any other problem occurs.

Below is some of the diverge messege in terminal window at 1.75s, does anyone know what lead to the below error?

PHP Code:
--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 184 to 185 to distinguish between timeNames at time 1.75539
diagonal
:  Solving for alpha.airInitial residual 0, Final residual 0No Iterations 0
smoothSolver
:  Solving for alpha.otherInitial residual 2.3696e-08, Final residual 2.3696e-08No Iterations 0
Air phase volume fraction 
0.324605  Min(alpha1) = -4.25589e-08  Max(alpha1) = 1.00001
Liquid phase volume fraction 
0.669053  Min(alpha2) = 8.56707e-08  Max(alpha2) = 0.994056
--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 185 to 186 to distinguish between timeNames at time 1.75539
diagonal
:  Solving for alpha.airInitial residual 0, Final residual 0No Iterations 0
smoothSolver
:  Solving for alpha.otherInitial residual 2.48385e-08, Final residual 2.48385e-08No Iterations 0
Air phase volume fraction 
0.324605  Min(alpha1) = -4.25589e-08  Max(alpha1) = 1.00001
Liquid phase volume fraction 
0.669053  Min(alpha2) = 8.56707e-08  Max(alpha2) = 0.994056
--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 186 to 187 to distinguish between timeNames at time 1.75539
DICPCG
:  Solving for p_rghInitial residual 1, Final residual 0.0498738No Iterations 7
DICPCG
:  Solving for p_rghInitial residual 0.406991, Final residual 0.0133185No Iterations 7
DICPCG
:  Solving for p_rghInitial residual 0.215131, Final residual 0.00506178No Iterations 9
time step continuity errors 
sum local 1.17047e-07, global = 2.02128e-16cumulative 0.000179792
DICPCG
:  Solving for p_rghInitial residual 0.297671, Final residual 0.0142627No Iterations 6
DICPCG
:  Solving for p_rghInitial residual 0.189953, Final residual 0.00501451No Iterations 6
DICPCG
:  Solving for p_rghInitial residual 0.104014, Final residual 0.00271376No Iterations 7
time step continuity errors 
sum local 6.37272e-08, global = 1.14816e-16cumulative 0.000179792
DICPCG
:  Solving for p_rghInitial residual 0.0767467, Final residual 0.00226053No Iterations 7
DICPCG
:  Solving for p_rghInitial residual 0.0474809, Final residual 0.00236649No Iterations 5
DICPCG
:  Solving for p_rghInitial residual 0.027869, Final residual 5.98415e-08No Iterations 162
time step continuity errors 
sum local 1.58523e-12, global = 4.33374e-14cumulative 0.000179792
time step continuity errors 
sum local 1.58523e-12, global = 4.33374e-14cumulative 0.000179792
ExecutionTime 
3137.39 s  ClockTime 3152 s

Courant Number mean
1.94382e-06 max0.0831665
Interface Courant Number mean4.44921e-07 max0.0235118
deltaT 
1.51092e-98
--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 187 to 188 to distinguish between timeNames at time 1.75539
Time 
1.75538779576766668100162860355339944362640380859375

--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 188 to 189 to distinguish between timeNames at time 1.75539
diagonal
:  Solving for alpha.airInitial residual 0, Final residual 0No Iterations 0
smoothSolver
:  Solving for alpha.otherInitial residual 1.74632e-08, Final residual 1.74632e-08No Iterations 0
Air phase volume fraction 
0.324605  Min(alpha1) = -4.25588e-08  Max(alpha1) = 1.00001
Liquid phase volume fraction 
0.669053  Min(alpha2) = 8.56707e-08  Max(alpha2) = 0.994056
--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 189 to 190 to distinguish between timeNames at time 1.75539
diagonal
:  Solving for alpha.airInitial residual 0, Final residual 0No Iterations 0
smoothSolver
:  Solving for alpha.otherInitial residual 1.77111e-08, Final residual 1.77111e-08No Iterations 0
Air phase volume fraction 
0.324605  Min(alpha1) = -4.25588e-08  Max(alpha1) = 1.00001
Liquid phase volume fraction 
0.669053  Min(alpha2) = 8.56707e-08  Max(alpha2) = 0.994056
--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 190 to 191 to distinguish between timeNames at time 1.75539
DICPCG
:  Solving for p_rghInitial residual 0.0402945, Final residual 0.0017694No Iterations 13
DICPCG
:  Solving for p_rghInitial residual 0.0373714, Final residual 0.0014712No Iterations 13
DICPCG
:  Solving for p_rghInitial residual 0.0242536, Final residual 0.000964903No Iterations 13
time step continuity errors 
sum local 2.23938e-07, global = -1.23237e-15cumulative 0.000179792
DICPCG
:  Solving for p_rghInitial residual 0.0156767, Final residual 0.000338695No Iterations 14
DICPCG
:  Solving for p_rghInitial residual 0.0309064, Final residual 0.000904599No Iterations 14
DICPCG
:  Solving for p_rghInitial residual 0.0423388, Final residual 0.00153878No Iterations 13
time step continuity errors 
sum local 2.39132e-07, global = -1.06414e-15cumulative 0.000179792
DICPCG
:  Solving for p_rghInitial residual 0.0430841, Final residual 0.00124343No Iterations 14
DICPCG
:  Solving for p_rghInitial residual 0.040241, Final residual 0.00181448No Iterations 16
DICPCG
:  Solving for p_rghInitial residual 0.0282049, Final residual 8.08715e-08No Iterations 188
time step continuity errors 
sum local 2.21815e-11, global = -1.40556e-12cumulative 0.000179792
time step continuity errors 
sum local 2.21815e-11, global = -1.40556e-12cumulative 0.000179792
ExecutionTime 
3137.71 s  ClockTime 3152 s

Courant Number mean
0.000108574 max6.09365
Interface Courant Number mean1.63388e-05 max0.91143
deltaT 
1.23975e-99
--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 191 to 192 to distinguish between timeNames at time 1.75539
Time 
1.75538779576766668100162860355339944362640380859375

--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 192 to 193 to distinguish between timeNames at time 1.75539
diagonal
:  Solving for alpha.airInitial residual 0, Final residual 0No Iterations 0
smoothSolver
:  Solving for alpha.otherInitial residual 3.89654e-08, Final residual 3.89654e-08No Iterations 0
Air phase volume fraction 
0.324605  Min(alpha1) = -4.25582e-08  Max(alpha1) = 1.00001
Liquid phase volume fraction 
0.669053  Min(alpha2) = 8.56707e-08  Max(alpha2) = 0.994056
--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 193 to 194 to distinguish between timeNames at time 1.75539
diagonal
:  Solving for alpha.airInitial residual 0, Final residual 0No Iterations 0
smoothSolver
:  Solving for alpha.otherInitial residual 4.01786e-08, Final residual 4.01786e-08No Iterations 0
Air phase volume fraction 
0.324605  Min(alpha1) = -4.25576e-08  Max(alpha1) = 1.00001
Liquid phase volume fraction 
0.669053  Min(alpha2) = 8.56707e-08  Max(alpha2) = 0.994056
--> FOAM Warning 
    
From function Time::operator++()
    
in file db/Time/Time.C at line 1055
    Increased the timePrecision from 194 to 195 to distinguish between timeNames at time 1.75539
DICPCG
:  Solving for p_rghInitial residual 0.811475, Final residual 0.0353421No Iterations 23
DICPCG
:  Solving for p_rghInitial residual 0.501662, Final residual 0.0248784No Iterations 54
DICPCG
:  Solving for p_rghInitial residual 0.213015, Final residual 0.00977048No Iterations 56
time step continuity errors 
sum local 0.000131555, global = -2.43499e-07cumulative 0.000179549
DICPCG
:  Solving for p_rghInitial residual 0.233642, Final residual 0.0090329No Iterations 120
DICPCG
:  Solving for p_rghInitial residual 0.0983941, Final residual 0.00459639No Iterations 56
DICPCG
:  Solving for p_rghInitial residual 0.238709, Final residual 0.01156No Iterations 56
time step continuity errors 
sum local 7.26265, global = -0.0134776cumulative = -0.013298
DICPCG
:  Solving for p_rghInitial residual 0.235003, Final residual 0.0102599No Iterations 120
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#6  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#7  
 
at ??:?
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  
 
at ??:?
Floating point exception (core dumped

Best;

Weicent
weicent is offline   Reply With Quote

Old   September 9, 2016, 04:15
Default
  #15
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
The simulation finished, I set the endTime equal to 2.5s, the diverge message did not appear again, and the result seems fine, while I still don't know why the diverge problem appears at 1.75s.
weicent is offline   Reply With Quote

Old   September 9, 2016, 05:06
Default
  #16
Member
 
Join Date: May 2016
Posts: 39
Rep Power: 8
dzordz is on a distinguished road
I'm glad that you have been able to solve it. For why divergence happens I also have no clue
dzordz is offline   Reply With Quote

Old   September 9, 2016, 05:18
Default
  #17
New Member
 
weicent
Join Date: Aug 2016
Posts: 13
Rep Power: 8
weicent is on a distinguished road
Quote:
Originally Posted by dzordz View Post
I'm glad that you have been able to solve it. For why divergence happens I also have no clue
thanks dzordz, thanks you very much for your help. :
weicent is offline   Reply With Quote

Old   March 8, 2018, 03:11
Default Specification of 'value' keyword
  #18
New Member
 
Saicharan
Join Date: Jan 2018
Location: Bangalore, India
Posts: 29
Rep Power: 6
wavefunction is on a distinguished road
Quote:
Originally Posted by weicent View Post
But I actually don't understand what does 'value' specified? In my opinion, in alpha.xxx file, the 'inletValue' command specify the volume phase fraction at inlet, 1 represents inlet is full of that phase, 0 means there is no that phase at inlet (I'm not sure if I am right), so what on earth the 'value'command specify?
I remember reading somewhere that the keyword 'value' refers to the value used for initialisation. The keyword 'inletValue' refers to the fixedValue that is assigned if inlet flow occurs.
wavefunction is offline   Reply With Quote

Reply

Tags
diverge, intermixingfoam, multiphase flow, openfoam 2.3.1, time-dependent

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem of calculating divergence due to wrong boundary vaule of finite area mesh lzhou OpenFOAM Programming & Development 2 June 18, 2015 06:10
divergence problem santu Main CFD Forum 0 October 31, 2014 01:58
Divergence detected in pressure correction term problem calculation lm_laxman FLUENT 0 September 23, 2014 01:28
Multiphase phase (gas-solid) flow using Eulerian-Granular medel ( divergence problem) jessie FLUENT 3 May 29, 2014 12:05
[ANSYS Meshing] Inflation - Airfoil Meshing - Divergence problem dalecooper ANSYS Meshing & Geometry 0 July 19, 2013 08:04


All times are GMT -4. The time now is 05:39.