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

Help with new viscosity modell

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 31, 2009, 08:54
Default Help with new viscosity modell
  #1
New Member
 
Nabil
Join Date: Mar 2009
Posts: 5
Rep Power: 8
Mikro is on a distinguished road
Hello everyone,

I'm trying to implement a new viscosity modell for blood (well several acctually) and i'm having a problem with this one: (Powell-Eyring modell)

nuInf_+(nu0_-nuInf_)*scalar(1)/(sinh(lambda_*strainRate())*lambda_*strainRate());


The linking of the library goes without any problem, it's when i run nonNewtonianIcoFoam that the solver crashes, I was hoping that someone with better knowledge in OpenFoam might tell me were the problem is...
This is errors i get when i try to run it, incase it helps anyone to se it better...


#0 Foam::error:rintStack(Foam::Ostream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0xb7fd7420]
#3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&)
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensioned<double> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&)
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)
#6 Foam::viscosityModels::PowellEyring::calcNu() const
#7 Foam::viscosityModels::PowellEyring::PowellEyring( Foam::word const&, Foam::dictionary const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::surfaceMesh> const&)
#8 Foam::viscosityModel::adddictionaryConstructorToTa ble<Foam::viscosityModels::PowellEyring>::New(Foam ::word const&, Foam::dictionary const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::surfaceMesh> const&)
#9 Foam::viscosityModel::New(Foam::word const&, Foam::dictionary const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::surfaceMesh> const&)
#10 Foam::singlePhaseTransportModel::singlePhaseTransp ortModel(Foam::GeometricField<Foam::Vector<double> , Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::surfaceMesh> const&)
#11 main
#12 __libc_start_main
#13 __gxx_personality_v0 at ../sysdeps/i386/elf/start.S:122

Also i might poit out that this is the sixth modell i'v done, and the other five work perfectly fine with nonNewtonianIcoFoam...
Mikro is offline   Reply With Quote

Old   March 31, 2009, 10:27
Default
  #2
Member
 
Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 80
Rep Power: 8
juho is on a distinguished road
It's a divide by zero.

#1 Foam::sigFpe::sigFpeHandler(int)
Fpe = floating point exception

#3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&)
"divide"


I'm not 100% about the syntax, but you should use something like:

dimensionedScalar small('small', dimensionSet(<correct dimensions>),1e-10);

.../max(small, ...)
juho is offline   Reply With Quote

Old   April 1, 2009, 08:21
Default
  #3
New Member
 
Nabil
Join Date: Mar 2009
Posts: 5
Rep Power: 8
Mikro is on a distinguished road
thanks for the help, i tried the way u'v suggested, but it's still not correct...
Mikro 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
Turbulence viscosity limited Madhukar Rapaka FLUENT 0 June 26, 2006 03:17
Turbulent viscosity limited to viscosity ratio... Cyril FLUENT 8 February 7, 2006 21:34
kinematic viscosity at diff temperatures,pressures Mecobio Main CFD Forum 0 November 7, 2005 13:55
Turbulent viscosity in a riser ap FLUENT 8 April 19, 2003 08:00
Problem of Turbulent Viscosity Ratio Limited David Yang FLUENT 3 June 3, 2002 06:13


All times are GMT -4. The time now is 08:30.