CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

error while running modified pimple solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 28, 2015, 06:59
Default error while running modified pimple solver
  #1
New Member
 
rafath
Join Date: Jun 2014
Location: mumbai
Posts: 24
Rep Power: 4
R_21 is on a distinguished road
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
"
R_21 is offline   Reply With Quote

Reply

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
OF 2.1.1 error in compiling modified solver anuragm OpenFOAM Programming & Development 2 November 28, 2014 10:37
Running a solver from a case directory stored anywhere mich_K OpenFOAM Running, Solving & CFD 3 August 8, 2014 03:57
[ANSYS Meshing] Help with element size sandri_92 ANSYS Meshing & Geometry 13 May 26, 2014 09:05
Error when I try a case with pimple instead of simple solver samiam1000 OpenFOAM 6 July 29, 2013 12:56
Running CFD solver on PC Reaz Hasan Main CFD Forum 16 March 14, 2002 06:55


All times are GMT -4. The time now is 11:14.