[chrisb2244]

Hey,

So I'm currently running a parallel (although this also happens in serial) case using interDyMFoam to dynamically refine my mesh. Ignoring the likely issues with restriction on where I can decompose my mesh to avoid horrible load balancing, the case happily runs for a few iterations, including mesh refinement.

After a little while however (around 170s clocktime), the deltaT value starts dropping alarmingly (before this point it is around 10^-6, which is expected, but drops approximately an order of magnitude each iteration beyond this point), and then timePrecision is increased to compensate for presumably the data being otherwise labelled at the same time.

I am using the adjustableRunTime setup with Courant values limiting the time steps, but the Courant values being output are not equal (or even particularly close to) the limits I have set (0.75 for maxCo, 3 for maxAlphaCo, with maxDeltaT at 1).

If anyone can suggest why my case suddenly goes a little crazy, I would appreciate the help. I have attached the output of the data run and have listed below a smaller section of the output. The first file has the relevant sections (I think) in, but I also attached the second half in case the exit messages are important.

statLog.txt.gz
statLog2.txt.gz

Interface Courant Number mean: 4.10359e-10 max: 0.0069883269
Courant Number mean: 1.4039533e-07 max: 1.894693
deltaT = 1.7523542e-15
Time = 0.002081566360261

Selected 0 cells for refinement out of 97364.
Selected 25 split points out of a possible 10256.
Unrefined from 97364 to 97189 cells.
Execution time for mesh.update() = 0.28 s
time step continuity errors : sum local = 7.3659824e-09, global = -2.3019959e-14, cumulative = -0.0021294309
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 7.7189268e-05, No Iterations 25
time step continuity errors : sum local = 5.6857479e-13, global = -3.0683831e-14, cumulative = -0.0021294309
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 1029
    Increased the timePrecision from 13 to 14 to distinguish between timeNames at time 0.0020815664
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -1.1350829e-22  Max(alpha1) = 1.0000562
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -2.2336362e-22  Max(alpha1) = 1.0000562
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -3.3413735e-23  Max(alpha1) = 1.0000562
GAMG:  Solving for p_rgh, Initial residual = 0.23603412, Final residual = 0.0092550713, No Iterations 4
time step continuity errors : sum local = 1.2120614e-09, global = -1.2161317e-11, cumulative = -0.0021294309
GAMG:  Solving for p_rgh, Initial residual = 0.27790866, Final residual = 0.010788375, No Iterations 5
time step continuity errors : sum local = 1.4657257e-09, global = 7.3197851e-13, cumulative = -0.0021294309
GAMGPCG:  Solving for p_rgh, Initial residual = 0.2120555, Final residual = 5.3478591e-06, No Iterations 21
time step continuity errors : sum local = 1.2998473e-12, global = 3.3009001e-14, cumulative = -0.0021294309
ExecutionTime = 191.1 s  ClockTime = 192 s

Interface Courant Number mean: 1.2856088e-10 max: 0.0018052402
Courant Number mean: 5.0710819e-08 max: 0.29595013
deltaT = 2.102825e-15
Time = 0.0020815663602635

Selected 25 cells for refinement out of 97189.
Refined from 97189 to 97364 cells.
Selected 0 split points out of a possible 10256.
Execution time for mesh.update() = 0.28 s
time step continuity errors : sum local = 7.6474471e-09, global = 3.9610802e-14, cumulative = -0.0021294309
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 9.1571304e-05, No Iterations 18
time step continuity errors : sum local = 7.0028671e-13, global = -9.141037e-14, cumulative = -0.0021294309
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -9.037548e-23  Max(alpha1) = 1.0000562
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -3.5017732e-22  Max(alpha1) = 1.0000562
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -2.4070177e-22  Max(alpha1) = 1.0000562
GAMG:  Solving for p_rgh, Initial residual = 0.17774109, Final residual = 0.0088413556, No Iterations 3
time step continuity errors : sum local = 6.5470445e-09, global = 2.6478718e-11, cumulative = -0.0021294309
GAMG:  Solving for p_rgh, Initial residual = 0.16833738, Final residual = 0.0056375979, No Iterations 4
time step continuity errors : sum local = 4.4154238e-09, global = -1.4618546e-11, cumulative = -0.0021294309
GAMGPCG:  Solving for p_rgh, Initial residual = 0.27444178, Final residual = 5.3881467e-06, No Iterations 21
time step continuity errors : sum local = 3.3351473e-12, global = -5.6207951e-14, cumulative = -0.0021294309
ExecutionTime = 192.2 s  ClockTime = 193 s

Interface Courant Number mean: 4.1342592e-10 max: 0.0070762091
Courant Number mean: 1.417313e-07 max: 1.9252826
deltaT = 8.191622e-16
Time = 0.0020815663602643

Selected 0 cells for refinement out of 97364.
Selected 25 split points out of a possible 10256.
Unrefined from 97364 to 97189 cells.
Execution time for mesh.update() = 0.29 s
time step continuity errors : sum local = 7.3543405e-09, global = -2.1895987e-14, cumulative = -0.0021294309
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 7.7400147e-05, No Iterations 25
time step continuity errors : sum local = 5.6922704e-13, global = -3.0784214e-14, cumulative = -0.0021294309
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -1.7487261e-23  Max(alpha1) = 1.0000562
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -7.6284817e-23  Max(alpha1) = 1.0000562
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -3.5465846e-23  Max(alpha1) = 1.0000561
GAMG:  Solving for p_rgh, Initial residual = 0.23653799, Final residual = 0.0092471198, No Iterations 4
time step continuity errors : sum local = 1.1919742e-09, global = -1.1890053e-11, cumulative = -0.0021294309
GAMG:  Solving for p_rgh, Initial residual = 0.27559742, Final residual = 0.010717895, No Iterations 5
time step continuity errors : sum local = 1.4469695e-09, global = 9.3509345e-13, cumulative = -0.0021294309
GAMGPCG:  Solving for p_rgh, Initial residual = 0.20976959, Final residual = 5.1609997e-06, No Iterations 21
time step continuity errors : sum local = 1.2498129e-12, global = 3.1235512e-14, cumulative = -0.0021294309
ExecutionTime = 193.42 s  ClockTime = 194 s

Interface Courant Number mean: 1.2819201e-10 max: 0.0018057993
Courant Number mean: 5.0481435e-08 max: 0.29269287
deltaT = 9.8299463e-16
Time = 0.0020815663602653

Selected 25 cells for refinement out of 97189.
Refined from 97189 to 97364 cells.
Selected 0 split points out of a possible 10256.
Execution time for mesh.update() = 0.28 s
time step continuity errors : sum local = 7.5998511e-09, global = 3.7482614e-14, cumulative = -0.0021294309
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 8.5199928e-05, No Iterations 18
time step continuity errors : sum local = 6.4750677e-13, global = -8.2134137e-14, cumulative = -0.0021294309
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -7.2610449e-23  Max(alpha1) = 1.0000561
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -1.3967872e-22  Max(alpha1) = 1.0000561
MULES: Solving for alpha1
Phase-1 volume fraction = 0.89999969  Min(alpha1) = -4.035089e-22  Max(alpha1) = 1.0000561
GAMG:  Solving for p_rgh, Initial residual = 0.17541847, Final residual = 0.0087247483, No Iterations 3
time step continuity errors : sum local = 6.4492676e-09, global = 2.6064436e-11, cumulative = -0.0021294309
GAMG:  Solving for p_rgh, Initial residual = 0.16630811, Final residual = 0.0055667451, No Iterations 4
time step continuity errors : sum local = 4.3435881e-09, global = -1.4441328e-11, cumulative = -0.0021294309
GAMGPCG:  Solving for p_rgh, Initial residual = 0.26960764, Final residual = 5.2469098e-06, No Iterations 21
time step continuity errors : sum local = 3.252421e-12, global = -5.7942249e-14, cumulative = -0.0021294309
ExecutionTime = 194.51 s  ClockTime = 195 s

Yours,

Christian Butcher


Edit: The case is decomposed for 4 processors, in case that is relevant for the fact that some of the Courant numbers are in fact becoming higher, whilst others sit around 0.2. I assume this might be different processor regions. They are split using simple (2 2 1) and the mesh is cylindrical with the z axis as the axis of symmetry. The refinement is occurring at one end of the cylinder, for the higher z values (and should be at least approximately axisymmetric)

[Artur]

Hi Christian,

I am currently experiencing similar problems with interPhaseChangeFoam in 2.2.2. Have you found the solution to your problem and want to share?

Peace,

Artur

[chrisb2244]

Dear Artur,

I'm afraid I don't remember now - I stopped working on that simulation at the end of August last year, and since then have been trying to set up and run other simulations/set-ups.

Some of my current work is using interDyMFoam, or similar, and these do not suffer from the same problem (OF-2.2.2, 2.3.0) so I can't really guess what might be the problem.

Sorry for not having a useful response but thought this might be better than leaving you waiting here indefinitely.

Good luck and best wishes,

Christian

[gautami]

Dear All,

I am using OpenFOAM version 2.2.0.

I was able to run the tutorial "floating Object" which is in 3D in interDyMFoam without any issues.

I am trying to run the same case in 2D.

The following changes have been made, to run the case in 2D

1.In the blockMeshDict file the y coordinates range is changed.
In 3D y axis range was from 0 to 1 and for 2D y axis range is changed as 0 to 0.1.
For the front and back faces,created patch frontandback as empty

2.In the topoSetDict file:
In 3D, box (0.35 0.35 0.44) (0.65 0.65 0.56) was changed in 2D as box (0.35 0.0 0.44) (0.65 0.1 0.56)

3.In the pointDisplacement file:

In 3D

floatingObject
{
type sixDoFRigidBodyDisplacement;
centreOfMass (0.5 0.5 0.5);
momentOfInertia (0.08622222 0.08622222 0.144);
mass 9.6;
rhoInf 1; // needed only for solvers solving for kinematic pressure
report on;
value uniform (0 0 0);
}

In 2D was changed as:

floatingObject
{
type sixDoFRigidBodyDisplacement;
centreOfMass (0.5 0.05 0.5);
momentOfInertia (0.01952 0.08352 0.08);
mass 9.6;
rhoInf 1; // needed only for solvers solving for kinematic pressure
report on;
value uniform (0 0 0);
}
NOTE: Mass not changed

With a change in mass,the changes are as follows:

floatingObject
{
type sixDoFRigidBodyDisplacement;
centreOfMass (0.5 0.05 0.5);
momentOfInertia (0.08133 0.0348 0.033);
mass 4.0;
rhoInf 1; // needed only for solvers solving for kinematic pressure
report on;
value uniform (0 0 0);
}


4.For all the boundary condition files in the folder 0,the patch frontandback is specified as empty.

I run the case and and the run blows with delta T decreasing rapidly.

--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1029
Increased the timePrecision from 12 to 13 to distinguish between timeNames at time 0.0444347581734
GAMG: Solving for p_rgh, Initial residual = 0.00606079623689, Final residual = 5.1801156616e-05, No Iterations 2
time step continuity errors : sum local = 4.69383557933e-05, global = 4.4792886392e-07, cumulative = -0.000506180874993
GAMGPCG: Solving for p_rgh, Initial residual = 1.52101895333e-05, Final residual = 6.35747944249e-09, No Iterations 3
time step continuity errors : sum local = 0.000245342044169, global = -1.53144705592e-06, cumulative = -0.000507712322048
smoothSolver: Solving for epsilon, Initial residual = 0.00132531232957, Final residual = 9.99802871216e-07, No Iterations 41
bounding epsilon, min: -1.04792539095e+26 max: 8.15921076625e+28 average: 3.4517817453e+25
smoothSolver: Solving for k, Initial residual = 0.910004972962, Final residual = 0.0190942329344, No Iterations 1000
ExecutionTime = 120.48 s ClockTime = 121 s

Interface Courant Number mean: 0.00473017730727 max: 14.3965768911
Courant Number mean: 2.42817054949 max: 10937.5681779
deltaT = 6.09445041823e-18
Time = 0.04443475817338


Kindly can somebody please let me know where I am making a mistake and what I could do, to run the case in 2D successfully.Thank you.
Regards,
Gautami.