|
[Sponsors] |
Problems running customized solver for diffusion |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 30, 2012, 07:49 |
Problems running customized solver for diffusion
|
#1 |
Member
M Mallikarjuna Reddy
Join Date: Jul 2012
Posts: 91
Rep Power: 14 |
Hi everyone,
I am fresher to OpenFOAM, and i developed my own solver. When i use the wmake command there is no error, (i think there is no error in the solver). But when i used run my case using solver the error message coming like as follows, error message:- Create time Create mesh for time = 0 Reading transportProperties Reading field p Reading field U Reading field T Reading/calculating face flux field phi No field sources present SIMPLE: convergence criteria field p tolerance 0.01 field U tolerance 0.001 field T tolerance 1e-05 Starting time loop Time = 1 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 4.00548e-06, No Iterations 4 DILUPBiCG: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0 #0 Foam::error:rintStack(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:ILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #4 Foam:ILUPreconditioner:ILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #5 Foam::lduMatrix:reconditioner::addasymMatrixCons tructorToTable<Foam:ILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #6 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #7 Foam:BiCG::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_reddy-2.1.1/platforms/linuxGccDPOpt/bin/myDiffusiveFoam" #10 in "/home/malli_reddy/OpenFOAM/malli_reddy-2.1.1/platforms/linuxGccDPOpt/bin/myDiffusiveFoam" #11 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6" #12 in "/home/malli_reddy/OpenFOAM/malli_reddy-2.1.1/platforms/linuxGccDPOpt/bin/myDiffusiveFoam" Floating point exception malli_reddy@ubuntu:~/OpenFOAM/malli_reddy-2.1.1/PROJECT/myDiffusive$ can anyone suggest me what is the error and how to rectify it. thank you. |
|
August 30, 2012, 15:01 |
about error..
|
#2 |
Member
M Mallikarjuna Reddy
Join Date: Jul 2012
Posts: 91
Rep Power: 14 |
Hi everyone,
I am fresher to OpenFOAM, and i developed my own solver for predicting velocity profile and particle concentration profile for suspension flow (for steady state). And in my momentum equation, the viscosity is not constant, it is function of particle volume fraction. And i written the momentum equation as follows:tmp<fvVectorMatrix> UEqn ( fvm::div(phi, U) - fvm::laplacian(nu*pow((1-T/0.68),-1.82), U) ); UEqn().relax(); sources.constrain(UEqn()); solve(UEqn() == -fvc::grad(p)); where T is the particle volume fraction (dimensionless scalar field). For predicting particle volume fraction, the governing equation written as: solve ( fvm::div(phi, T) - fvm::laplacian(0.62*pow(a,2)*mag(symm(fvc::grad(U) ))*pow(T,2)*pow(0.68-T,-1), T) ); And i used SIMPLE scheme for pressure velocity coupling. I checked dimensions. When i use the wmake command it is compiling with no errors (i think there is no error in the solver). But when i run my case using solver the error message coming like as follows, error message:- Create time Create mesh for time = 0 Reading transportProperties Reading field p Reading field U Reading field T Reading/calculating face flux field phi No field sources present SIMPLE: convergence criteria field p tolerance 0.01 field U tolerance 0.001 field T tolerance 1e-05 Starting time loop Time = 1 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 4.00548e-06, No Iterations 4 DILUPBiCG: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0 #0 Foam::errorprintStack(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 FoamDILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #4 FoamDILUPreconditionerDILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #5 Foam::lduMatrixpreconditioner::addasymMatrixCons tructorToTable<FoamDILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #6 Foam::lduMatrixpreconditioner::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_reddy-2.1.1/platforms/linuxGccDPOpt/bin/myDiffusiveFoam" #10 in "/home/malli_reddy/OpenFOAM/malli_reddy-2.1.1/platforms/linuxGccDPOpt/bin/myDiffusiveFoam" #11 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6" #12 in "/home/malli_reddy/OpenFOAM/malli_reddy-2.1.1/platforms/linuxGccDPOpt/bin/myDiffusiveFoam" Floating point exception malli_reddy@ubuntu:~/OpenFOAM/malli_reddy-2.1.1/PROJECT/myDiffusive$ can anyone suggest me what is the error and how to rectify it. thank you. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Segmentation fault running waveDyMFoam solver (mod. interDyMFoam solver - waves2Foam) | Ed R | OpenFOAM Running, Solving & CFD | 5 | July 2, 2013 12:36 |
problem in running a modified solver | adambarfi | OpenFOAM | 5 | August 10, 2012 16:52 |
PROBLEM IN SOLVER RUNNING PARALLEL | Íñigo | CFX | 3 | February 4, 2007 06:08 |
Problems in a Background running of Intel Fortran | Minsuk Choi | Main CFD Forum | 2 | February 11, 2005 05:55 |
CFX 5.5 | Roued | CFX | 1 | October 2, 2001 17:49 |