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

Turbomachinery Mass imbalance

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 31, 2019, 04:15
Default Turbomachinery Mass imbalance
  #1
Senior Member
 
sheaker's Avatar
 
Oskar
Join Date: Nov 2015
Location: Poland
Posts: 184
Rep Power: 11
sheaker is on a distinguished road
Dear All,

I am doing a simulation of a Axial Turbine stage in CFX Turbomachinery. Turbine stage have 15 stationary blades and 23 rotational blades. My geometry and mesh is 1 stationary blade and 2 rotational blades. Inlet and outlet BC's are Mass Flow at Inlet and Static Pressure at Outlet. Mesh is about 8mln elements. I am performing 200iterations while the residuals stabilizing after about 130 iterations.

After performing a lot of simulations with small rotational blade geometry modifications I found out that ~1.7% of mass was removed from inlet and ~3.5% mass was added at GGI interface. Mass flow is in kg/s.


Example:
MassFlow at S1Inlet = 0.023844
MassFlow at S1toR1 = -0.0238434
MassFlow at R1toS1 = 0.0322807
MassFlow at R1Outlet = -0.0322806


Then to calculate mass flow at interface through entire machine:
S1toR1: -0.0238434 * 15 = -0.357651

R1toS1: 0.0322807 * 23 / 2 = 0.37122805


Mass flow at inlet should be 0.3619



Mass flow at inlet is:
0.023844 * 15 = 0.35766


Imbalance at Inlet:
0.3619 - 0.35766 = 0.00424
0.00424/0.3619*100% = 1.1716% of mass loss at S1 Inlet



Imbalance at GGI interface:
0.37122805 - 0.357651 = 0.01357705
0.01357705/0.37122805*100% = 3.657% of mass was added.


Could You, please explain why the mass flow through inlet is not what I set and why the mass is added on GGI interface?
How could I possibly prevent it?


Best regards,
Oskar

Last edited by sheaker; August 31, 2019 at 05:47.
sheaker is offline   Reply With Quote

Old   September 1, 2019, 01:59
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,854
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
It sounds like your simulation is not converged. In an unconverged simulation the results can show anything.

Your comment ".... the residuals stabilizing after about 130 iterations" suggests this FAQ describes what you need to consider: https://www.cfd-online.com/Wiki/Ansy...gence_criteria
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
ghorrocks is offline   Reply With Quote

Old   September 1, 2019, 04:11
Default
  #3
Senior Member
 
sheaker's Avatar
 
Oskar
Join Date: Nov 2015
Location: Poland
Posts: 184
Rep Power: 11
sheaker is on a distinguished road
Dear Glenn,


Thank You for Your reply.
I checked maximum Mass residuals by volume rendering and by point for both domains and I think they are pretty good. Please, see the attachments.
I think it is not related to residuals because mass is well preserved within domain. Changes are present at inlet and GGI interface.



Best regards,
Oskar
Attached Images
File Type: png S1Volume.png (86.4 KB, 36 views)
File Type: png R1Volume.png (76.5 KB, 24 views)
File Type: png S1MaxPoint.png (47.1 KB, 21 views)
File Type: png R1MaxPoint.png (48.6 KB, 14 views)
sheaker is offline   Reply With Quote

Old   September 2, 2019, 14:16
Default
  #4
Senior Member
 
Join Date: Jun 2009
Posts: 1,873
Rep Power: 33
Opaque will become famous soon enough
I think you have confused residuals with mass imbalance.

The mass imbalance is reported in the output file separately from residuals in the flow summary.

You could probably post those here.
Opaque is offline   Reply With Quote

Old   September 2, 2019, 15:03
Default
  #5
Senior Member
 
sheaker's Avatar
 
Oskar
Join Date: Nov 2015
Location: Poland
Posts: 184
Rep Power: 11
sheaker is on a distinguished road
Thank You for Your reply.


Quote:
+--------------------------------------------------------------------+
| U-Mom-R1 |
+--------------------------------------------------------------------+
Boundary : R1 Blade -1.1978E+01
Boundary : R1 Default -3.2658E-01
Boundary : R1 Hub -4.1550E+00
Boundary : R1 Outlet -3.4936E-01
Boundary : R1 Shroud 6.9161E+00
Domain Src (Neg) : R1 -5.3682E-01
Domain Src (Pos) : R1 2.0972E-01
Domain Interface : R1 to R1 Periodic 1 (Side 1) -2.2746E+01
Domain Interface : R1 to R1 Periodic 1 (Side 2) 2.1026E+01
Domain Interface : R1 to S1 (Side 1) 1.1941E+01
-----------
Domain Imbalance : 2.0906E-05

+--------------------------------------------------------------------+
| V-Mom-R1 |
+--------------------------------------------------------------------+
Boundary : R1 Blade -8.9358E+00
Boundary : R1 Default -6.2172E-03
Boundary : R1 Hub 8.3788E-01
Boundary : R1 Outlet -3.4083E+01
Boundary : R1 Shroud -6.6731E-02
Domain Interface : R1 to R1 Periodic 1 (Side 1) -9.7446E+00
Domain Interface : R1 to R1 Periodic 1 (Side 2) 9.7414E+00
Domain Interface : R1 to S1 (Side 1) 4.2257E+01
-----------
Domain Imbalance : 1.2632E-04

+--------------------------------------------------------------------+
| W-Mom-R1 |
+--------------------------------------------------------------------+
Boundary : R1 Blade -6.2494E-03
Boundary : R1 Default 2.0886E+00
Boundary : R1 Hub 1.4441E+01
Boundary : R1 Outlet 1.3347E-01
Boundary : R1 Shroud -3.3029E+01
Domain Src (Neg) : R1 -9.6557E-02
Domain Src (Pos) : R1 2.4863E+00
Domain Interface : R1 to R1 Periodic 1 (Side 1) 3.0604E+00
Domain Interface : R1 to R1 Periodic 1 (Side 2) 9.2045E+00
Domain Interface : R1 to S1 (Side 1) 1.7174E+00
-----------
Domain Imbalance : -1.5751E-06

+--------------------------------------------------------------------+
| P-Mass-R1 |
+--------------------------------------------------------------------+
Boundary : R1 Outlet -3.1318E-02
Domain Interface : R1 to R1 Periodic 1 (Side 1) 4.2258E-03
Domain Interface : R1 to R1 Periodic 1 (Side 2) -4.2258E-03
Domain Interface : R1 to S1 (Side 1) 3.1318E-02
-----------
Domain Imbalance : 1.6682E-07

+--------------------------------------------------------------------+
| U-Mom-S1 |
+--------------------------------------------------------------------+
Boundary : S1 Blade 8.6572E+00
Boundary : S1 Hub -2.1917E+01
Boundary : S1 Inlet -1.0853E-06
Boundary : S1 Shroud 3.6294E+01
Domain Interface : R1 to S1 (Side 2) -9.1612E+00
Domain Interface : S1 to S1 Periodic 1 (Side 1) 7.3396E+01
Domain Interface : S1 to S1 Periodic 1 (Side 2) -8.7269E+01
-----------
Domain Imbalance : 8.6873E-05

+--------------------------------------------------------------------+
| V-Mom-S1 |
+--------------------------------------------------------------------+
Boundary : S1 Blade -9.4145E+00
Boundary : S1 Hub -1.8097E+00
Boundary : S1 Inlet 4.6025E+01
Boundary : S1 Shroud -2.8937E+00
Domain Interface : R1 to S1 (Side 2) -3.2152E+01
Domain Interface : S1 to S1 Periodic 1 (Side 1) -2.7227E+01
Domain Interface : S1 to S1 Periodic 1 (Side 2) 2.7471E+01
-----------
Domain Imbalance : 1.1898E-04

+--------------------------------------------------------------------+
| W-Mom-S1 |
+--------------------------------------------------------------------+
Boundary : S1 Blade 3.5618E+00
Boundary : S1 Hub 4.6615E+01
Boundary : S1 Inlet 2.2477E-06
Boundary : S1 Shroud -8.3417E+01
Domain Interface : R1 to S1 (Side 2) -8.8898E-01
Domain Interface : S1 to S1 Periodic 1 (Side 1) 4.9656E+01
Domain Interface : S1 to S1 Periodic 1 (Side 2) -1.5527E+01
-----------
Domain Imbalance : 1.2597E-05

+--------------------------------------------------------------------+
| P-Mass-S1 |
+--------------------------------------------------------------------+
Boundary : S1 Inlet 2.3842E-02
Domain Interface : R1 to S1 (Side 2) -2.3840E-02
Domain Interface : S1 to S1 Periodic 1 (Side 1) -3.8574E-04
Domain Interface : S1 to S1 Periodic 1 (Side 2) 3.8574E-04
-----------
Domain Imbalance : 1.4838E-06

+--------------------------------------------------------------------+
| H-Energy-R1 |
+--------------------------------------------------------------------+
Boundary : R1 Outlet 5.4540E+03
Bnd Src/Visc Work: R1 Outlet 3.2580E-03
Bnd Src/Visc Work: R1 Shroud 2.4551E+01
Bnd Src/Visc Work: R1 to R1 Periodic 1 Side 1 -4.9736E-01
Bnd Src/Visc Work: R1 to R1 Periodic 1 Side 2 7.1249E-01
Bnd Src/Visc Work: R1 to S1 Side 1 4.2363E-01
Domain Src (Neg) : R1 -2.6186E+02
Domain Src (Pos) : R1 2.6186E+02
Domain Interface : R1 to R1 Periodic 1 (Side 1) -7.5045E+02
Domain Interface : R1 to R1 Periodic 1 (Side 2) 7.5044E+02
Domain Interface : R1 to S1 (Side 1) -5.4793E+03
-----------
Domain Imbalance : -7.1069E-02

+--------------------------------------------------------------------+
| H-Energy-S1 |
+--------------------------------------------------------------------+
Bnd Src/Visc Work: R1 to S1 Side 2 6.2703E-01
Boundary : S1 Inlet -7.9155E+02
Bnd Src/Visc Work: S1 Inlet 3.9477E-05
Bnd Src/Visc Work: S1 to S1 Periodic 1 Side 1 -8.8488E-01
Bnd Src/Visc Work: S1 to S1 Periodic 1 Side 2 5.7055E-01
Domain Src (Neg) : S1 -7.9268E+01
Domain Src (Pos) : S1 7.9268E+01
Domain Interface : R1 to S1 (Side 2) 7.9119E+02
Domain Interface : S1 to S1 Periodic 1 (Side 1) 1.2919E+01
Domain Interface : S1 to S1 Periodic 1 (Side 2) -1.2919E+01
-----------
Domain Imbalance : -4.9278E-02
In fact those numbers are pretty high. Is refining interfaces a solution?
Attached Images
File Type: jpg residuals.jpg (130.2 KB, 21 views)
sheaker is offline   Reply With Quote

Old   September 3, 2019, 08:53
Default
  #6
Senior Member
 
Join Date: Jun 2009
Posts: 1,873
Rep Power: 33
Opaque will become famous soon enough
If the concern is a mass imbalance, could you please post the flow summary for P-Mass in S1, and R1? The summary for R1 is missing.

Equation imbalances have nothing to do with mesh refinement. If the solution is converged, the imbalances must be satisfied. The solution may not be accurate because of a coarse grid; however, the imbalance must still be correct.
Opaque is offline   Reply With Quote

Old   September 3, 2019, 14:59
Default
  #7
Senior Member
 
sheaker's Avatar
 
Oskar
Join Date: Nov 2015
Location: Poland
Posts: 184
Rep Power: 11
sheaker is on a distinguished road
Thank Your for Your reply.
I think P-Mass for S1 and for R1 is in my previous post. Am I wrong? Are You talking about something different?

Quote:
[...]

+--------------------------------------------------------------------+
| P-Mass-R1 |
+--------------------------------------------------------------------+
Boundary : R1 Outlet -3.1318E-02
Domain Interface : R1 to R1 Periodic 1 (Side 1) 4.2258E-03
Domain Interface : R1 to R1 Periodic 1 (Side 2) -4.2258E-03
Domain Interface : R1 to S1 (Side 1) 3.1318E-02
-----------
Domain Imbalance : 1.6682E-07
[...]


[...]
+--------------------------------------------------------------------+
| P-Mass-S1 |
+--------------------------------------------------------------------+
Boundary : S1 Inlet 2.3842E-02
Domain Interface : R1 to S1 (Side 2) -2.3840E-02
Domain Interface : S1 to S1 Periodic 1 (Side 1) -3.8574E-04
Domain Interface : S1 to S1 Periodic 1 (Side 2) 3.8574E-04
-----------
Domain Imbalance : 1.4838E-06
[...]
Please, forgive me if I misunderstand something.




Best regards,
Oskar
sheaker is offline   Reply With Quote

Old   September 3, 2019, 17:18
Default
  #8
Senior Member
 
Join Date: Jun 2009
Posts: 1,873
Rep Power: 33
Opaque will become famous soon enough
My apologies.. I missed among all the others.

Two things come to mind:
1 - The solution is far from converged
2 - What pitch ratio are you running? The fact that there is such a change when crossing the rotor-stator interface is a sign of pitch change, and the mass flow must be scaled to match the full wheel conservation. Do I make sense

Mass flow at inlet = Mass Flow interface side inlet = Mass Flow interface side outlet = Mass Flow at outlet

Now divide by the number of passages in a given row, and multiply by the number of passage in the model. You should get your scaling factor.
Opaque is offline   Reply With Quote

Old   September 3, 2019, 17:58
Default
  #9
Senior Member
 
sheaker's Avatar
 
Oskar
Join Date: Nov 2015
Location: Poland
Posts: 184
Rep Power: 11
sheaker is on a distinguished road
Quote:
2 - What pitch ratio are you running? The fact that there is such a change when crossing the rotor-stator interface is a sign of pitch change, and the mass flow must be scaled to match the full wheel conservation. Do I make sense

I think this is a solution.
First improvement was to set inlet mass flow 'for component' not 'for entire machine'. Inlet mass flow became good after this. The issue with mass imbalance in between domains is most likely due to automatic pitch ratio and specific connection between stator and rotor domains. See attachment.



I think setting the right pitch ratio between domains will resolve the issue.



Thank You very much.
Attached Images
File Type: png GGI.png (16.6 KB, 25 views)
sheaker is offline   Reply With Quote

Old   September 3, 2019, 18:15
Default
  #10
Senior Member
 
Join Date: Jun 2009
Posts: 1,873
Rep Power: 33
Opaque will become famous soon enough
Your model seems to have a problem on the hub to shroud space.

I can see the pitch ratio scaling from periodic to periodic, but the lip/overlap at the bottom is a wall; therefore, the mass flow either do not fit, or you will get an artificial acceleration due to the area change.

Say you scale the model to full wheel, the cross-sectional area at the interface must be identical, correct? or your machine has a step/lip/contraction when at the crossing?
Opaque is offline   Reply With Quote

Old   September 4, 2019, 12:42
Default
  #11
Senior Member
 
sheaker's Avatar
 
Oskar
Join Date: Nov 2015
Location: Poland
Posts: 184
Rep Power: 11
sheaker is on a distinguished road
Thank You for Your reply.


Yes, in fact there is a step change in between rotor and stator domains. I think I need to set Pitch Ratio:

CFX Modelling Guide: "This option enables you to specify the pitch angle on side 1 and side 2 of the domain interface.You can use this option when the precise pitch angle is needed or when the area ratio of the two sides is not equal to the pitch ratio. In most cases, the automatic option is sufficient for determining the pitch ratio."


I need at least 2 days for simulation. I need time to check how those changes influences the final solution. Moreover I need to understand where and why my solution is far from convergence.
sheaker is offline   Reply With Quote

Old   September 4, 2019, 19:59
Default
  #12
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,854
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Quote:
Moreover I need to understand where and why my solution is far from convergence.
I commented on that in post #2. See FAQ: https://www.cfd-online.com/Wiki/Ansy...gence_criteria
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
ghorrocks is offline   Reply With Quote

Old   September 5, 2019, 09:09
Default
  #13
Senior Member
 
sheaker's Avatar
 
Oskar
Join Date: Nov 2015
Location: Poland
Posts: 184
Rep Power: 11
sheaker is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
I commented on that in post #2. See FAQ: https://www.cfd-online.com/Wiki/Ansy...gence_criteria



Thank You! That was what I meant by saying that I need to figure out how to improve convergence.
sheaker is offline   Reply With Quote

Reply


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
mass imbalance in Transient simulation sina_mech FLUENT 3 May 1, 2018 06:26
mass imbalance at the inlet, interface, outlet CFDturbomachinery CFX 2 October 30, 2014 12:54
mass imbalance problem iknowjesus FLUENT 0 September 2, 2014 04:22
mass imbalance in wind tunnel problem Robert2011 FLUENT 0 January 25, 2011 07:46
Residual Problem -Mass Imbalance Jie FLUENT 0 February 19, 2002 11:49


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