CFD Online Discussion Forums
Page 21 of 24 « First 111920 21 2223 Last »
Show 40 post(s) from this thread on one page

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Viscoelastic Fluid Flows using OpenFOAM The solver viscoelasticFluidFoam (https://www.cfd-online.com/Forums/openfoam-solving/57846-viscoelastic-fluid-flows-using-openfoam-solver-viscoelasticfluidfoam.html)

alqayssi69 June 13, 2016 19:55
Giesekus model
 
HI member
Depending in my experience the using of Giesekus model reveal divergence result with large shear rate value and this is one of drawback in using this model. I really search in the source of divergence and try to find suitable approach .

I would like to know detail information about your geometry so I can more helpful in this issue.

you are free to message me on my G Mail .
engineer.alqayssi@gmail.com

Best Regard

Julya_Shmit June 15, 2016 02:10
1 Attachment(s)
Thanks for your reply, Ahmed.

I am very glad to know that anybody tries to help me because I am at a nonplus.
 What detail information about geometry do you need?

I have already specified geometry, model's parameters and boundary conditions in previous reply (file Geometry.jpg).

alimea February 18, 2017 10:54
viscoelastic model + free convection
 
Dear jovani
Thanks for developping OF for viscoelastic fluids.
I'm a new foamer.
I want to solve "free convection flow and heat transfer of viscoelastic fluid by Buossinesq approximation"
Could you plz help me how to do that?

pippo2013 June 8, 2017 18:18
viscoelasticFluidFoam for dynamic mesh
 
Quote:
Originally Posted by Hiroshiman (Post 528160)
Hi every one,
I'm working on multiphasic flows of viscoelastic fluids. I used the interDyMFoam and viscoelasticFluidFoam from of OE 3.1 to create a interViscoDyMFoam, using the work here.
However, the viscosity from the constant/transportProperties is still used by the twoPhaseMixture operation :
Code:
   
    nuModel1_
    (
    viscosityModel::New
        (
            "nu1",
            subDict(phase1Name_),
            U,
            phi
        )
    ),
    nuModel2_
    (
        viscosityModel::New
        (
            "nu2",
            subDict(phase2Name_),
            U,
            phi
        )
    ),
How can I replace these values read in a dictionnary by values computed from the flow field ? Or maybe there is a simpler approach..


edit:
actually I'd only need to read the viscoelascticProperties from the twoPhaseMixture files. However, the definitions are quite different from reading from the transportProperties files as usual. How could I, for example, read etap1 from constant/viscoelascticProperties in twoPhaseMixture.C ?


Best regards,
Florian

Hi everyone!
I am interested in applying viscoelasticFluidFoam to a dynamic mesh.
Frorian, have you succeded in developing interViscoDyMFoam? The code (cited in the quoted post) is no more linked. Could you please make it available?
Any other help in combining viscoelasticFluidFoam and interDyMFoam would be highly appreciated! :)

Thank you!!

tdof December 3, 2017 05:22
Hi, I've got a question concerning the calculated stresses of the Leonov model. Since OF uses kinematic viscosity and kinematic pressure, are the displayed stresses from tau and sigma also divided by rho?

bfzah December 25, 2017 02:35
vortex in 3D channel
 
Hello
How can I form a vortex in a square-shaped channel?
Or how can I see vortices?
In all the papers, a vortex pair is seen at the corners of the channel with the Giesekus model.
How can this be modeled with open foam?
what is the properties flow ?
B.C ?

please help me?

my Email: farhadhoseinbor@gmail.com
Telegram id: @bfzah

adambarfi December 25, 2017 03:10
Quote:
Originally Posted by bfzah (Post 676031)
Hello
How can I form a vortex in a square-shaped channel?
Or how can I see vortices?
In all the papers, a vortex pair is seen at the corners of the channel with the Giesekus model.
How can this be modeled with open foam?
what is the properties flow ?
B.C ?

please help me?

my Email: farhadhoseinbor@gmail.com
Telegram id: @bfzah
Greetings,

First of all, welcome to CFD-online.

Let me explain a little bit when some vortices are formed in fully developed laminar viscoelastic flow. Having corner in the cross section of a channel in fully developed viscoelstic fluid flow results in a pair of elastic vortices. So, in a channel with square cross section we can see eight elastic vortices, very similar to the vortices in the fully developed turbulent flow of a Newtonian fluid in the same geometry. The intensity of these eight vortices in fully developed laminar viscoelastic flow is around 3 percent while in fully developed turbulent flow of a Newtonian fluid it is around 30%.

I cannot remember if the Weissenberg number may affect the intensity of these vortices, but I certainly know that non-newtonian model such as Oldroyd-B, power-law, etc cannot simulate these vortices.

the boundary conditions are simple. An uniform velocity intet, pressure outlet, and non-slip boundary condition for the walls. The flow should reach the fully developed state in the channel. Just use the viscoelastic properties provided by Jovani in the tutorial case.

Have fun ;)
Mostafa

bfzah December 25, 2017 03:36
Quote:
Originally Posted by adambarfi (Post 676041)
Greetings,

First of all, welcome to CFD-online.

Let me explain a little bit when some vortices are formed in fully developed laminar viscoelastic flow. Having corner in the cross section of a channel in fully developed viscoelstic fluid flow results in a pair of elastic vortices. So, in a channel with square cross section we can see eight elastic vortices, very similar to the vortices in the fully developed turbulent flow of a Newtonian fluid in the same geometry. The intensity of these eight vortices in fully developed laminar viscoelastic flow is around 3 percent while in fully developed turbulent flow of a Newtonian fluid it is around 30%.

I cannot remember if the Weissenberg number may affect the intensity of these vortices, but I certainly know that non-newtonian model such as Oldroyd-B, power-law, etc cannot simulate these vortices.

the boundary conditions are simple. An uniform velocity intet, pressure outlet, and non-slip boundary condition for the walls. The flow should reach the fully developed state in the channel. Just use the viscoelastic properties provided by Jovani in the tutorial case.

Have fun ;)
Mostafa
tnx my dear friend .
i'm very happy that answer my question's.
But I can not do it.
All the details that the rest of the researchers did in their previous work, I did. I do not know why I can not see the vortex

bfzah December 25, 2017 03:45
3 Attachment(s)
this is my file.

adambarfi December 25, 2017 04:04
Quote:
Originally Posted by bfzah (Post 676047)
this is my file.
Your boundary conditions seems ok to me as well as the fluid properties. What is the Reynolds number? How long is the length of your channel? Does the flow get fully developed though it?
How do you visualize the secondary flows?

bfzah December 26, 2017 02:45
Re=rho*u*D/etha
rho=1000, u=0.2, D=0.02, etha=etha(s)+etha(p)=4

Assume a long channel, but I'm not sure that the flow has been developed!

adambarfi December 26, 2017 03:29
Quote:
Originally Posted by bfzah (Post 676128)
Re=rho*u*D/etha
rho=1000, u=0.2, D=0.02, etha=etha(s)+etha(p)=4

Assume a long channel, but I'm not sure that the flow has been developed!
So Re=1; creeping flow! I think that's a little tricky.

Did u assume the flow is steady state? If yes, it is not a concern how the outlet velocity is changing during your simulations. Just let it solve the geometry completely and then draw the secondary flows on a surface at the fully developed section of your duct (NOT at the outlet surface).

On the other hand, plotting the secondary flows is tricky, too. You can google it how to draw secondary flows in paraFoam.

Good luck!

bfzah December 26, 2017 04:00
Quote:
Originally Posted by adambarfi (Post 676129)
So Re=1; creeping flow! I think that's a little tricky.

Did u assume the flow is steady state? If yes, it is not a concern how the outlet velocity is changing during your simulations. Just let it solve the geometry completely and then draw the secondary flows on a surface at the fully developed section of your duct (NOT at the outlet surface).

On the other hand, plotting the secondary flows is tricky, too. You can google it how to draw secondary flows in paraFoam.

Good luck!
yes flow is steady state.
thanx

tdof January 5, 2018 08:34
Hi, I've got a question about the PISO algorithm used in viscoelasticFluidFoam. The official statement in the OF documentation is that it does not use underrelaxation, but the example cases all implement relaxation factors. I therefore assume that the PISO algorithm has been modified to allow their usage, but is a description available somewhere describing this exactly? Edit: I've found it in the source code, but is there a diagram or other explanation available?

Also, does the solver use the complete Navier Stokes equations or only Stokes?

hjasak January 7, 2018 10:16
The solver solves the Navier-Stokes equations

I do not think you need an entire book just for the fact that PISO under-relaxaes the momentum equation to improve stability.

In any case, we are moving to coupled solvers - see recent work from Prof Nobrega.

Hrv

tdof January 8, 2018 06:19
Great, thank you for your reply Prof. Jasak.

bfzah March 4, 2018 10:52
How to create a transient fluid flow in the channel?
What conditions are required in a viscoelasticfluidfoam solver to generate a transient fluid flow?

alimea March 4, 2018 11:00
Quote:
Originally Posted by bfzah (Post 683690)
How to create a transient fluid flow in the channel?
What conditions are required in a viscoelasticfluidfoam solver to generate a transient fluid flow?
Hi farhad
It's like what they do in tutorials. what do you mean about "conditions"?

bfzah March 6, 2018 08:33
Quote:
Originally Posted by alimea (Post 683693)
Hi farhad
It's like what they do in tutorials. what do you mean about "conditions"?
boundary condition for walls and symmetry palne in square channel and (tau,u,p) define in PTT-EXP single mode!

alimea March 24, 2018 11:34
relaxation factors in viscoelasticfluidfoam
 
Hi foamers

I'm working with viscoelasticfluidfoam. We know that this solver is written based on PISO algorithm and this algorithm doesn't have relaxation factor. So why did they use relaxation factors in many examples of tutorials?

What should I do?
I'v solved some problems with relaxation factors (e.g. 0.3 , 0.5) and now some of my friends tell me that they are wrong!

please help me if you are confident about this issue.


All times are GMT -4. The time now is 08:28.
Page 21 of 24 « First 111920 21 2223 Last »
Show 40 post(s) from this thread on one page