CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

3D Computation through a radial impeller

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

LinkBack Thread Tools Display Modes
Old   December 15, 1999, 11:52
Default 3D Computation through a radial impeller
Posts: n/a

I'm using a 3D solver to compute the flow through a radial impeller passage. To make the computation, I have the choice between 2 numerical schemes: a second order central sheme with second and fourth order numerical dissipation and a first order upwind scheme.

The problem is that I have convergence problem when I use the central scheme. The massflow and all other values oscillate.

I find a solution with the upwind scheme but the massflow seems to be smaller than with the other scheme. I asked myself whether the massflow difference between the two shemes comes from the numerical viscosity of the first order scheme. What do you think about that?

Do the oscillations I found using the second order scheme come from a supersonic region?

I've checked and in certain regions the velocity is (M=0.8 1). Near impeller entry and in diffusor outlet.

Thanks for your answer

  Reply With Quote

Old   December 15, 1999, 13:36
Default Re: 3D Computation through a radial impeller
John C. Chien
Posts: n/a
(1). I think you should be very happy with the code you use, because it can produce different answers using different options. (2). I would say that it is likely different user will also get different answers. (3). That is fairly consistent. (4). For the person who is interested in saving money, he can use the upwind option. (5). For the person who is interested in DNS, the transient nature of the second-order central difference solution would be very attractive. (6). So, there is no guarantee that the second-order method with a second-order and a fourth-order artificial viscosity terms is better than the first-order difference. (7). This is important, because when you run into the shock waves, it is likely that everything will be down graded to the first-order accuracy. (8). So, the conclusion is, they are not supposed to be the same, and you are all right. (9). Well, about the relationship between the oscillation and the supersonic region, you can easily find the answer by lowering the Mach number in the flow field. (by simply increasing the inlet temperature will lower the Mach number) (10). It is hard to know the answer unless you actually try it out. (11). The oscillation itself is not a bad thing though, because some people like to think that they are getting higher-order solutions. For them, a converged first-order upwind solution is just like venela icecream, it is not very exciting.
  Reply With Quote

Old   December 15, 1999, 15:49
Default Re: 3D Computation through a radial impeller
Adrin Gharakhani
Posts: n/a
>For the person who is interested in DNS, the transient nature of the second-order central difference solution would be very attractive.

Huh? What does this mean?

As regards to losing mass when using upwind, I'm not sure what you mean by "losing" mass; however, if you have a (consistent) scheme that is conservative then you shouldn't be losing mass, no matter what order differencing you use. You may/will get a wrong (unconverged) velocity profile, but its integral (of course including density effects) should satisfy continuity. Theoretically, the upwind scheme introduces (numerical) diffusion. And all diffusion is supposed to do is to spread the diffused variable out. Period. It shouldn't create or destroy mass!

As for the oscillations with central differencing - these are quite famous. Without knowing the details of your code and simulation, I can guess that your CFL condition is not satisfied properly. That is, you are using too coarse a grid or too large a timestep. A simulation at a finer level should reduce the oscillations. If you (can) reduce your Mach number, you "may" see a correction to your problem as well. That's because you have effectively changed your CFL condition.

Adrin Gharakhani
  Reply With Quote

Old   December 15, 1999, 20:56
Default Re: 3D Computation through a radial impeller
John C. Chien
Posts: n/a
(1). I shold have used the words "less attractive, relative to the upwind method". (2). And I am sure that people in DNS use much higher order methods, or spectral method. The solution is forced into a series of waves and oscillations as well, long before one even started the solution procedure.
  Reply With Quote


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
how to extract the components of radial and tangential velocities from pump impeller renyun0511 OpenFOAM Running, Solving & CFD 2 November 1, 2011 23:09
problem in the CFX12.1 parallel computation BalanceChen ANSYS 2 July 7, 2011 10:26
Setting Pressure outlet boundary for simulation of Radial impeller Mukund Pondkule Main CFD Forum 1 April 5, 2011 09:24
Two phase flashing flow through radial expander impeller Allan OpenFOAM 0 April 16, 2009 01:00
Computation through radial impeller Ursenbacher Main CFD Forum 5 January 12, 2000 02:27

All times are GMT -4. The time now is 07:27.