# Bounding epsilon and convergence

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

 February 13, 2013, 06:06 Bounding epsilon and convergence #1 New Member   Join Date: Jan 2013 Posts: 4 Rep Power: 6 Hello Foamers. I am kind of new in OpenFoam.. I am trying to run SimpleFoam turbulence simulation on a 3D geometry. It is a kind of pipe with flow splitter. The fluid is water (one phase). I got this kind of error: DILUPBiCG: Solving for Ux, Initial residual = 1.96286e-14, Final residual = 1.96286e-14, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 3.56233e-14, Final residual = 3.56233e-14, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.38946e-13, Final residual = 1.38946e-13, No Iterations 0 DICPCG: Solving for p, Initial residual = 1, Final residual = 2.63715e+06, No Iterations 1001 time step continuity errors : sum local = 7.66448e+73, global = 1.60172e+66, cumulative = 1.60172e+66 [3] #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [3] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [3] #2 in "/lib/libc.so.6" [3] #3 Foam::PBiCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" What is the reason for the negative K? Thanks

February 13, 2013, 08:39
#2
New Member

Join Date: Jan 2013
Posts: 4
Rep Power: 6
I believe that the problem is in the mesh.

I have made a CAD model (solidWorks) and then mesh it in salome. Import it in UNV format to OpenFoam
Attached Images
 Capture2.JPG (32.8 KB, 30 views) Capture1.JPG (78.0 KB, 32 views)

 June 5, 2013, 07:29 #3 New Member   shyam prasad Join Date: Mar 2009 Posts: 25 Rep Power: 10 Hi Kirli, Most of us go through this problem of bounding epsilon or omega. Some of things which helped me are the following, which are taken from various threads of CFD-online. check your mesh using the command checkMesh and look for non-orthogonality, if its between a. 0 to 50 full corrected scheme is applicable, b. 50 to 70 limited correction is required, c. 70 to 80 stability possible, accuracy compromised, d. above 80 stability very difficult to attain For case b use the following settings in fvSchemes gradSchemes { default faceLimited leastSquares 0.5; } laplacianSchemes { default Gauss linear limited 0.33; } snGradSchemes { default limited 0.33; } for case c and d I would suggest to revisit your mesh. Romarius, kostnermo and andreachr like this.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post sagnikmazumdar OpenFOAM Running, Solving & CFD 24 March 1, 2015 08:16 tH3f0rC3 OpenFOAM 7 February 23, 2013 06:37 jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24 swahono OpenFOAM 10 November 26, 2010 06:38 Carlos FLUENT 4 August 27, 2007 11:22

All times are GMT -4. The time now is 15:54.