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

nagetive time marching, ggi patch, rotating frame.

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 9, 2020, 19:29
Default nagetive time marching, ggi patch, rotating frame.
  #1
Member
 
Bidesh Sengupta
Join Date: Sep 2018
Posts: 67
Rep Power: 3
BSengupta is on a distinguished road
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.
Code:
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.
Attached Images
File Type: bmp Figures.bmp (176.9 KB, 2 views)
BSengupta is offline   Reply With Quote

Old   February 11, 2020, 13:17
Default
  #2
HPE
Senior Member
 
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 619
Rep Power: 7
HPE is on a distinguished road
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?
HPE is offline   Reply With Quote

Old   February 12, 2020, 02:18
Default
  #3
Member
 
Bidesh Sengupta
Join Date: Sep 2018
Posts: 67
Rep Power: 3
BSengupta is on a distinguished road
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.


Code:
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.
Attached Images
File Type: jpg pressure.jpg (34.7 KB, 2 views)
File Type: jpg Velocity.jpg (34.5 KB, 1 views)
BSengupta is offline   Reply With Quote

Reply

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
laplacianFoam with source term Herwig OpenFOAM Running, Solving & CFD 17 November 19, 2019 13:47
[mesh manipulation] Importing Multiple Meshes thomasnwalshiii OpenFOAM Meshing & Mesh Conversion 18 December 19, 2015 18:57
Compressor Simulation using rhoPimpleDyMFoam Jetfire OpenFOAM Running, Solving & CFD 107 December 9, 2014 13:38
mixerVesselAMI2D's mass is not balancing sharonyue OpenFOAM Running, Solving & CFD 6 June 10, 2013 09:34
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03


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