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

Viscoelastic Fluid Flows using OpenFOAM The solver viscoelasticFluidFoam

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

Like Tree21Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   November 9, 2012, 06:48
Default
  #221
Member
 
ali jafari
Join Date: Sep 2012
Posts: 49
Rep Power: 4
ali jafari is on a distinguished road
Hi
I have two problem in viscoelastic:

1)I confused about stress boundary conditions in 3D models (walls) , zerogradient or fixedvalue ?

2)I want to solve taylor problem ( rotational concentric cylinders ) but it did not converge , I changed fvschemes fvsolution , but it just cause to changed ill conditions
ali jafari is offline   Reply With Quote

Old   November 9, 2012, 07:24
Default
  #222
Member
 
Sita Drost
Join Date: Mar 2009
Location: Arnhem, The Netherlands
Posts: 46
Rep Power: 7
sita is on a distinguished road
Dear Ali,

1) You can use zerogradient, Jovani Favero explained why somewhere on this forum
2) Maybe your gradients are too sharp? Does mesh refinement help, or does this make it worse? You could try a gradual start-up of the flow, don't exactly remember how this is called in OpenFOAM, something like timeVaryingFixedValue (the OpenFOAM User Manual in combination with this forum and Google will probably give you sufficient information)

Good luck!
sita is offline   Reply With Quote

Old   November 14, 2012, 12:40
Default having a problem with fixedFluxPressure in the viscoelasticFluidFoam package
  #223
New Member
 
Joana Cruz
Join Date: Nov 2012
Posts: 2
Rep Power: 0
Joana is on a distinguished road
Hi!
I'm a recent user of OpenFoam and I need your help about the viscoelasticFluidFoam package. I'm facing an error when trying to use the fixedFluxPressure boundary condition to define the channel inlet/outlet pressure.

I defined the "OpenFOAM/channel/0/p" as follow:

«««««««««««««««««««««««««««««««««««««««««««««««««« «««««««««««««««««««««««««««««««««««««««««««««««««« ««««««««««
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open Source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{

inlet
{
type fixedFluxPressure;
adjoint no;
U U;
phi phi;
rho rho;
}

wall
{
type zeroGradient;
}

outlet
{
type fixedFluxPressure;
adjoint no;
U U;
phi phi;
rho rho;

}

symmetryPlane
{
type symmetryPlane;
}

}
«««««««««««««««««««««««««««««««««««««««««««««««««« «««««««««««««««««««««««««««««««««««««««««««««««««« «««««««««««««

when trying to run my case, the program gives the following error:


«««««««««««««««««««««««««««««««««««««««««««««««««« «««««««««««««««««««««««««««««««««««««««««««««««««« ««««««««««««
--> FOAM FATAL ERROR:

request for volScalarField (1|A(U)) from objectRegistry region0 failed
available objects of type volScalarField are

1
(
p
)

From function objectRegistry::lookupObject<Type>(const word&) const
in file /home/joanacruz/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 131.

FOAM aborting
«««««««««««««««««««««««««««««««««««««««««««««««««« «««««««««««««««««««««««««««««««««««««««««««««««««« ««««

I don't understand what's happen, because I'm using p as the object type in volscalarfield and I already use this boundary condition with icoFoam solver and it works fine.

Thanks in advance
Joana
Joana is offline   Reply With Quote

Old   November 15, 2012, 04:46
Default
  #224
New Member
 
Emil Baric
Join Date: May 2012
Location: Graz. Austria
Posts: 12
Rep Power: 4
emil is on a distinguished road
Hello everybody can someone help me about how can I plot shear components in paraView since the stressComponents postprocessing is not working in this case.


Thank you
emil is offline   Reply With Quote

Old   November 20, 2012, 11:03
Default
  #225
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 283
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Hi everybody
I have some problems with my new solver that is a modified viscoelasticFluidFoam which I added to it temperature and Boussinesq assumption. when I ran it the below error appeared:

Code:
Courant Number mean: 2925.82 max: 9889.16
deltaT = 3.51873e-106
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 982
    Increased the timePrecision from 20 to 21 to distinguish between timeNames at time 2
Time = 2.00000003865963549288

DILUPBiCG:  Solving for Ux, Initial residual = 0.0020599, Final residual = 6.32235e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00485814, Final residual = 3.28543e-09, No Iterations 2
DILUPBiCG:  Solving for T, Initial residual = 0.00258186, Final residual = 1.20616e-08, No Iterations 4
#0  Foam::error::printStack(Foam::Ostream&) in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Uninterpreted: 
#3  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#5  Foam::fvMatrix<double>::solve() in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/bin/nbbvff2"
#6  
 in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/bin/nbbvff2"
#7  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#8  
 in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/bin/nbbvff2"
Floating point exception
I reduced the deltaT even to 1e-15 but this error didn't solved.
is here anybody knows what should I do?

thank you
adambarfi is offline   Reply With Quote

Old   November 20, 2012, 20:52
Thumbs up
  #226
New Member
 
k zhang
Join Date: Jun 2012
Location: China
Posts: 9
Rep Power: 4
skykzhang is on a distinguished road
hello, every formers
i am a newcomer of the viscoelasticFluidFoam. In the tutorial of Leonov constitutive, I am confused with the four parameters(rho,etaS, etaP, lambda). I have seeked some papers to want to make them certain, but ....
About the Leonov constitute equation, there are about 11 parameters to describe the equation in a software.
So I want to know about the parameters in viscoelasticFluidFoam, how can i find the four-parameter model of Leonov (rho, etaS, etaP, lambda)? In some papers?can someone kindly help me?
skykzhang is offline   Reply With Quote

Old   November 21, 2012, 09:24
Default
  #227
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 8
ata is on a distinguished road
Hi
To adambarfi:
It seems you have convergence problem.
To skykzhang
See the Bird's Rheology of Polymers book.
ata is offline   Reply With Quote

Old   November 21, 2012, 11:37
Default
  #228
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 283
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Dear Ata,
thank you for your reply, I fixed it.
But I have some problem with my T contour. it seems to be wrong. I hope that I can fix it, otherwise I will be here again
adambarfi is offline   Reply With Quote

Old   November 25, 2012, 23:50
Default
  #229
New Member
 
k zhang
Join Date: Jun 2012
Location: China
Posts: 9
Rep Power: 4
skykzhang is on a distinguished road
Thank u! Ata.
I will use the solver to simulate polymer injection.
skykzhang is offline   Reply With Quote

Old   December 8, 2012, 09:49
Default
  #230
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 283
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
hi everybody
again I have the same problem I posted 2 weeks ago!!! this error:

Code:
Courant Number mean: 244.174 max: 825.398
deltaT = 3.09674e-105
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 982
    Increased the timePrecision from 21 to 23 to distinguish between timeNames at time 2.64
Time = 2.6400000319743162080499

DILUPBiCG:  Solving for Ux, Initial residual = 0.00206474, Final residual = 6.42745e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00481591, Final residual = 3.30023e-09, No Iterations 2
DILUPBiCG:  Solving for T, Initial residual = 0.00259811, Final residual = 1.11192e-08, No Iterations 4
#0  Foam::error::printStack(Foam::Ostream&) in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Uninterpreted: 
#3  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#5  Foam::fvMatrix<double>::solve() in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/bin/nbbvff2"
#6  
 in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/bin/nbbvff2"
#7  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#8  
 in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/bin/nbbvff2"
Floating point exception
I'm using Giesekus model and my lambda=1e08, rho=4000 and etha0=1e24.
I added temperature equation and Boussinesq assumption to viscoelasticFluidFoam. it works for normal viscoelastic fluids. but when I want to solve my problem with above data, it failed.
anybody knows where is wrong?

thank you
adambarfi is offline   Reply With Quote

Old   January 8, 2013, 10:24
Default Simulation of thixotropic fluids
  #231
New Member
 
Stefan Hildenbrand
Join Date: Aug 2009
Location: Stuttgart, Germany
Posts: 20
Rep Power: 7
shildenbrand is on a distinguished road
Dear FOAMers,

I would like to simulate a process like this in the future:
- Spray injection of a fluid (I guess we can exclude the injection itself, the beginning point would be a given particle distribution)
- Spray hits vertical wall, a wall film develops and the wall film starts moving downwards
- after a while the wall film stops moving due to increasing viscosity

The fluid to be used for the simulation is a thixotropic fluid which is highly viscous under high shear but then gets less viscous after the period of high shear ends.

The original models implemented in OpenFOAM (power law, Cross, Carreau, etc.) are in my opinion not suitable as they do not take care of the temporal change of viscosity. That's why I came across this thread which seems to handle this problem.

So my questions are:
1) Is possible to make such simulations with OF (Well, I guess the answer to this question is: "everything is possible with OF" ) -- with a reasonable effort?
2) Are the models in viscoelasticFluidFoam generally suitable?
3) What is the state of the viscoelasticInterFoam solver already mentioned in this thread?

Thank you for your help - any help is highly appreciated.

Stefan
shildenbrand is offline   Reply With Quote

Old   February 22, 2013, 04:34
Exclamation dear adambarfi,
  #232
Senior Member
 
Tushar C.
Join Date: Jul 2009
Location: India
Posts: 249
Blog Entries: 1
Rep Power: 8
Tushar@cfd is on a distinguished road
Try running coarse mesh for your problem case.

Also, please elaborate some more iterations which you ran in order to better understand the problem with your solver..
Tushar@cfd is offline   Reply With Quote

Old   February 22, 2013, 15:21
Default
  #233
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 283
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Quote:
Originally Posted by Tushar@cfd View Post
Try running coarse mesh for your problem case.

Also, please elaborate some more iterations which you ran in order to better understand the problem with your solver..
thank you Tushar,
I've done it before and my problems were solved.
adambarfi is offline   Reply With Quote

Old   March 4, 2013, 12:11
Default
  #234
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 283
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Hello everybody,

Is here anybody knows how can I define a temperature dependent viscosity in viscoelasticFluidFoam. I before add temperature equation to the solver.
mu(T)=mu1*exp(-gamma*(T-TRef))

thank you
Mostafa

Last edited by adambarfi; March 4, 2013 at 13:46.
adambarfi is offline   Reply With Quote

Old   March 4, 2013, 13:46
Default
  #235
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 283
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
hi everybody again,

I started it with adding a volScalarField& T to the solver and make some modifications. also, since I need gamma and T1, I added them in viscolesticLaws/GiesekusTDependent/GiesekusTDependent.C and GiesekusTDependent.H

when I want to make it the following error appeared:

Code:
Making dependency list for source file viscoelasticLaws/GiesekusTDependent/GiesekusTDependent.C
SOURCE=viscoelasticLaws/GiesekusTDependent/GiesekusTDependent.C ;  g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3  -DNoRepository -ftemplate-depth-100 -I/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/finiteVolume/lnInclude -IlnInclude -I. -I/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude -I/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linuxGccDPOpt/GiesekusTDependent.o
viscoelasticLaws/GiesekusTDependent/GiesekusTDependent.C: In member function ‘virtual Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::GiesekusTDependent::divTau(Foam::volVectorField&) const’:
viscoelasticLaws/GiesekusTDependent/GiesekusTDependent.C:78:50: error: no match for ‘operator-’ in ‘((const Foam::GiesekusTDependent*)this)->Foam::viscoelasticLawTDependent::T - ((const Foam::GiesekusTDependent*)this)->Foam::GiesekusTDependent::Tn_’
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionSet.H:255:29: note: candidates are: Foam::dimensionSet Foam::operator-(const Foam::dimensionSet&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionSet.H:263:29: note:                 Foam::dimensionSet Foam::operator-(const Foam::dimensionSet&, const Foam::dimensionSet&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionedScalar.H:52:19: note:                 Foam::dimensionedScalar Foam::operator-(const Foam::dimensionedScalar&, Foam::scalar)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionedScalar.H:53:19: note:                 Foam::dimensionedScalar Foam::operator-(Foam::scalar, const Foam::dimensionedScalar&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::scalar&, const Foam::UList<double>&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::scalar&, const Foam::tmp<Foam::Field<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::UList<double>&, const Foam::scalar&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::tmp<Foam::Field<double> >&, const Foam::scalar&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::label&, const Foam::UList<int>&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::label&, const Foam::tmp<Foam::Field<int> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::UList<int>&, const Foam::label&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::tmp<Foam::Field<int> >&, const Foam::label&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::UList<Foam::DiagTensor<double> >&, const Foam::UList<Foam::Tensor<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::UList<Foam::DiagTensor<double> >&, const Foam::tmp<Foam::Field<Foam::Tensor<double> > >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::tmp<Foam::Field<Foam::DiagTensor<double> > >&, const Foam::UList<Foam::Tensor<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::tmp<Foam::Field<Foam::DiagTensor<double> > >&, const Foam::tmp<Foam::Field<Foam::Tensor<double> > >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::diagTensor&, const Foam::UList<Foam::Tensor<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::diagTensor&, const Foam::tmp<Foam::Field<Foam::Tensor<double> > >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::UList<Foam::DiagTensor<double> >&, const Foam::tensor&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::tmp<Foam::Field<Foam::DiagTensor<double> > >&, const Foam::tensor&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/vectorTensorTransformI.H:283:36: note:                 Foam::vectorTensorTransform Foam::operator-(const Foam::vectorTensorTransform&, const Foam::vector&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/quaternionI.H:273:25: note:                 Foam::quaternion Foam::operator-(const Foam::quaternion&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/quaternionI.H:279:25: note:                 Foam::quaternion Foam::operator-(const Foam::quaternion&, const Foam::quaternion&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/septernionI.H:190:25: note:                 Foam::septernion Foam::operator-(const Foam::septernion&, const Foam::vector&)
viscoelasticLaws/GiesekusTDependent/GiesekusTDependent.C:84:53: error: no match for ‘operator-’ in ‘((const Foam::GiesekusTDependent*)this)->Foam::viscoelasticLawTDependent::T - ((const Foam::GiesekusTDependent*)this)->Foam::GiesekusTDependent::Tn_’
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionSet.H:255:29: note: candidates are: Foam::dimensionSet Foam::operator-(const Foam::dimensionSet&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionSet.H:263:29: note:                 Foam::dimensionSet Foam::operator-(const Foam::dimensionSet&, const Foam::dimensionSet&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionedScalar.H:52:19: note:                 Foam::dimensionedScalar Foam::operator-(const Foam::dimensionedScalar&, Foam::scalar)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionedScalar.H:53:19: note:                 Foam::dimensionedScalar Foam::operator-(Foam::scalar, const Foam::dimensionedScalar&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::scalar&, const Foam::UList<double>&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::scalar&, const Foam::tmp<Foam::Field<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::UList<double>&, const Foam::scalar&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::tmp<Foam::Field<double> >&, const Foam::scalar&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::label&, const Foam::UList<int>&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::label&, const Foam::tmp<Foam::Field<int> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::UList<int>&, const Foam::label&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::tmp<Foam::Field<int> >&, const Foam::label&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::UList<Foam::DiagTensor<double> >&, const Foam::UList<Foam::Tensor<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::UList<Foam::DiagTensor<double> >&, const Foam::tmp<Foam::Field<Foam::Tensor<double> > >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::tmp<Foam::Field<Foam::DiagTensor<double> > >&, const Foam::UList<Foam::Tensor<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::tmp<Foam::Field<Foam::DiagTensor<double> > >&, const Foam::tmp<Foam::Field<Foam::Tensor<double> > >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::diagTensor&, const Foam::UList<Foam::Tensor<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::diagTensor&, const Foam::tmp<Foam::Field<Foam::Tensor<double> > >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::UList<Foam::DiagTensor<double> >&, const Foam::tensor&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::tmp<Foam::Field<Foam::DiagTensor<double> > >&, const Foam::tensor&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/vectorTensorTransformI.H:283:36: note:                 Foam::vectorTensorTransform Foam::operator-(const Foam::vectorTensorTransform&, const Foam::vector&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/quaternionI.H:273:25: note:                 Foam::quaternion Foam::operator-(const Foam::quaternion&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/quaternionI.H:279:25: note:                 Foam::quaternion Foam::operator-(const Foam::quaternion&, const Foam::quaternion&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/septernionI.H:190:25: note:                 Foam::septernion Foam::operator-(const Foam::septernion&, const Foam::vector&)
viscoelasticLaws/GiesekusTDependent/GiesekusTDependent.C: In member function ‘virtual void Foam::GiesekusTDependent::correct()’:
viscoelasticLaws/GiesekusTDependent/GiesekusTDependent.C:106:27: error: no match for ‘operator-’ in ‘((Foam::GiesekusTDependent*)this)->Foam::viscoelasticLawTDependent::T - ((Foam::GiesekusTDependent*)this)->Foam::GiesekusTDependent::Tn_’
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionSet.H:255:29: note: candidates are: Foam::dimensionSet Foam::operator-(const Foam::dimensionSet&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionSet.H:263:29: note:                 Foam::dimensionSet Foam::operator-(const Foam::dimensionSet&, const Foam::dimensionSet&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionedScalar.H:52:19: note:                 Foam::dimensionedScalar Foam::operator-(const Foam::dimensionedScalar&, Foam::scalar)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionedScalar.H:53:19: note:                 Foam::dimensionedScalar Foam::operator-(Foam::scalar, const Foam::dimensionedScalar&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::scalar&, const Foam::UList<double>&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::scalar&, const Foam::tmp<Foam::Field<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::UList<double>&, const Foam::scalar&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/scalarField.H:78:1: note:                 Foam::tmp<Foam::Field<double> > Foam::operator-(const Foam::tmp<Foam::Field<double> >&, const Foam::scalar&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::label&, const Foam::UList<int>&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::label&, const Foam::tmp<Foam::Field<int> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::UList<int>&, const Foam::label&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:54:1: note:                 Foam::tmp<Foam::Field<int> > Foam::operator-(const Foam::tmp<Foam::Field<int> >&, const Foam::label&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::UList<Foam::DiagTensor<double> >&, const Foam::UList<Foam::Tensor<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::UList<Foam::DiagTensor<double> >&, const Foam::tmp<Foam::Field<Foam::Tensor<double> > >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::tmp<Foam::Field<Foam::DiagTensor<double> > >&, const Foam::UList<Foam::Tensor<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:61:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::tmp<Foam::Field<Foam::DiagTensor<double> > >&, const Foam::tmp<Foam::Field<Foam::Tensor<double> > >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::diagTensor&, const Foam::UList<Foam::Tensor<double> >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::diagTensor&, const Foam::tmp<Foam::Field<Foam::Tensor<double> > >&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::UList<Foam::DiagTensor<double> >&, const Foam::tensor&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/diagTensorField.H:64:1: note:                 Foam::tmp<Foam::Field<Foam::Tensor<double> > > Foam::operator-(const Foam::tmp<Foam::Field<Foam::DiagTensor<double> > >&, const Foam::tensor&)
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/vectorTensorTransformI.H:283:36: note:                 Foam::vectorTensorTransform 
.
.
.
the .C and .H files are attached.
any body knows what should i do?

thanks
Mostafa
Attached Files
File Type: c GiesekusTDependent.C (3.4 KB, 20 views)
File Type: h GiesekusTDependent.H (3.8 KB, 11 views)

Last edited by adambarfi; March 4, 2013 at 14:33.
adambarfi is offline   Reply With Quote

Old   March 11, 2013, 07:51
Post Adambarfi
  #236
Senior Member
 
Tushar C.
Join Date: Jul 2009
Location: India
Posts: 249
Blog Entries: 1
Rep Power: 8
Tushar@cfd is on a distinguished road
Try adding the term effectively... store & update the term in a loop...

Last edited by Tushar@cfd; March 12, 2013 at 03:08.
Tushar@cfd is offline   Reply With Quote

Old   March 11, 2013, 11:11
Default
  #237
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 283
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Quote:
Originally Posted by Tushar@cfd View Post
Try adding the term effectively... store & update the term in a loop...
my problem was solved.
The idea was: don't put T into the argument lists. Instead look up the T volScalarField from the "registry"!!!

thanks
adambarfi is offline   Reply With Quote

Old   March 20, 2013, 03:36
Angry
  #238
New Member
 
salame ama
Join Date: Dec 2012
Posts: 28
Rep Power: 3
salame is on a distinguished road
HI, Formers:
i want to combine the viscoelasticFluidFoam and InterFoam.
Now, about the transportProperties file, write the code as it in the InterFoam tutorial:

phase1
{

type multiMode;

models
(
first
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 803.87097;
etaS etaS [1 -1 -1 0 0 0 0] 0.002;
etaP etaP [1 -1 -1 0 0 0 0] 0.04;
lambda lambda [0 0 1 0 0 0 0] 0.6855;
alpha alpha [0 0 0 0 0 0 0] 0.5;
}

second
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 803.87097;
.....
phase2
{

type multiMode;

models
(
first
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 1;
etaS etaS [1 -1 -1 0 0 0 0] 0;
etaP etaP [1 -1 -1 0 0 0 0] 0.0000148;
lambda lambda [0 0 1 0 0 0 0] 0;
alpha alpha [0 0 0 0 0 0 0] 0;
}
.....


but it said:
Reading transportProperties



--> FOAM FATAL IO ERROR:
keyword transportModel is undefined in dictionary "/home/caelinux/zk/try/constant/transportProperties:hase1"

file: /home/caelinux/zk/try/constant/transportProperties:hase1 from line 20 to line 64.

From function dictionary::lookupEntry(const word&, bool, bool) const
in file db/dictionary/dictionary.C at line 396.


can u give me a hand?
salame is offline   Reply With Quote

Old   March 20, 2013, 04:31
Default
  #239
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 283
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Greeting, welcome and happy your Nowruz day!!!!

did you modify the transportModel and your new solver files so that they contain both viscoelasticFluidFoam and interFoam information???
adambarfi is offline   Reply With Quote

Old   March 20, 2013, 07:38
Post Hello salame,
  #240
Senior Member
 
Tushar C.
Join Date: Jul 2009
Location: India
Posts: 249
Blog Entries: 1
Rep Power: 8
Tushar@cfd is on a distinguished road
Add these lines:
In EXE_INC = \
"-I$(LIB_SRC)/transportModels/viscoelastic/lnInclude \

EXE_LIBS = \
-viscoelasticTransportModels \
Tushar@cfd 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
VOF simulation of a viscoelastic fluid sinah OpenFOAM Running, Solving & CFD 10 November 25, 2010 11:02
FREE SURFACE VISCOELASTIC FLOWS Valdemir G. Ferreira Main CFD Forum 6 December 18, 2009 06:14
Viscoelastic flow modeling in OpenFOAM vulda OpenFOAM Running, Solving & CFD 1 March 17, 2008 08:32
Polyflow & OpenFoam on Viscoelastic flow modeling Sumeshen Main CFD Forum 0 March 14, 2008 09:29
Viscoelastic fluid codes joel davison Main CFD Forum 0 November 6, 2001 05:09


All times are GMT -4. The time now is 23:50.