error while running modified pimple solver
 

Hi,
Can someone Please help me with this Error. No idea why my courant number is increasing at this rate. I had modified PimpleFoam by adding these two terms , solving for these and then calculating the force term which is then added into the UEqn

solve
(
fvm::laplacian(epsilonR, phiE)
);
solve
(
fvm::laplacian(epsilonR,rhoC) == (rhoC/sqr(lamdaD))
);
force = rhoC * -fvc::grad(phiE) * phiEmax * rhoCmax / rho;

ERROR:

"PIMPLE: no residual control data found. Calculations will employ 2 corrector loops


Starting time loop

forces forces:

forceCoeffs forceCoeffs:

Courant Number mean: 0.024509014 max: 0.026000006
Time = 0.001

PIMPLE: iteration 1
HeSre

GAMG: Solving for phiE, Initial residual = 1, Final residual = 0.0038661331, No Iterations 1
phiE

GAMG: Solving for rhoC, Initial residual = 1, Final residual = 0.0071713282, No Iterations 1
problem

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 3.6621659e-05, No Iterations 1
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0055900888, No Iterations 5
time step continuity errors : sum local = 2.276262e-09, global = -4.5947473e-10, cumulative = -4.5947473e-10
PIMPLE: iteration 2
HeSre

GAMG: Solving for phiE, Initial residual = 0.0011334459, Final residual = 5.3259934e-08, No Iterations 8
phiE

GAMG: Solving for rhoC, Initial residual = 0.99997241, Final residual = 4.0807083e-08, No Iterations 18
problem

DILUPBiCG: Solving for Ux, Initial residual = 0.99999588, Final residual = 1.984607e-10, No Iterations 2
GAMG: Solving for p, Initial residual = 0.99821094, Final residual = 9.1749927e-08, No Iterations 21
time step continuity errors : sum local = 4.1502946e-10, global = -1.6314663e-10, cumulative = -6.2262136e-10
ExecutionTime = 6.03 s ClockTime = 6 s

forces forces output:
forces(pressure,viscous,porous) = ((5629322.9 -4.3812668e+09 0),(-135427.57 -69123.817 0),(0 0 0))
moment(pressure,viscous,porous) = ((2.1906334e+09 2814661.4 -3.7213028e+08),(34561.909 -67713.784 7125671.4),(0 0 0))

forceCoeffs forceCoeffs output:
Cm = -71569.532
Cd = 10772.344
Cl = -8590854.7
Cl(f) = -4366996.9
Cl(r) = -4223857.8

Courant Number mean: 0.4743086 max: 35.507646
Time = 0.002

PIMPLE: iteration 1
HeSre

GAMG: Solving for phiE, Initial residual = 5.3727614e-08, Final residual = 5.3727614e-08, No Iterations 0
phiE

GAMG: Solving for rhoC, Initial residual = 0.99996636, Final residual = 0.0031849954, No Iterations 6
problem

DILUPBiCG: Solving for Ux, Initial residual = 0.99999946, Final residual = 0.054384783, No Iterations 1
GAMG: Solving for p, Initial residual = 0.99966943, Final residual = 0.0087385975, No Iterations 7
time step continuity errors : sum local = 28.348433, global = 0.048852483, cumulative = 0.048852482
PIMPLE: iteration 2
HeSre

GAMG: Solving for phiE, Initial residual = 5.3727614e-08, Final residual = 5.3727614e-08, No Iterations 0
phiE

GAMG: Solving for rhoC, Initial residual = 0.9999664, Final residual = 4.1562921e-08, No Iterations 18
problem

DILUPBiCG: Solving for Ux, Initial residual = 0.99637735, Final residual = 4.3410571e-07, No Iterations 169
GAMG: Solving for p, Initial residual = 0.91262788, Final residual = 8.6100457e-08, No Iterations 30
time step continuity errors : sum local = 1.3939272, global = 0.43977182, cumulative = 0.4886243
ExecutionTime = 13.8 s ClockTime = 14 s

forces forces output:
forces(pressure,viscous,porous) = ((3.4949523e+22 -2.1121571e+25 0),(-1.2209863e+15 -7.832371e+13 0),(0 0 0))
moment(pressure,viscous,porous) = ((1.0560785e+25 1.7474761e+22 -5.606328e+24),(3.9161855e+13 -6.1049315e+14 6.1527873e+16),(0 0 0))

forceCoeffs forceCoeffs output:
Cm = -1.09928e+21
Cd = 6.8528474e+19
Cl = -4.1414844e+22
Cl(f) = -2.1806702e+22
Cl(r) = -1.9608142e+22

Courant Number mean: 4.6288647e+09 max: 6.4880419e+10
Time = 0.003

PIMPLE: iteration 1
HeSre

GAMG: Solving for phiE, Initial residual = 5.3727614e-08, Final residual = 5.3727614e-08, No Iterations 0
phiE

GAMG: Solving for rhoC, Initial residual = 0.99996645, Final residual = 0.003103529, No Iterations 6
problem

DILUPBiCG: Solving for Ux, Initial residual = 0.26573861, Final residual = 0.016023816, No Iterations 1
GAMG: Solving for p, Initial residual = 6.8243353e-08, Final residual = 6.8243353e-08, No Iterations 0
time step continuity errors : sum local = 2.5742666e+08, global = 808722.85, cumulative = 808723.34
PIMPLE: iteration 2
HeSre

GAMG: Solving for phiE, Initial residual = 5.3727614e-08, Final residual = 5.3727614e-08, No Iterations 0
phiE

GAMG: Solving for rhoC, Initial residual = 0.99996648, Final residual = 4.8058389e-08, No Iterations 18
problem

DILUPBiCG: Solving for Ux, Initial residual = 0.9999915, Final residual = 5.9945125e-09, No Iterations 3
GAMG: Solving for p, Initial residual = 0.28940382, Final residual = 8.8857892e-08, No Iterations 20
time step continuity errors : sum local = 6.7338046e+08, global = 26600070, cumulative = 27408793
ExecutionTime = 18.05 s ClockTime = 18 s

forces forces output:
forces(pressure,viscous,porous) = ((5.1864076e+33 -1.3159236e+37 0),(-1.4750133e+23 -8.950878e+21 0),(0 0 0))
moment(pressure,viscous,porous) = ((6.5796178e+36 2.5932038e+33 -5.713251e+36),(4.475439e+21 -7.3750664e+22 7.4316262e+24),(0 0 0))

forceCoeffs forceCoeffs output:
Cm = -1.1202453e+33
Cd = 1.0169427e+31
Cl = -2.5802423e+34
Cl(f) = -1.4021457e+34
Cl(r) = -1.1780966e+34

Courant Number mean: 3.893362e+17 max: 1.5326777e+19
Time = 0.004

PIMPLE: iteration 1
HeSre

GAMG: Solving for phiE, Initial residual = 5.3727614e-08, Final residual = 5.3727614e-08, No Iterations 0
phiE

GAMG: Solving for rhoC, Initial residual = 0.9999665, Final residual = 0.0031738288, No Iterations 6
problem

DILUPBiCG: Solving for Ux, Initial residual = 0.079714094, Final residual = 0.0051008615, No Iterations 1
GAMG: Solving for p, Initial residual = 1.9584815e-12, Final residual = 1.9584815e-12, No Iterations 0
time step continuity errors : sum local = 6.2253282e+15, global = 8.1698749e+13, cumulative = 8.1698776e+13
PIMPLE: iteration 2
HeSre

GAMG: Solving for phiE, Initial residual = 5.3727614e-08, Final residual = 5.3727614e-08, No Iterations 0
phiE

GAMG: Solving for rhoC, Initial residual = 0.99996651, Final residual = 4.7719294e-08, No Iterations 18
problem

DILUPBiCG: Solving for Ux, Initial residual = 0.99992586, Final residual = 1.0915491e-08, No Iterations 2
GAMG: Solving for p, Initial residual = 4.5868743e-07, Final residual = 7.0974347e-08, No Iterations 3
time step continuity errors : sum local = 3.2206052e+20, global = -1.0586468e+19, cumulative = -1.0586387e+19
ExecutionTime = 21.85 s ClockTime = 22 s

forces forces output:
forces(pressure,viscous,porous) = ((3.6172817e+46 -3.6394692e+49 0),(-3.3877048e+28 -4.9447829e+25 0),(0 0 0))
moment(pressure,viscous,porous) = ((1.8197346e+49 1.8086409e+46 1.0197544e+50),(2.4723915e+25 -1.6938524e+28 1.7017471e+30),(0 0 0))

forceCoeffs forceCoeffs output:
Cm = 1.9995185e+46
Cd = 7.0927093e+43
Cl = -7.1362142e+46
Cl(f) = -1.5685886e+46
Cl(r) = -5.5676256e+46

Courant Number mean: 1.2891281e+23 max: 1.403584e+25
Time = 0.005

PIMPLE: iteration 1
HeSre

GAMG: Solving for phiE, Initial residual = 5.3727614e-08, Final residual = 5.3727614e-08, No Iterations 0
phiE

GAMG: Solving for rhoC, Initial residual = 0.99996646, Final residual = 0.0045614595, No Iterations 5
problem

DILUPBiCG: Solving for Ux, Initial residual = 0.11153223, Final residual = 0.0037393102, No Iterations 1
GAMG: Solving for p, Initial residual = 1.2854517e-19, Final residual = 1.2854517e-19, No Iterations 0
time step continuity errors : sum local = 3.7452618e+21, global = -6.9557705e+13, cumulative = -1.0586456e+19
PIMPLE: iteration 2
HeSre

GAMG: Solving for phiE, Initial residual = 5.3727614e-08, Final residual = 5.3727614e-08, No Iterations 0
phiE

GAMG: Solving for rhoC, Initial residual = 0.99996638, Final residual = 6.792861e-08, No Iterations 16
problem

DILUPBiCG: Solving for Ux, Initial residual = 0.99999904, Final residual = 6.0169959e-09, No Iterations 2
GAMG: Solving for p, Initial residual = 1.8666026e-12, Final residual = 1.8666026e-12, No Iterations 0
time step continuity errors : sum local = 7.7637328e+28, global = -3.2681334e+15, cumulative = -1.0589724e+19
ExecutionTime = 25.45 s ClockTime = 26 s

forces forces output:
forces(pressure,viscous,porous) = ((3.6172817e+46 -3.6394692e+49 0),(-5.3520264e+36 1.1110674e+34 0),(0 0 0))
moment(pressure,viscous,porous) = ((1.8197346e+49 1.8086409e+46 1.0197544e+50),(-5.5553372e+33 -2.6760132e+36 2.6818729e+38),(0 0 0))

forceCoeffs forceCoeffs output:
Cm = 1.9995185e+46
Cd = 7.0927093e+43
Cl = -7.1362142e+46
Cl(f) = -1.5685886e+46
Cl(r) = -5.5676256e+46

Courant Number mean: 7.574677e+30 max: 1.3363822e+32
"