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

Lift amplitude error in viscoelastic Fluid

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 16, 2018, 05:40
Default Lift amplitude error in viscoelastic Fluid
  #1
Senior Member
 
A. Min
Join Date: Mar 2015
Posts: 279
Rep Power: 8
alimea is on a distinguished road
Dear Foamers

I'm simulating the shedding flow around a circular cylinder in viscoelastic fluid with rheoFoam.
At first I was working with viscoelasticFluidFoam that is just allow me to use upwind scheme for div(phi,tau) that gave me 30% error in CL amplitude. If I use any other scheme insted of that, the solution would be diverged!
I reduced Courant number to 0.01 and ​the grid became smaller. but again with upwind scheme I couldn't get appropriate result and I had 25% error in CL amplitude!

So I installed openFOAM v 4.0 and compile rheoFoam and read the related part of user manual. Then I tried to get results of Fig. 4 of Oliviera [1] with this condition:
FENE-CR, Re=100 , De=80 , L2=100
ddt: CN 1 - div(phi,U) , div(phi,tau) , div(phi,A) : GaussDefCmpw cubista
SIMPLE algorithm - URF=1 - nInIter=1 - Courant Number = 0.01

But the results are very strange: CLmax=0.4 (Oliviera: 0.15) .
Also there are some strange fluctuations that I don't knw what they are!
1- Could you please tell me what is wrong in my simulation?



Note: To calculate CL and CD I have written this function:
stress = tau+etta_S*(grad(U)+gradT(U)) - p
F= - stress & A[i]
Drag=F.x();
Lift=F.y();
CD=Drag/ (0.5*rho*U_inf^2*Aref)
CL=Lift/ (0.5*rho*U_inf^2*Aref)

3- what is difference between FENE-CR and FENE-CRlog?

Validating this problem has been my problem for a long period of time.

I appreciate if any answers me.
--------------------------------------------
[1] Method for time-dependent simulations of viscoelastic flows: vortex shedding behind cylinder, Paulo J. Oliveira, J. Non-Newtonian Fluid Mech. 101 (2001) 113137
alimea is offline   Reply With Quote

Old   April 17, 2018, 05:26
Default
  #2
Senior Member
 
A. Min
Join Date: Mar 2015
Posts: 279
Rep Power: 8
alimea is on a distinguished road
Quote:
Originally Posted by alimea View Post
Dear Foamers

I'm simulating the shedding flow around a circular cylinder in viscoelastic fluid with rheoFoam.
At first I was working with viscoelasticFluidFoam that is just allow me to use upwind scheme for div(phi,tau) that gave me 30% error in CL amplitude. If I use any other scheme insted of that, the solution would be diverged!
I reduced Courant number to 0.01 and ​the grid became smaller. but again with upwind scheme I couldn't get appropriate result and I had 25% error in CL amplitude!

So I installed openFOAM v 4.0 and compile rheoFoam and read the related part of user manual. Then I tried to get results of Fig. 4 of Oliviera [1] with this condition:
FENE-CR, Re=100 , De=80 , L2=100
ddt: CN 1 - div(phi,U) , div(phi,tau) , div(phi,A) : GaussDefCmpw cubista
SIMPLE algorithm - URF=1 - nInIter=1 - Courant Number = 0.01

But the results are very strange: CLmax=0.4 (Oliviera: 0.15) .
Also there are some strange fluctuations that I don't knw what they are!
1- Could you please tell me what is wrong in my simulation?



Note: To calculate CL and CD I have written this function:
stress = tau+etta_S*(grad(U)+gradT(U)) - p
F= - stress & A[i]
Drag=F.x();
Lift=F.y();
CD=Drag/ (0.5*rho*U_inf^2*Aref)
CL=Lift/ (0.5*rho*U_inf^2*Aref)

3- what is difference between FENE-CR and FENE-CRlog?

Validating this problem has been my problem for a long period of time.

I appreciate if any answers me.
--------------------------------------------
[1] Method for time-dependent simulations of viscoelastic flows: vortex shedding behind cylinder, Paulo J. Oliveira, J. Non-Newtonian Fluid Mech. 101 (2001) 113137
any answer?
alimea is offline   Reply With Quote

Reply

Tags
cylinder, lift, rheofoam, viscoelastic, viscoelasticfluidfoam

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
Intl Conf Computational Methods in Fluid Power Jacek Stecki Main CFD Forum 0 November 10, 2002 06:49
Fluid mechanics / drag & lift Jamcam Main CFD Forum 0 October 31, 2002 05:39
My Revised "Time Vs Energy" Article For Review Abhi Main CFD Forum 2 July 9, 2002 10:08
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11
Viscoelastic fluid codes joel davison Main CFD Forum 0 November 6, 2001 06:09


All times are GMT -4. The time now is 09:57.