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

New geometry in tutorial mixer2d unphysical solution for fine mesh

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 14, 2009, 13:00
Default Starting from the tutorial mix
  #1
New Member
 
Christina Smuda
Join Date: Mar 2009
Location: Germany
Posts: 12
Rep Power: 8
christinasmuda is on a distinguished road
Starting from the tutorial mixer2d I tried to do simulations of a different geometry with moving mesh and sliding interfaces. The calculations ran well but unfortunately I get very different flow fields for a coarse and a fine mesh. While the calculation with the coarse mesh results in the expected flow field (first picture below), the simulation with a finer mesh ends up with a wrong solution (second picture). I used the same settings for solver and numerical schemes as in mixer 2d. With adjustable time step Co is < 0.5.

I tried to change some solver settings, numerical schemes and a smaller Co, but I didn't get to a satisfying result with a finer mesh.

Could anyone please give me a hint, which settings I might have to change in order to get the physically right solution?

Thanks a lot,
Christina

coarse mesh:
fine mesh:
christinasmuda is offline   Reply With Quote

Old   January 14, 2009, 13:04
Default Sorry, my first shot with pict
  #2
New Member
 
Christina Smuda
Join Date: Mar 2009
Location: Germany
Posts: 12
Rep Power: 8
christinasmuda is on a distinguished road
Sorry, my first shot with pictures didn't work. Second try:

coarse mesh:
fine mesh:
christinasmuda is offline   Reply With Quote

Old   January 14, 2009, 13:13
Default Again I got a Server Time Out.
  #3
New Member
 
Christina Smuda
Join Date: Mar 2009
Location: Germany
Posts: 12
Rep Power: 8
christinasmuda is on a distinguished road
Again I got a Server Time Out. So here's a link to the pictures: link
christinasmuda is offline   Reply With Quote

Old   January 14, 2009, 13:22
Default I do not see anything terribly
  #4
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,763
Rep Power: 21
hjasak will become famous soon enough
I do not see anything terribly wrong with the solution: it seems the fine mesh is showing transient behaviour (try making a movie) and it would be brave to say this is wrong. Looks to me like you are having moving vortices in some of the cavities...

Could you check the boundary conditions (are you uisng movingWallVelocity?). Also, what happens if you let it run longer?

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   January 15, 2009, 11:55
Default Dear Hrv, thanks for your q
  #5
New Member
 
Christina Smuda
Join Date: Mar 2009
Location: Germany
Posts: 12
Rep Power: 8
christinasmuda is on a distinguished road
Dear Hrv,

thanks for your quick reply. I thought about it and tried some more simulations - but I'm still not really convinced. I agree, there might be transient vortices arising. But shouldn't they appear periodically at every pin of the rotor? The viscosity is very high (10 Pas), so the flow is completely laminar.

I put two films at this Link. The first one showing the velocity distribution and the second one the pressure distribution.

Looking at the pressure distribution, in my opinion the pressure is supposed to be the same in the whole flow field using an incompressible fluid (maybe minor changes close to the pins due to the rotation). But during some time steps the pressure in one half of the flow field is very different to the other half.

I also tried the same simulations as a steady state simulation with simpleSRFFOAM (fine mesh). It converged without any problem to the same flow field, the coarse mesh with icoDyMFoam showed. Afterwards I built icoSRFFoam from simpleSRFFoam in order to do a transient calculation with moving reference frame. Here the same problem arises as with the icoDyMFoam: the transient simulation shows a "strange" flow field with no periodicity (neither in space, nor in time).

Best regards,
Christina
christinasmuda is offline   Reply With Quote

Old   January 15, 2009, 21:48
Default Looks like spikes in the press
  #6
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 386
Rep Power: 15
deepsterblue will become famous soon enough
Looks like spikes in the pressure. Is this animation plotted for every time-step? Or is it at intermittent time intervals? Since you're using icoDyMFoam, I'm assuming that there are topo-changes in the mesh. Do you see the pressure variations immediately after a topology change?
Also, how many PISO correctors/non-ortho correctors are you using?
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   January 16, 2009, 08:11
Default Thank you very much for your a
  #7
New Member
 
Christina Smuda
Join Date: Mar 2009
Location: Germany
Posts: 12
Rep Power: 8
christinasmuda is on a distinguished road
Thank you very much for your answer. I tried increasing the number of correctors and now it's converging perfectly to the expected flow field.

Thanks for your help,
Christina
christinasmuda 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
Create fine mesh that grows to coarse mesh (Urgent CZ FLUENT 1 January 3, 2009 11:36
Has anyone created geometry for HVAC tutorial cfx user CFX 2 July 16, 2008 17:26
steady solution on fine mesh Flo Main CFD Forum 2 May 31, 2008 14:55
Variable inletoutlet for dynamic mesh mixer2D case soeren OpenFOAM Running, Solving & CFD 0 May 11, 2008 17:22
how to converge the solution for fine mesh kathiravan CD-adapco 5 August 11, 2006 01:30


All times are GMT -4. The time now is 05:15.