Trouble with Mesh nonorthogonality
Hi Foamers,
I prepared my solver which is nonlinear PDE. I am working on suspension flow (Fluid+Particles) through different geometries. When i compiled my solver then it is giving no errors, but when i run my case file by using my solver it is giving some error message. Error message: Starting time loop Time = 1nl [0 0 1 0 0 0 0] 0.02716 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 4.70429e06, No Iterations 5 DILUPBiCG: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0 swak4Foam: Allocating new repository for sampledGlobalVariables #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<Foam::DILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libfiniteVolume.so" #9 Foam::fvMatrix<double>::solve() in "/home/malli_reddy/OpenFOAM/malli_reddy2.1.1/platforms/linuxGccDPOpt/bin/twoFoam" #10 in "/home/malli_reddy/OpenFOAM/malli_reddy2.1.1/platforms/linuxGccDPOpt/bin/twoFoam" #11 __libc_start_main in "/lib/i386linuxgnu/libc.so.6" #12 in "/home/malli_reddy/OpenFOAM/malli_reddy2.1.1/platforms/linuxGccDPOpt/bin/twoFoam" Floating point exception My governing equation involves nonlinear equations. When i give good initial guess (i.e by running the case by using linear solver instead of nonlinear solver at first for the guess) then the case is running properly and giving good results for simple channel flow.(in which for Mesh nonorthogonality Max: 0 average: 0) When i repeat the above same procedure for non zero Mesh nonorthogonality then after two iterations the solution is diverging. Please let me know what is the problem and what above error trying to explain. Thanks Regards Mallikarjuna 
Only Ux and Uy?..You need to consider a Control Volume for OF.It appears there is some problem with the geometry specification.If you are going for solely 2D case, create a unit mesh in the z direction...
hope it helps 
Quote:
If you're running a 2D case OF will only solve for the x and y components in the segregated solution so there's no problem there. If introducing mesh nonorthogonality is killing your stability, try increasing the number of nNonOrthogonalCorrectors in the /system/fvSolution dictionary. If you mesh is too skewed, however, you might be out of luck and want to try some different meshing strategies. 
O my......totally forgot about segregated solver situation. Thank you very much for the correction.

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