|
[Sponsors] |
March 31, 2009, 09:54 |
Help with new viscosity modell
|
#1 |
New Member
Nabil
Join Date: Mar 2009
Posts: 5
Rep Power: 17 |
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... |
|
March 31, 2009, 11:27 |
|
#2 |
Member
Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 89
Rep Power: 17 |
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, ...) |
|
April 1, 2009, 09:21 |
|
#3 |
New Member
Nabil
Join Date: Mar 2009
Posts: 5
Rep Power: 17 |
thanks for the help, i tried the way u'v suggested, but it's still not correct...
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Turbulence viscosity limited | Madhukar Rapaka | FLUENT | 0 | June 26, 2006 04: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 09:00 |
Problem of Turbulent Viscosity Ratio Limited | David Yang | FLUENT | 3 | June 3, 2002 07:13 |