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

Question About Arrhenius viscosity

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 23, 2014, 10:02
Default Question About Arrhenius viscosity
  #1
New Member
 
Janaína de Andrade Silva
Join Date: Mar 2014
Posts: 8
Rep Power: 12
nainasjdr is on a distinguished road
I'm trying implment Arrhenius viscosity.

Foam::tmp<Foam::volScalarField>
Foam::viscosityModels::ArrheniusDin::calcNu() const
{
const volScalarField& T= U_.mesh().lookupObject<volScalarField>("T");
return
nu0_* exp(E_/T);
}

But when I sove my problem this messagem appear



/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.2 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.2-9240f8b967db
Exec : TemppimpleFoam
Date : Jul 23 2014
Time : 09:55:19
Host : "cp-jana"
PID : 6073
Case : /home/jana/OpenFOAM/jana-2.2.2/run/Simulacao_Temperatura/temp/Vaso50micro
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading field T

Reading/calculating face flux field phi

Selecting incompressible transport model ArrheniusDin


--> FOAM FATAL ERROR:
Argument of trancendental function not dimensionless

From function trans(const dimensionSet&)
in file dimensionSet/dimensionSet.C at line 424.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::trans(Foam::dimensionSet const&) at ??:?
#3 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::exp<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#4 Foam::viscosityModels::ArrheniusDin::calcNu() const at ??:?
#5 Foam::viscosityModels::ArrheniusDin::ArrheniusDin( Foam::word const&, Foam::dictionary const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
#6 Foam::viscosityModel::adddictionaryConstructorToTa ble<Foam::viscosityModels::ArrheniusDin>::New(Foam ::word const&, Foam::dictionary const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
#7 Foam::viscosityModel::New(Foam::word const&, Foam::dictionary const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
#8 Foam::singlePhaseTransportModel::singlePhaseTransp ortModel(Foam::GeometricField<Foam::Vector<double> , Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
#9
at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
at ??:?
Someone Can I help me?
nainasjdr is offline   Reply With Quote

Old   July 24, 2014, 00:04
Default
  #2
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18
Cyp is on a distinguished road
Hi!

It seems that your variable E_ is not declared are either a dimensionedScalar or a volScalarField. So you tried to compute the exponential of a variable that has the dimension of the inverse of the temperature. Which is no sense since the exponential is dimensionless...

By the way, I also suggest you to use Foam::exp() instead of exp() in order to specify the compiler that you are going to use the overloaded exp function.

Cheers,
Cyp
Cyp is offline   Reply With Quote

Old   July 27, 2014, 19:24
Default
  #3
New Member
 
Janaína de Andrade Silva
Join Date: Mar 2014
Posts: 8
Rep Power: 12
nainasjdr is on a distinguished road
Quote:
Originally Posted by Cyp View Post
Hi!

It seems that your variable E_ is not declared are either a dimensionedScalar or a volScalarField. So you tried to compute the exponential of a variable that has the dimension of the inverse of the temperature. Which is no sense since the exponential is dimensionless...

By the way, I also suggest you to use Foam::exp() instead of exp() in order to specify the compiler that you are going to use the overloaded exp function.

Cheers,
Cyp
E_ is a dimensionedScalar...
nainasjdr is offline   Reply With Quote

Old   July 28, 2014, 01:21
Default
  #4
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18
Cyp is on a distinguished road
have you tried with Foam::exp() instead of exp() ?
Cyp is offline   Reply With Quote

Old   July 28, 2014, 11:20
Default
  #5
New Member
 
Janaína de Andrade Silva
Join Date: Mar 2014
Posts: 8
Rep Power: 12
nainasjdr is on a distinguished road
Quote:
Originally Posted by Cyp View Post
have you tried with Foam::exp() instead of exp() ?
Yes. I tried it. But not worked



Foam::tmp<Foam::volScalarField>
Foam::viscosityModels::ArrheniusDin::calcNu() const
{
const volScalarField& T= U_.mesh().lookupObject<volScalarField>("T");
return
nu0_* Foam::exp(E_/T);
}

// Variables

dimensionedScalar nu0_;
dimensionedScalar E_;
volScalarField nu_;
nainasjdr is offline   Reply With Quote

Old   July 28, 2014, 17:43
Default
  #6
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18
Cyp is on a distinguished road
and are you sure that E_ and T have the same dimension ?
Cyp is offline   Reply With Quote

Old   July 28, 2014, 19:16
Default
  #7
New Member
 
Janaína de Andrade Silva
Join Date: Mar 2014
Posts: 8
Rep Power: 12
nainasjdr is on a distinguished road
Quote:
Originally Posted by Cyp View Post
and are you sure that E_ and T have the same dimension ?
Thank you for help me.
It's worked
nainasjdr 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
Too low temperature at combustor outlet romekr FLUENT 2 February 6, 2012 11:02
Dynamic Viscosity at transportproperties dictionary NickolasPl Main CFD Forum 0 June 2, 2011 15:06
Question to turbulence viscosity ratio Claud FLUENT 0 October 17, 2006 02:18
Question...Turbulence Intensity & Viscosity ratio Jay FLUENT 1 October 6, 2005 05:41
Problem of Turbulent Viscosity Ratio Limited David Yang FLUENT 3 June 3, 2002 07:13


All times are GMT -4. The time now is 04:44.