[BSengupta]

Hi Foamers!


I am trying to simulate a rotating wing with a fixed rpm. Schematic of my problem is attached.


When I am trying to run the simulation, it giving negative time.

Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh turboFvMesh
Initializing the GGI interpolator between master/shadow patches: RotorGGI/StatorGGI
Turbomachine Mixer mesh:
    origin: (0 0 0)
    axis  : (0 0 1)
Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
Allocating field rho

Reading field U

Reading/calculating face flux field phi

Creating MRF model

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model laminar
Create Riemann solver

Allocating field rhoU

Allocating field rhoE

Allocating physDeltaT list for RK and Dual-Time Stepping

Creating oldTime fields for dual-time stepping


Starting time loop

Creating ggi check function object
Moving Cell Zone Name: cellRegion0 rpm: 1250
Moving Cell Zone Name: cellRegion1 rpm: 0
Moving Face Zone Name: RotorGGI_faces rpm: 1250

 Time = -0.0008517870533

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding rho, min: -8002.780332 max: 1126.854835 average: 0.2519080825
new min: 0.8002967477
new max: 1.999596239
bounding rhoE, min: -2413391697 max: 338960313.3 average: -24780.12515
new min: 50034.0114
new max: 999650.5291
bounding U max: 16551.47357
bounding h, min: -710.4404189 max: 907465.0083 average: 306656.9902
new min: 10023.24515
bounding p, min: 2675.587422 max: 401158.8261 average: 104114.4592
new min: 10000.03539
rho L2 Residual: 4.806956853
 rho LInf Residual: 6.972972733

My boundary file is:
RotorGGI
{
type ggi;
nFaces 21232;
startFace 36777457;
rotationAxis (0 0 1);
nCopies 30;
shadowPatch StatorGGI;
zone RotorGGI_faces;
bridgeOverlap true;
}
blade
{
type wall;
nFaces 1135638;
startFace 36798689;
}
farfield
{
type patch;
nFaces 60000;
startFace 37934327;
}
StatorGGI
{
type ggi;
nFaces 38736;
startFace 37994327;
rotationAxis (0 0 1);
nCopies 30;
shadowPatch RotorGGI;
zone StatorGGI_faces;
bridgeOverlap true;
}


Have anyone encountered this kind of problem before. If anyone can give me a clue about this, I shall be really helped.

Thanks a lot for your time.

[HPE]

wow, time travel with OpenFOAM! OpenFOAM is very powerful as always.

Could you please provide a test case, so that we can replicate the reported issue?

[BSengupta]

HI,


Many thanks for your reply.



Just few observations that I made.
The case is running well with rhoPimpleDyMFoam.
But I wish to run this case with dbns solver so I tried transonicMRFDyMFoam, this creates all the problem, although I could run a simple case withis solver. (Figure attached.)


I am trying to verify Caradonna Tung case. The solver ran well for very small time and suddenly went to negative time and crashed.

GGI pair (cylr, cyls)
Area: 1.275277793 1.277828438 Diff = -0.002550645763 or 0.1996078414 % 
Flux: 5.127987472 5.377813765 Diff = -0.3272316814 or 6.381288627 %
GGI pair (discrl, discsl)
Area: 5.30890084 5.31952395 Diff = -0.01062311045 or 0.199700397 % 
Flux: 56.36113422 56.36637665 Diff = -0.01146315011 or 0.02033874987 %
GGI pair (discru, discsu)
Area: 5.30890084 5.31952395 Diff = -0.01062311045 or 0.199700397 % 
Flux: 50.45643701 50.43524362 Diff = -0.01267533051 or 0.02512133488 %

 Time = 4.508921598e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
rho L2 Residual: 3.46313128
rho LInf Residual: 4.763521822

BiCGStab:  Solving for epsilon, Initial residual = 0.9999950781, Final residual = 2.212067388e-14, No Iterations 1
BiCGStab:  Solving for k, Initial residual = 0.008327140338, Final residual = 2.044108189e-16, No Iterations 1

    ExecutionTime = 78.21 s

GGI pair (cylr, cyls)
Area: 1.275277793 1.277828438 Diff = -0.002550645763 or 0.1996078414 % 
Flux: 5.131371036 5.386445304 Diff = -0.337659898 or 6.580305646 %
GGI pair (discrl, discsl)
Area: 5.30890084 5.31952395 Diff = -0.01062311045 or 0.199700397 % 
Flux: 56.36112165 56.36636457 Diff = -0.01147948593 or 0.0203677386 %
GGI pair (discru, discsu)
Area: 5.30890084 5.31952395 Diff = -0.01062311045 or 0.199700397 % 
Flux: 50.45640114 50.43519843 Diff = -0.01269114308 or 0.02515269182 %

 Time = -0.0002009053193

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding rho, min: -19.10611878 max: 40.90783651 average: 0.9802035543
new min: 1e-06
new max: 9.999680957
bounding rhoE, min: -7694045.878 max: 17575530.49 average: 217490.5868
new min: 10000
new max: 1000000
bounding U max: 1.105151028e+10
bounding h, min: -693476.5548 max: 9.102065628e+11 average: 59528188.91
new min: 10000
new max: 100000000
rho L2 Residual: 3.456650242
rho LInf Residual: 4.787876414

BiCGStab:  Solving for epsilon, Initial residual = 0.9999999522, Final residual = 4.955930645e-12, No Iterations 1
BiCGStab:  Solving for k, Initial residual = 0.1883215581, Final residual = 6.680483789e-12, No Iterations 1

    ExecutionTime = 88.02 s

GGI pair (cylr, cyls)
Area: 1.275277793 1.277828438 Diff = -0.002550645763 or 0.1996078414 % 
Flux: 5.598183521 5.294696389 Diff = -0.2341575821 or 4.182742156 %
GGI pair (discrl, discsl)
Area: 5.30890084 5.31952395 Diff = -0.01062311045 or 0.199700397 % 
Flux: 56.35205222 56.35510544 Diff = -0.01107282278 or 0.01964936918 %
GGI pair (discru, discsu)
Area: 5.30890084 5.31952395 Diff = -0.01062311045 or 0.199700397 % 
Flux: 50.26947062 50.25606473 Diff = -0.01603900949 or 0.03190606405 %

 Time = -0.006116840835

[bidesh:26460] *** Process received signal ***
[bidesh:26460] Signal: Floating point exception (8)
[bidesh:26460] Signal code:  (-6)
[bidesh:26460] Failing at address: 0x3e80000675c
[bidesh:26460] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f6deb858f20]
[bidesh:26460] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f6deb858e97]
[bidesh:26460] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f6deb858f20]
[bidesh:26460] [ 3] /home/bidesh/foam/bidesh-4.0/lib/linux64GccDPOpt/libGodunovFlux.so(_ZNK4Foam10roeALEFlux12evaluateFluxERdRNS_6VectorIdEES1_ddS3_S3_ddddddS3_dS3_d+0x440)[0x7f6ded86d140]
[bidesh:26460] [ 4] transonicMRFDyMFoam(+0xac444)[0x563f19f38444]
[bidesh:26460] [ 5] transonicMRFDyMFoam(+0x3b68c)[0x563f19ec768c]
[bidesh:26460] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f6deb83bb97]
[bidesh:26460] [ 7] transonicMRFDyMFoam(+0x40aba)[0x563f19eccaba]
[bidesh:26460] *** End of error message ***
[bidesh:26456] *** Process received signal ***
[bidesh:26456] Signal: Floating point exception (8)
[bidesh:26456] Signal code:  (-6)
[bidesh:26456] Failing at address: 0x3e800006758
[bidesh:26456] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f21dbd3cf20]
[bidesh:26456] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f21dbd3ce97]
[bidesh:26456] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f21dbd3cf20]
[bidesh:26456] [ 3] /home/bidesh/foam/bidesh-4.0/lib/linux64GccDPOpt/libGodunovFlux.so(_ZNK4Foam10roeALEFlux12evaluateFluxERdRNS_6VectorIdEES1_ddS3_S3_ddddddS3_dS3_d+0x440)[0x7f21ddd51140]
[bidesh:26456] [ 4] transonicMRFDyMFoam(+0xac444)[0x556f1e0b2444]
[bidesh:26456] [ 5] transonicMRFDyMFoam(+0x3b68c)[0x556f1e04168c]
[bidesh:26456] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f21dbd1fb97]
[bidesh:26456] [ 7] transonicMRFDyMFoam(+0x40aba)[0x556f1e046aba]
[bidesh:26456] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 4 with PID 26456 on node bidesh exited on signal 8 (Floating point exception).

May I have your mail id as I cannot upload my case here. My mail id is bideshsengupta.08@gmail.com


Thanks a lot.