
[Sponsors] 
September 30, 2013, 20:52 
Diverging solution in transonicMRFDyMFoam

#1 
New Member
Trent
Join Date: May 2013
Posts: 2
Rep Power: 0 
Sponsored Links
I am attempting to replicate a 2D compressor profile test case using the transonicMRFDyMFoam solver from the Density Based Turbo package for OpenFoam1.6extend. I have run the case in sonicFoam and have a reasonable result, which I am attempting to replicate in transonicMRFDyMFoam. sonicFoam results: sonicFoam_p_converged.png sonicFoam_U_converged.png The pressure and velocity transonicMRFDyMFoam results diverge very quickly. transonicMRFDyMFoam_p.png transonicMRFDyMFoam_U.png In the attached case, I used the keyword "cyclic" for the periodic boundary conditions (to model an infinite array of blades) at the top and bottom (using createPatch to merge two patches into one periodic patch), however I have also used 'cyclicGgi' periodic boundary conditions to the same result. I have been stuck on this problem for a while now, any help would be greatly appreciated. I have attached the set up files, however the full polymesh files were too large to attach, so only the 'boundary' polyMesh file has been included. sc10_tranconicMRFDyMFoam.tar.gz Thanks heaps 

Sponsored Links 
April 29, 2014, 01:46 
Hi

#2 
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
Hi Trent,
Have you been able to run transonicMRFDyMFoam case of your's successfully ? I need some help in this regard. Can you please help me ? 

May 9, 2014, 17:06 

#3 
New Member
Join Date: May 2011
Posts: 28
Rep Power: 8 
hi sam
please list more detail about what you want to know. rencently i am doing rans simulation by using transonicMRFDyMFoam for compressor. 

May 12, 2014, 23:53 
Floating point exception (core dumped) in transonicMRFDyMFoam

#4  
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
Quote:
I have set up the case for radial inflow turbine passage. I have used cyclic boundary condition for periodic surfaces and overlap GGI between the rotating and stationary parts. When I run the case I am getting following error. Code:
sml5kor@BMH301562:~/OpenFOAM/sml5kor1.6ext/run/Turbine_CASE_29.04.2014$ transonicMRFDyMFoam /**\  =========    \\ / F ield  OpenFOAM Extend Project: Open source CFD   \\ / O peration  Version: 1.6ext   \\ / A nd  Web: www.extendproject.de   \\/ M anipulation   \**/ Build : 1.6ext Exec : transonicMRFDyMFoam Date : Apr 30 2014 Time : 10:02:43 Host : BMH301562 PID : 28172 Case : /home/sml5kor/OpenFOAM/sml5kor1.6ext/run/Turbine_CASE_29.04.2014 nProcs : 1 SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create dynamic mesh for time = 0 Selecting dynamicFvMesh staticFvMesh 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 kOmegaSST kOmegaSSTCoeffs { Prt 0.9; alphaK1 0.85034; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.85616; gamma1 0.5532; gamma2 0.4403; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; c1 10; } Create Riemann solver Allocating field rhoU Allocating field rhoE Allocating physDeltaT list for RK and DualTime Stepping Starting time loop Time = 5e05 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: 5.22287757622644 rho LInf Residual: 6.44791071235071 DILUPBiCG: Solving for omega, Initial residual = 0.00130432164252392, Final residual = 2.96795116630587e09, No Iterations 1 bounding omega, min: 8706879.03922889 max: 81910182030.3853 average: 143455780.414936 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 6.44126182533454e11, No Iterations 2 bounding k, min: 0.10155708303169 max: 3128.93273261909 average: 5.22430823460557 ExecutionTime = 22.19 s performance data: Head (m) = 153590.297255837 TOmega (W) = 105.444029188369 Eff (%) = 9.05754907371286 Forces = (0.116067242676771 0.448441765109099 1.03710250567382) Moments = (0.0109810410962604 0.0145957547140102 0.00728641017051233) Fluid power output: dEm (W) = 1164.15631127401 Head (m) = 153590.297255837 forces output: forces(pressure, viscous)((0.0246120695715455 0.0593173620106552 0.064901448211144) (0.0914551731052256 0.389124403098444 0.972201057462674)) moment(pressure, viscous)((0.00701764721585918 3.55306881895442e05 0.000938169101871696) (0.00396339388040122 0.0146312854021998 0.00634824106864064)) Time = 0.0001 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: 3571.52349710333 max: 135.300557117857 average: 0.0624592646242397 bounding rhoE, min: 5840210714.15816 max: 285770183.417902 average: 16821.1460355255 bounding U max: 9163515.75862727 bounding h, min: 559187.987560337 max: 1329471.19462314 average: 570184.378251595 rho L2 Residual: 5.75472688663313 rho LInf Residual: 7.10596896265805 DILUPBiCG: Solving for omega, Initial residual = 0.000870710888200336, Final residual = 4.13051427693362e09, No Iterations 1 bounding omega, min: 485017.969038368 max: 82333537116.776 average: 144331930.942993 DILUPBiCG: Solving for k, Initial residual = 0.740046810060922, Final residual = 7.44167813561397e09, No Iterations 3 ExecutionTime = 32.92 s performance data: Head (m) = 79014.871433648 TOmega (W) = 113.872177860226 Eff (%) = 0.263868113135903 Forces = (0.132477788867598 0.487273650008056 1.06569241218417) Moments = (0.011384295242643 0.0150292348287509 0.00786881344809862) Fluid power output: dEm (W) = 43154.9596906303 Head (m) = 79014.871433648 forces output: forces(pressure, viscous)((0.0453001853033088 0.115529386787037 0.125766544339912) (0.0871776035642893 0.371744263221019 0.939925867844254)) moment(pressure, viscous)((0.00758893570416305 0.000882449356738547 0.00179562705778156) (0.00379535953847992 0.0141467854720123 0.00607318639031707)) Time = 0.00015 Floating point exception (core dumped) Could you please help me in this regard. Regards, Sam 

May 14, 2014, 04:09 

#5 
New Member
Join Date: May 2011
Posts: 28
Rep Power: 8 
Hi, Sam
1. I have a look at the attached files. For the inlet boundary, total pressure, total temperature and flow direction are usually specified. I find that you did not set the inlet flow direction, but zeroGradient in 0/U file is set in the inlet patch instead. I think you can try the typical boundary condition provided by Borm. 2. According to my experience, you can first turn off turbulence, set a higher outlet pressure (e.g. 190000Pa) and set the rotating speed to zero. If this calculation works, then you can gradually improve the rotating speed and lower the outlet pressure. 3. For the interrow treatment, which kind of method did you use? (mixing plane or frozen rotor) Last edited by dowlee; May 14, 2014 at 04:12. Reason: wrong words 

May 15, 2014, 02:01 

#6  
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
Quote:
1. Can you provide me the typical boundary condition provided by Borm. ? Means where I can find them. 3. frozen rotor. means MRF method as Borm used in transonicMRFDyMFoam solver . 

May 15, 2014, 04:45 

#7 
New Member
Join Date: May 2011
Posts: 28
Rep Power: 8 
Hi, you will download the folder named densityBasedSolver in http://sourceforge.net/projects/openfoamextend/.
There are some tutorials provided by Brom in the directory "densityBasedSolver/run" 

May 17, 2014, 05:14 

#8 
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
I was able to run the case but after 1000 iterations I am getting wrong flow rate .
Code:
Time = 1000 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 U max: 46939.3687296215 rho L2 Residual: 5.19492515331512 rho LInf Residual: 6.6959350398147 DILUPBiCG: Solving for omega, Initial residual = 3.95477114014134e06, Final residual = 2.94743732008884e11, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.000135496572893564, Final residual = 4.27530898304972e09, No Iterations 1 bounding k, min: 0.00138428982961019 max: 258480.340654354 average: 1085.527034435 ExecutionTime = 7105.86 s MassFlows: INLET = 2.69804743578668e09 MassFlows: OUTLET = 0.123375043979185 Averages of p : INLET = 200000 OUTLET = 100203.673206441 Averages of rho : INLET = 0.700070745782449 OUTLET = 0.696048056366656 Averages of T : INLET = 992.999999999999 OUTLET = 500.385591367588 Averages of U : INLET = (1.50903942957102e06 0 0) OUTLET = (4.30635301362846 4.48022752259606 98.9735954313029) end Pressure and Temperature are fine . What could be the reason for this ? 

May 19, 2014, 01:28 
BC and Constant files

#9 
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
Hi,
Please find the attached BC and Initial values of the case which has run but wrong mass flow. Let me know your feed back upon this . Regards, Sam 

May 19, 2014, 04:25 
"libsimpleFunctionObjects.so"

#10 
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
Code:
functions ( inletMassFlow { type patchMassFlow; functionObjectLibs ( "libsimpleFunctionObjects.so" ); verbose true; patches ( inlet_0 ); factor 23; } Whats the meaning of "factor 23;". How to calculate this ? I feel this factor only making some problem for my simulation . 

May 19, 2014, 05:19 

#11  
New Member
Join Date: May 2011
Posts: 28
Rep Power: 8 
Quote:


May 19, 2014, 05:28 
'outputControlMode' not found in inletMassFlow

#12 
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
Code:
Starting time loop > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 108 'outputControlMode' not found in inletMassFlow Assuming: timestep > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 116 'outputInterval' not found in inletMassFlow Assuming: 1 phi: phi Compressible: 1 Turbulent: 1 LES: 0 > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 108 'outputControlMode' not found in outletMassFlow Assuming: timestep > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 116 'outputInterval' not found in outletMassFlow Assuming: 1 phi: phi Compressible: 1 Turbulent: 1 LES: 0 > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 108 'outputControlMode' not found in areaAverage Assuming: timestep > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 116 'outputInterval' not found in areaAverage Assuming: 1 Time = 1 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: 6.44313258180784 rho LInf Residual: 7.65994499810326 DILUPBiCG: Solving for omega, Initial residual = 0.00244809190957999, Final residual = 2.65944354682262e09, No Iterations 1 bounding omega, min: 24070579.7366927 max: 182073584542.249 average: 319780497.816931 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 5.11096081209236e11, No Iterations 3 ExecutionTime = 20.82 s MassFlows: INLET = 1.76248019141362 MassFlows: OUTLET = 0.0595373845935599 Averages of p : INLET = 200000 OUTLET = 100000.341228185 Averages of rho : INLET = 0.710387031990251 OUTLET = 0.695171381547665 Averages of T : INLET = 978.579590078318 OUTLET = 500.000892595263 Averages of U : INLET = (1.50903942957102e06 0 0) OUTLET = (4.99964468397136 4.99890645781862 99.9962223316277) Time = 2 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: 6.09112121652158 rho LInf Residual: 7.21831669034706 DILUPBiCG: Solving for omega, Initial residual = 0.000171867025540015, Final residual = 6.85362954444674e10, No Iterations 1 bounding omega, min: 8373293.02948838 max: 182228019464.836 average: 321729239.314316 DILUPBiCG: Solving for k, Initial residual = 0.0497553547344649, Final residual = 1.76120797453748e10, No Iterations 2 ExecutionTime = 28.15 s MassFlows: INLET = 1.30938549343049e09 MassFlows: OUTLET = 0.0595352322647667 Averages of p : INLET = 200000 OUTLET = 100000.68756607 Averages of rho : INLET = 0.700070745782449 OUTLET = 0.695172520026074 Averages of T : INLET = 992.999999999999 OUTLET = 500.001805227856 Averages of U : INLET = (1.50903942957102e06 0 0) OUTLET = (4.99910923914695 4.99772011195967 99.9895354653937) Time = 3 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: 6.04557557426488 rho LInf Residual: 7.1978292209289 DILUPBiCG: Solving for omega, Initial residual = 0.000170850205219554, Final residual = 5.30721748237609e10, No Iterations 1 bounding omega, min: 766813.296739794 max: 182364930358.078 average: 323558740.303541 DILUPBiCG: Solving for k, Initial residual = 0.021230682776069, Final residual = 7.83053164320755e11, No Iterations 2 ExecutionTime = 35.19 s MassFlows: INLET = 1.29037051276754e09 MassFlows: OUTLET = 0.059531348423872 Averages of p : INLET = 200000 OUTLET = 100001.033089627 Averages of rho : INLET = 0.700070745782449 OUTLET = 0.6951736690745 Averages of T : INLET = 992.999999999999 OUTLET = 500.002706121094 Averages of U : INLET = (1.50903942957102e06 0 0) OUTLET = (4.99857272900048 4.99659077737812 99.9830652506599) 

May 19, 2014, 05:31 

#13 
New Member
Join Date: May 2011
Posts: 28
Rep Power: 8 
Hi, I am sorry that I am very busy recently to reply to you late. There seem to be no problems for your case. for the outlet pressure boundary, you can use "inletOutlet"instead of "zeroGradient" to avoid reverse flows.


May 19, 2014, 05:35 

#14 
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
Hi,
Really that's not a problem. You are replying and guiding me so i am very thankful to u. Thank you very much. Can you share your personnel mail ID if you don't mind ? 

May 19, 2014, 05:38 

#15  
New Member
Join Date: May 2011
Posts: 28
Rep Power: 8 
Quote:


May 19, 2014, 06:45 

#16  
New Member
Join Date: May 2011
Posts: 28
Rep Power: 8 
Quote:


May 19, 2014, 06:49 

#17  
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
Quote:
I will contact you through this mail ID. 

May 19, 2014, 09:57 
mass flow Out of range

#18 
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
Hi ,
I tried with inletOutlet BC at OUTLET patch and laminar case but the output of mass flow remains the same .. which is as shown below. Code:
sml5kor@BMH301562:~/OpenFOAM/sml5kor1.6ext/run/Turbine_CASE_17.05.2014$ transonicMRFDyMFoam /**\  =========    \\ / F ield  OpenFOAM Extend Project: Open source CFD   \\ / O peration  Version: 1.6ext   \\ / A nd  Web: www.extendproject.de   \\/ M anipulation   \**/ Build : 1.6ext Exec : transonicMRFDyMFoam Date : May 19 2014 Time : 19:24:29 Host : BMH301562 PID : 26940 Case : /home/sml5kor/OpenFOAM/sml5kor1.6ext/run/Turbine_CASE_17.05.2014 nProcs : 1 SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create dynamic mesh for time = 0 Selecting dynamicFvMesh staticFvMesh Reading thermophysical properties Selecting thermodynamics package hPsiThermo<pureMixture<sutherlandTransport<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 laminar Create Riemann solver Allocating field rhoU Allocating field rhoE Allocating physDeltaT list for RK and DualTime Stepping Starting time loop > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 108 'outputControlMode' not found in inletMassFlow Assuming: timestep > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 116 'outputInterval' not found in inletMassFlow Assuming: 1 > FOAM Warning : From function SolverInfo::SolverInfo(const dictionary& dict,const objectRegistry &obr) in file SolverInfo/SolverInfo.C at line 87 Neither LES nor RAS found. Assuming no turbulence phi: phi Compressible: 1 Turbulent: 0 LES: 0 > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 108 'outputControlMode' not found in outletMassFlow Assuming: timestep > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 116 'outputInterval' not found in outletMassFlow Assuming: 1 > FOAM Warning : From function SolverInfo::SolverInfo(const dictionary& dict,const objectRegistry &obr) in file SolverInfo/SolverInfo.C at line 87 Neither LES nor RAS found. Assuming no turbulence phi: phi Compressible: 1 Turbulent: 0 LES: 0 > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 108 'outputControlMode' not found in areaAverage Assuming: timestep > FOAM Warning : From function simpleFunctionObject::simpleFunctionObject in file simpleFunctionObject/simpleFunctionObject.C at line 116 'outputInterval' not found in areaAverage Assuming: 1 Time = 1 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: 6.40372710604553 rho LInf Residual: 7.49081222690984 ExecutionTime = 14.52 s MassFlows: INLET = 1.76248019141362 MassFlows: OUTLET = 0 Averages of p : INLET = 200000 OUTLET = 100000 Averages of rho : INLET = 0.710387031990251 OUTLET = 0.695170250561966 Averages of T : INLET = 978.579590078318 OUTLET = 500.000000000001 Averages of U : INLET = (1.50903942957102e06 0 0) OUTLET = (9.31784332147014e18 7.85744018292775e18 1.7403312614265e17) Time = 2 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: 5.83719515087281 rho LInf Residual: 7.11130368620766 ExecutionTime = 18.25 s MassFlows: INLET = 1.30938549343049e09 MassFlows: OUTLET = 2.58934377791194e20 Averages of p : INLET = 200000 OUTLET = 100000 Averages of rho : INLET = 0.700070745782449 OUTLET = 0.695170250561966 Averages of T : INLET = 992.999999999999 OUTLET = 500.000000000001 Averages of U : INLET = (1.50903942957102e06 0 0) OUTLET = (1.17465609866511e17 3.25777301233301e18 3.34564026928883e16) Time = 3 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: 5.71596933487263 rho LInf Residual: 7.0872781206772 ExecutionTime = 22.05 s MassFlows: INLET = 1.29037051276754e09 MassFlows: OUTLET = 2.40670548093986e19 Averages of p : INLET = 200000 OUTLET = 100000 Averages of rho : INLET = 0.700070745782449 OUTLET = 0.695170250561966 Averages of T : INLET = 992.999999999999 OUTLET = 500.000000000001 Averages of U : INLET = (1.50903942957102e06 0 0) OUTLET = (7.42904399900281e18 8.66558587599736e18 1.96335447031351e15) 

May 21, 2014, 07:36 
mass flow problem

#19 
Senior Member
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 5 
Hi,
Now the case is running but gives some error after some 100's of iterations. The error is as follows Code:
Time = 2000 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: 1057.88 max: 3688.91 average: 0.756053 bounding rhoE, min: 1.08367e+09 max: 1.20384e+10 average: 385827 bounding U max: 9.60207e+07 bounding h, min: 3.98146e+07 max: 7.099e+07 average: 620124 From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/sml5kor/OpenFOAM/OpenFOAM1.6ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73 Maximum number of iterations exceeded. Rescue by HJ rho L2 Residual: 6.93647 rho LInf Residual: 9.05835 DILUPBiCG: Solving for omega, Initial residual = 3.88886e06, Final residual = 7.15152e09, No Iterations 1 bounding omega, min: 3.69309e+09 max: 1.75975e+11 average: 3.65849e+08 DILUPBiCG: Solving for k, Initial residual = 0.00153104, Final residual = 2.9086e10, No Iterations 3 bounding k, min: 61245.7 max: 3.02127e+07 average: 22009.7 ExecutionTime = 7454.68 s MassFlows: INLET = 0.123706 MassFlows: OUTLET = 0.0734246 Averages of p : INLET = 175946 OUTLET = 105000 Averages of rho : INLET = 0.636106 OUTLET = 0.847162 Averages of T : INLET = 961.373 OUTLET = 495.48 Averages of U : INLET = (159.772 0 0) OUTLET = (27.5507 9.71856 17.1995) end Inlet pressure 2 bar Inlet temperature 973 K What would be the problem ? Please help me . Regards, Sam 

May 22, 2014, 03:40 

#20  
New Member
Join Date: May 2011
Posts: 28
Rep Power: 8 
Quote:


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
underrelaxation factors > level of residuals  Zigainer  FLUENT  19  July 21, 2017 16:53 
grid dependancy  gueynard a.  Main CFD Forum  19  June 27, 2014 21:22 
diverging solution in bubble rise  shash  FLUENT  0  August 11, 2012 16:39 
IcoFoam parallel woes  msrinath80  OpenFOAM Running, Solving & CFD  9  July 22, 2007 02:58 
Wall functions  Abhijit Tilak  Main CFD Forum  6  February 5, 1999 02:16 
Sponsored Links 