[anon_c]

this error can be

a writing error

in the file

src/transportmodel/viscoelastic/make

there is a data calld files

you have to write there the real path line´s name without wiriting mistake

and then you go inside

src/transportmodel/viscoelastic/viscoelastcilaws

and there the name of the (((target))) the folder/directory must be the same like n the data called files (of the make directory)

after that go out side to

src/transportmodels/viscoelastic/

and write in your shell

wmake libso

(another way is go out to src/transportmodells/ and there typ in shell ./Allmake)

but you really need to be shure that in your new models directory under

src/transportmodel/viscoelastic/viscoelastcilaws

the name of folder is the same like (let us say newmodell)
is in the
path

src/transportmodel/viscoelastic/viscoelastcilaws/newmodell

with

newmodell.C
and
newmodell.H

[ata]

Hi Tajoooko
Thank you very much for you quick reply. I'll check that and if I have problem I'll inform you.
Thank you very much.
Best regards

Ata

[purushotam85]

Hello everyone,

I want to simulate die swell and jet impact on moving surface. The fluid under consideration is viscoelastic with shear thinning behaviour. I was reading through these posts and got few viscosity models developed by Jovani.

For my case I want to add new viscosity model that will read one viscoelastic fluid and one Newtonian fluid like air or water together. The question is how do I combine these two very different models together.

And finally I want to solve this using viscoelasticInterFoam. So if anybody has done this before please help.

Any help is appreciated.

Thanks
Purushotam Kumar

[ata]

Hi Purushotam Kumar
I hope you are doing very well. I have been done you desired code in OpenFOAM but because I am PhD student I can not release it before end of my PhD. However you can use interFoam and viscoelasticFluidFoam solvers and combine them to achieve your purpose. I'll be very glad if I can help you.
Best regards

Ata

[purushotam85]

Thanks for your concern ata and I do understand your point.

Just tell me how did you add those multi mode viscosity model with newtonian viscosity model.

I was trying to follow "incompressibleTwoPhaseMixture" viscosity and write it for viscoelastic case but I could not get it going. any help on this.

Regards,
Purushotam Kumar

[ata]

Hi Purushotam Kumar
For your problem see icoFoam and viscoelasticFluidFoam solvers. And in future I'll explain this much more in my paper.
Best regards

Ata

[ata]

Hi Antonio Martins
Thank you very much for your message.
As I know Jovani is working on his PhD thesis and as I know his PhD thesis is on two phase viscoelastic fluid flow. It seems that he can not release his code until end of his PhD thesis.
However is Jovani released his code?
Thank you very much
Best regards

Ata

[desert_1250]

Hello Jovani,how are u?

My name is rasoul and I need PTT model for solving the viscoelastic Fluid,if you help me ,i am very glad,
i am waiting...

[ata]

Hi s.rasoul_varedi
You can find a complete reply to your question in one post that I were asked from Jovani a question like your question.
Best regards

Ata

[jovani]

Hello,

I think there is anything about it in older posts. The line command bellow can download all OF-dev:

svn checkout http://openfoam-extend.svn.sourcefor...nFOAM-1.5-dev/ OpenFOAM-1.5-dev/

and then it need be compiled. For only LPTT the follow line command can be used:

svn checkout http://openfoam-extend.svn.sourcefor...sticLaws/LPTT/ LPTT/


In ubuntu 10.04, there is an easy way to get OF-dev, developed by CAE-Team, see:

http://openfoamwiki.net/index.php/In...ucid_10.04_LTS


Ata,

I am not working with viscoelastic fluids in my PhD thesis, I am working with population balances. The solver based on interFoam, to take account viscoelastic effects, was developded soon I finnished my Mc. There is about one year I dont take it again. I think the solver is not yet closed once I can see for make it better. If you can contribute in this subject in your PhD thesys this will be very good . A small paper with some results for viscoelasticInterFoam you can see here:

http://ayumi106.co.cc/science/articl...ca8189dde0d11f

or,

http://www.aidic.it/escape20/webpapers/87Favero.pdf

Enjoy,

Jovani

[ata]

Hi jovani
Thank you very much and excuse me for my fault about your PhD thesis subject. I have been prepare a code for two phase viscoelastic fluid flows based on your excellent code and interFoam solver for my PhD thesis .Abviously because of your abilities and your relation with Prof. Jasak your code is a better one.
Thank you very much
Best regards

Ata

[ata]

Dear Jovani
Hello
I hope you are doing very well.
I have a question. Did you release your viscoelasticInterFoam code?
Thank you very much
Best regards

Ata

[jovani]

Hello Ata,

The solver was not yet released. People from my old work group are using it to proceed with some simulations.

Best regards,
Jovani

[ata]

Hi Jovani
Thank you very much and congratulation for your excellent football team (Brazil is the first team in the word).
Best regards

Ata

[mohsenkh599]

Hi every body,
I have a problem when I want to compile a model.
the error is:

************************************************** *********************
************************************************** *********************

viscoelasticLaws/mohsen/mohsen.C:159: error: no match for ‘operator==’ in ‘Foam:perator+(const Foam::tmp<Foam::fvMatrix<Type> >&, const Foam::tmp<Foam::fvMatrix<Type> >&) [with Type = Foam::SymmTensor<double>](((const Foam::tmp<Foam::fvMatrix<Foam::SymmTensor<double> > >&)((const Foam::tmp<Foam::fvMatrix<Foam::SymmTensor<double> > >*)(& Foam::fvm::div(const Foam::surfaceScalarField&, Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&) [with Type = Foam::SymmTensor<double>](((Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>&)(&((Foam::mohsen*)this)->Foam::mohsen::a_))))))) == Foam:perator+(const Foam::tmp<Foam::GeometricField<TypeR, PatchField, GeoMesh> >&, const Foam::tmp<Foam::GeometricField<Type1, PatchField, GeoMesh> >&) [with Type1 = Foam::Tensor<double>, Type2 = Foam::SymmTensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >&)((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >*)(& Foam:perator*(const Foam::dimensioned<double>&, const Foam::tmp<Foam::GeometricField<Type, PatchField, GeoMesh> >&) [with Type = Foam::SymmTensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >&)((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >*)(& Foam:perator&(const Foam::tmp<Foam::GeometricField<TypeR, PatchField, GeoMesh> >&, const Foam::tmp<Foam::GeometricField<Type1, PatchField, GeoMesh> >&) [with Type1 = Foam::SymmTensor<double>, Type2 = Foam::SymmTensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >&)((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >*)(& Foam:perator-(const Foam::dimensioned<Type>&, const Foam::tmp<Foam::GeometricField<Type1, PatchField, GeoMesh> >&) [with Form = Foam::SymmTensor<double>, Type = Foam::SymmTensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >&)((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >*)(& Foam:perator*(const Foam::scalar&, const Foam::GeometricField<Type, PatchField, GeoMesh>&) [with Type = Foam::SymmTensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>&)((const Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>*)(&((Foam::mohsen*)this)->Foam::mohsen::a_))))))))))))))))))))’

************************************************** ************
************************************************** ************

and the model is :

// Velocity gradient tensor
volTensorField L = fvc::grad( U() );

// Convected derivate term
// volTensorField C = tau_ & L;

// vorticity term
volTensorField Vor = L - L.T();

// Twice the rate of deformation tensor
volSymmTensorField twoD = twoSymm( L );


// fc
scalar fc = 0.5; //1 - 27 * det( a_ );

//double dot term

// volSymmTensorField mk = -( (2 / 35) * ( 1 - fc ) * twoD - (2 / 7) * ( a_ & twoD + twoD & a_ + (1 / 2) * tr( twoD & a_ ) * I_ ) - fc * tr( twoD & a_ ) * a_ );

// Stress transport equation
tmp<fvSymmTensorMatrix> aEqn
(
fvm::ddt(a_)
+ fvm::div(phi(), a_)
==
( 1 / 2 ) * ( Vor & a_ - a_ & Vor ) + ( 1 / 2 ) * keisi_ * ( twoD & a_ + a_ & twoD + 2 * ( ( (2 / 35) * ( 1 - fc ) * twoD - (2 / 7) * ( 1 - fc ) * ( a_ & twoD + twoD & a_ + (1 / 2) * tr( twoD & a_ ) * I_ ) - fc * tr( twoD & a_ ) * a_ ) ) + 2 * ci_ * twoD & ( I_ - 3 * a_ )
//+ twoSymm( C )
//- zeta_ / 2 * ( (tau_ & twoD) + (twoD & tau_) )
//- fvm::Sp( epsilon_ / etaP_ * tr(tau_) + 1/lambda_, tau_ )
);

aEqn().relax();
solve(aEqn);
// Viscoelastic stress
tau_ = (etaP_ / keisi_) * (a_ - I_);
}

************************************************** ****************
************************************************** ****************

where "a" , "tau" and "I" are symmetric tensors and "Vor" is asymmetric tensor. and the others are constant parameters.

best

[makaveli_lcf]

You don't have an operator "+" defined for the operation, which you perform at the line number 159 of your source code. Unfortunately, it is not possible to verify from your post, which line has number 159. Check if you operate there with the variables of the correct dimensions.
Good luck!

[T.D.]

Hi prof. jasak,

My name is Talib, and i am trying to implement a new solver for the Rheology of Suspensions of non-Brownian particles that i'll hope put online as soon as i'm finished.

I'll need some help please?

I have already implemented a Model for the effective viscosity of suspensions which is known as Krieger-Dougherty Model.

I started from using simpleFoam solver to implement the new solver.

As we know in simpleFoam there exists the nuEff() that to be calculated from the defined viscosityModel by user (in my case) the Krieger-Daugherty Model.

My problem is the following:

I need to use in the same solver a new nuEff2() that will accounts for the stress due to the particles and that must be defined by another viscosityModel Law? I am facing a problem in how to imlplement this nuEff2() knowing that it should be called from the turbulence->laminar...etc directory in simpleFoam?

Simply and more clearly, how can i make to use two different effective viscosities "nu1" and "nu2" in the same solver where "nu1" uses a "viscosityModel1" and "nu2" uses a "viscosityModel2" ? ("viscosityModel1" and "viscosityModel2" will be defined by Me)


Thanks a lot

[ata]

Hi T.D.
How are you bro.?
I think you must define a new abstract class for viscosity model and then derive your model from that.
Best regards

[T.D.]

Hi ata
yes i know that but the problem is i am so poor with c++ programming
can you give a small example pleaaaaase

For example let's say, a simple class for nuEff2() to be equal to like that in crosspowerLaw Model, and how to call and use it in the solver pleaaaase.



My Best Regards,

yours,

T.D.

[ata]

Hi T.D
Use viscosityModel folder (as template for your abstract class) and one of the existing models in src/transportModels/incompressible/viscosityModels as your derived class for start.
Is it OK or you need more help?
Best regards