CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Problems about the selection among "turbFoam", "rhoTurbFoam" and "buoyantFoam"

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 16, 2009, 06:00
Default
  #21
Member
 
Jagmohan Meena
Join Date: May 2009
Posts: 30
Rep Power: 8
jmmeena is on a distinguished road
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
jmmeena is offline   Reply With Quote

Old   June 16, 2009, 07:41
Default
  #22
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
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
chiven is offline   Reply With Quote

Old   June 16, 2009, 08:10
Default
  #23
Member
 
Jagmohan Meena
Join Date: May 2009
Posts: 30
Rep Power: 8
jmmeena is on a distinguished road
Can you please tell how and where you implemented "rho = rho0 + psi(p - p0)" ??
jmmeena is offline   Reply With Quote

Old   June 16, 2009, 08:35
Default
  #24
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
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
chiven is offline   Reply With Quote

Old   June 16, 2009, 08:44
Default
  #25
Member
 
Jagmohan Meena
Join Date: May 2009
Posts: 30
Rep Power: 8
jmmeena is on a distinguished road
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
jmmeena is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Run time Selection Mechanism - Some help required to understand jaswi OpenFOAM Programming & Development 2 April 18, 2014 07:44
Needed Benchmark Problems for FSI Mechstud Main CFD Forum 4 July 26, 2011 12:13
Help required to solve Hydraulic related problems aero CFX 0 October 30, 2006 12:00
Some problems with Star CD Micha CD-adapco 0 August 6, 2003 13:55
unstructured grid sreekanth Main CFD Forum 1 August 6, 2001 15:09


All times are GMT -4. The time now is 03:37.