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

please help me to solve this problem

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 17, 2011, 21:21
Default please help me to solve this problem
  #1
New Member
 
qi bin
Join Date: Jan 2011
Posts: 5
Rep Power: 15
qibintj is on a distinguished road
hello, Foamers

I am trying to modify the scalarTransportFoam solver TO solve these following equation.



# include "readSIMPLEControls.H"
volScalarField nu1=DO2*pow(DT*(1-S),Tor);
volScalarField nu2=Kp*S*S*S/Nu*(-1*Sigma)*0.5*pow((DT/Kp),0.5)*(-4.24*S+3.789*S*S);
volScalarField jc=(1-S)*Tcu*(T/CO2)*exp(-alpha*F*Eta/(R*T1));
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
{
solve
(
//fvm::ddt(T)
fvm::div(nu1,T)
+(1/(4*F))*jc
);
solve
(
//fvm::ddt(S)
fvm::laplacian(nu2,S)
-(1/(2*F))*jc*M
);
}

}


there is always problem.

[COLOR="rgb(139, 0, 0)"]
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::diagonalPreconditioner::diagonalPrecondition er(Foam::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::lduMatrix:reconditioner::addsymMatrixConst ructorToTable<Foam::diagonalPreconditioner>::New(F oam::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam171/lib/linuxGccDPOpt/libfiniteVolume.so"
#8
in "/opt/openfoam171/applications/bin/linuxGccDPOpt/myFoam"
#9
in "/opt/openfoam171/applications/bin/linuxGccDPOpt/myFoam"
#10 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#11
in "/opt/openfoam171/applications/bin/linuxGccDPOpt/myFoam"
Floating point exception
[/COLOR]

please help me to solve this problem. very appreciated.
qibintj 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
Could CFX solve hydraulic jump problem? Andy Chen CFX 0 August 18, 2009 11:13
Can I use OpenFOAM to solve unsteady diffusion problem yongshenglian OpenFOAM Running, Solving & CFD 1 September 17, 2008 13:03
Someone use Icepak to solve microchannel problem ? Bernie FLUENT 1 October 31, 2006 15:29
Can CFX solve problem with structured mesh? Peted CFX 2 June 1, 2004 00:45
how to solve this problem in CFX5 calculation? cfxbeginer CFX 2 May 1, 2003 09:55


All times are GMT -4. The time now is 07:33.