CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

Different results with Symmetry

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 30, 2014, 05:52
Default Different results with Symmetry
  #1
New Member
 
Florian
Join Date: Jan 2014
Posts: 9
Rep Power: 3
Floschi is on a distinguished road
Hallo Everyone,

I simulated a rectangular duct one time with and one time without symmetry boundary. I just wanted to obtain the velocity profile at the cross-section of the channel. But the results are different. (See the pictures below)

Do you have an idea why there is an such a difference?

Kind regards!

With Symmetry



Without Symmetry Boundary
Floschi is offline   Reply With Quote

Old   January 30, 2014, 07:20
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,556
Rep Power: 83
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Before you start thinking about any physical reason for this, first you must make sure your result is an accurate representation of the simulation you set up. Is the simulation fully converged? Is the mesh fine enough (have you done a mesh sensitivity study)?
ghorrocks is offline   Reply With Quote

Old   January 30, 2014, 07:33
Default
  #3
New Member
 
Florian
Join Date: Jan 2014
Posts: 9
Rep Power: 3
Floschi is on a distinguished road
Thank you for your reply.

All simulations converged with a RMS Residual Target below 10^-5. I've done a mesh sensitivity study. I set up different meshes with unstructured tetrahedral objects with three mesh sizes and compared the velocity profiles of the simulations. After the second refinement the max velocity in the middle of the duct does not change anymore. Is that the correct way?
Floschi is offline   Reply With Quote

Old   January 30, 2014, 07:36
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,556
Rep Power: 83
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
It sounds like you have covered the basic numerical accuracy stuff.

Can you explain what you are modelling? Especially what is happening upstream and downstream, with what boundary conditions you have applied.
ghorrocks is offline   Reply With Quote

Old   January 30, 2014, 08:37
Default
  #5
New Member
 
Florian
Join Date: Jan 2014
Posts: 9
Rep Power: 3
Floschi is on a distinguished road
This simulation is just about a straight rectangular duct to obtain the velocity profile. The profile will be used as initial condition for an inlet of a microchannel with different geometries, so you can save on the inlet zone in the following simulations.

The boundary condition at the inlet is a constant velocity over the cross section. The outlet is an OUTLET with average static pressure of 0. The wall is a non slip wall.

The profile is not changing upstream and downstream.

I don't know, why the use of symmetry ends in different results here.

Floschi is offline   Reply With Quote

Old   January 30, 2014, 17:47
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,556
Rep Power: 83
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Have you checked the secondary flows? The symmetry plane might be stopping a secondary flow.

Plot vectors in the YZ plane on your cross sections and I bet the secondary flows are different.
ghorrocks is offline   Reply With Quote

Old   January 31, 2014, 05:10
Default
  #7
New Member
 
Florian
Join Date: Jan 2014
Posts: 9
Rep Power: 3
Floschi is on a distinguished road
Hmm, you mean the vectors in the pictures below? The flows are different but 10^-3 to 10^-4 less than the primary flow. Does this have a effect on the primary flow?

Maybe I don't understand the concept of the symmetry boundary, but if the flow is the same on each side of the plane, a symmetry boundary can be applied.


Symmetry


Without symmetry
Floschi is offline   Reply With Quote

Old   January 31, 2014, 05:19
Default
  #8
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,556
Rep Power: 83
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
You have not shown the one which looks different. The no symmetry and the symmetry plane on the XZ plane look the same in your first post. It is the one with the symmetry plane in the XY plane which looks different.
ghorrocks is offline   Reply With Quote

Old   January 31, 2014, 05:29
Default
  #9
New Member
 
Florian
Join Date: Jan 2014
Posts: 9
Rep Power: 3
Floschi is on a distinguished road
I actually placed the wrong picture in the first post. The profile is different for both symmetry boundaries.

Symmetry

No symmetry
Floschi is offline   Reply With Quote

Old   January 31, 2014, 06:10
Default
  #10
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,556
Rep Power: 83
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
That changes things - now I am thinking of a different error

Did you use a mass flow rate to specify the flow? When you use a symmetry you halve the inlet area so that would result in double the velocity due to the area being halved.
ghorrocks is offline   Reply With Quote

Old   January 31, 2014, 06:21
Default
  #11
New Member
 
Florian
Join Date: Jan 2014
Posts: 9
Rep Power: 3
Floschi is on a distinguished road
I specified the flow as a constant velocity over the cross section. In this case, I think I don't have to halve it.
Floschi is offline   Reply With Quote

Old   January 31, 2014, 07:02
Default
  #12
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,556
Rep Power: 83
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
OK, fair enough.

Can you post your output file?
ghorrocks is offline   Reply With Quote

Old   January 31, 2014, 07:17
Default
  #13
New Member
 
Florian
Join Date: Jan 2014
Posts: 9
Rep Power: 3
Floschi is on a distinguished road
Sure Symmetry.zip
Floschi is offline   Reply With Quote

Old   February 1, 2014, 06:19
Default
  #14
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,556
Rep Power: 83
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Your convergence on the symmetry plane one is quite a bit looser. Can you converge it tighter?

Also, use larger psuedo-time steps than auto time step. Run maybe 20 iterations with a time step multiplier of 10x or 100x. On long thin domains like this it can take a lot of time to get to a true converged solution when the convection time is so long.
ghorrocks is offline   Reply With Quote

Old   February 4, 2014, 04:47
Default
  #15
New Member
 
Florian
Join Date: Jan 2014
Posts: 9
Rep Power: 3
Floschi is on a distinguished road
I did the simulation again with a RMS residual target below 10^-6 and auto timescale with a timescale factor of 10. The .out files are attended Inlet path_symmetry over hight_10-6.zip

It seems that the solution did not change a lot.

Floschi is offline   Reply With Quote

Old   February 4, 2014, 18:12
Default
  #16
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,556
Rep Power: 83
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
This seems very strange.

Can you get the flow rates over these areas? If the flow is faster in one than the other then I cannot see how they are both conserving mass. Also integrate the velocity over the patch to check that as well.
ghorrocks 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
Oscillating Airfoil Poor Results at High k (reduced frequency) dancfd OpenFOAM Running, Solving & CFD 3 November 4, 2013 09:32
Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam) cfdonline2mohsen OpenFOAM 3 October 21, 2013 09:28
[ANSYS Meshing] Effect of Symmetry BC on results ali.qasim ANSYS Meshing & Geometry 4 May 19, 2012 12:22
Transient Run - Output "Time" in partial results? evcelica CFX 2 May 16, 2012 21:36
Properly using symmetry with both CFX and ANSYS workbench for a FSI analysis, help! Cirion0000 CFX 0 July 6, 2009 14:26


All times are GMT -4. The time now is 00:39.