[jmmeena]

Hi Chevin,

Sry for not replying for a long time. In this time I was using hit and trial method to solve my problem and end up with these results:

1 main problem in understanding rho value is that it defines rho value at 2 points !! 1st in specie.H and 2nd in basic classes where I am not able to change the value of psi. So I did not change anything in basic classes. Finally in buoyantSimpleFoam solver I used "rho = rho0 - psi*p0 + thermo->rho" but it is not same as "rho = rho0 + psi(p - p0)" because psi is not same as called by thermo class and defined for water. However, variation from rho0 is very small. What you did to solve this problem. Can you please post your codes of basic thermo classes ? I just want to know your approach for this to make solver better.

As far as divergence of my approximate solver is concerned it is due to bad mesh and I am working on mesh currently. Because I tried everything is fvsolution and fvschemes files but it diverged sooner or later !!

As far as H2O class is concerned, I am working on it. And will definitely post useful information.

JM

[chiven]

Hello, JM, sorry about my slow progress. This time I am still struggling with the problem of,


[3] Maximum number of iterations exceeded#0 _ZN4Foam5error10printStackERNS_7OstreamE-0xb48690
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.5/lib/linuxIA64GccDPOpt/libOpenFOAM.so"
#1 _ZN4Foam5error5abortEv-0x1011710
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.5/lib/linuxIA64GccDPOpt/libOpenFOAM.so"
#2 _ZN4Foam7hThermoINS_11pureMixtureINS_14constTransp ortINS_12specieThermoINS_12hConstThermoINS_11const antRhoEEEEEEEEEE9calculateEv-0x1590590
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.5/lib/linuxIA64GccDPOpt/libbasicThermophysicalModels.so"
#3 _ZN4Foam7hThermoINS_11pureMixtureINS_14constTransp ortINS_12specieThermoINS_12hConstThermoINS_11const antRhoEEEEEEEEEE7correctEv-0x15714c0
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.5/lib/linuxIA64GccDPOpt/libbasicThermophysicalModels.so"
#4 main in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxIA64GccDPOpt/buoyantFoam"
#5 __libc_start_main-0x7492e0
in "/lib/tls/libc.so.6.1"
#6 _start in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxIA64GccDPOpt/buoyantFoam"
[3]
[3]
[3] From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const
[3] in file /home/g2/e090012/OpenFOAM/OpenFOAM-1.5/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 78.
[3]

Chiven

[jmmeena]

Can you please tell how and where you implemented "rho = rho0 + psi(p - p0)" ??

[chiven]

Hi, JM, I revised the code of constantRhoI.H in the directory of /src/thermophysicalModels/specie/equationOfState/constantRho as follows.

// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
//- Return density [kg/m^3]
inline scalar constantRho::rho(scalar p, scalar psi) const
{
return 766.52 + psi * (p - 15.0e+06);
}

Chiven

[jmmeena]

Hi,

What about function of psi or you used some constant value ? After this constantrho.H, it again compute rho in basic thermo files, have you changed something there ?

do u have gmail id where we can discuss(chat) and sort it out fast ? mine is jmmeenaiit@gmail.com !!

JM