|
[Sponsors] |
Why does velocity increase unexpectedly in density based solver for compressible flow |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 13, 2021, 11:10 |
Why does velocity increase unexpectedly in density based solver for compressible flow
|
#1 |
New Member
stalukar
Join Date: Aug 2016
Posts: 13
Rep Power: 9 |
I am trying to simulate a mixing process of Helium jet in atmospheric air. He is coming out from a tube of 4.6 mm inner dia with a Re number of 2500 where the air is considered stagnant.
To set the model, I use an inlet velocity BC with a velocity of 66.5 m/s for He jet inflow in the domain. I also set a zero pressure outlet BC at a sufficient distance from the inlet to allow enough space for mixing. The model is in a steady state. To capture the mixing, I use the species transport model available in Fluent. It looks like the model works fine. However, I found the velocity goes unexpectedly high (~1200 m/s max) in the domain when the solution is converged. If anyone has any suggestions, it would be highly appreciated. |
|
April 13, 2021, 11:19 |
|
#2 |
Senior Member
Sayan Bhattacharjee
Join Date: Mar 2020
Posts: 495
Rep Power: 8 |
Wouldn't the zero pressure outlet cause the air in the duct to accelerate out?
The setup seems to be basically a rocket exhaust nozzle in space vacuum. High exhaust velocity makes sense for such a rocket nozzle. Though, I'm not sure about 1200 m/s. Am I misunderstanding something? Last edited by aerosayan; April 13, 2021 at 11:21. Reason: fix |
|
April 13, 2021, 11:27 |
|
#3 |
New Member
stalukar
Join Date: Aug 2016
Posts: 13
Rep Power: 9 |
aerosayan,
Thanks for your reply. Zero Pressure outlet: that's a valid point. Though the velocity profile looks like rocket propulsion, the physics of my case is simple - helium is coming from a tube like a jet and being mixed into the air. To my understanding, the velocity at the outlet should be much smaller than the inlet velocity as the flow decelerates as it gets diffused to air. In that case, could you please suggest any other BC type for the outlet? |
|
April 13, 2021, 11:37 |
|
#4 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71 |
If I understand, the model is for full compressible flows, isn't it?
With such a high velocity the flow is supersonic and the outlet condition does not enter in the computation What about the other inflow conditions? Being subsonic you need to fix all but one boundary condition. Have you verified the other variable? |
|
April 13, 2021, 11:45 |
|
#5 |
New Member
stalukar
Join Date: Aug 2016
Posts: 13
Rep Power: 9 |
the model is for full compressible flows, isn't it?
Yes The inlet velocity is only 66.5 m/s (M ~ 0.2), and I am not sure why it is becoming supersonic as it moves forward. I am not sure if the inflow condition is the crucial factor. Have you verified the other variable? I am not sure what you meant. I made sure the model setup was okay except the BCs as the results are not expected. |
|
April 13, 2021, 12:03 |
|
#6 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71 |
Quote:
First of all, your simulation should be unsteady in a deterministic sense, thus are you using a RANS formulation? In compressible flows you need to prescribe the inlet condition for the full set of equations, that is density, momentum and energy. For subsonic flows you have to let a condition coming from the interior, you cannot prescribe all conditions in terms of Dirichlet. Hence, what about the prescription of the other variables? |
||
April 13, 2021, 12:10 |
|
#7 | |
New Member
stalukar
Join Date: Aug 2016
Posts: 13
Rep Power: 9 |
Quote:
If you could write a bit detail, it would be really helpful for me. |
||
April 13, 2021, 12:33 |
|
#8 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71 |
Quote:
You need to know the theory of the chararacteristic for hyperbolic equations. As example, consider a 2D case, you have to consider the conditions for 4 variables at the inlet, that is density, x-velocity, y-velocity, total energy. For subsonic flow the characteristic line related to the eigenvalue u-a will come from the interior towards the inlet section. Therefore you have to set one variable free to come from the interior. Have a look to the theory in some textbook. |
||
April 14, 2021, 02:11 |
|
#9 | |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,274
Rep Power: 34 |
Quote:
Are you using pressure based solver or density based solver? |
||
April 14, 2021, 11:18 |
|
#10 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71 |
||
April 15, 2021, 01:18 |
|
#11 |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,274
Rep Power: 34 |
Ohh thank you. I did not carefully read the title.
I recently implemented two phase comparessible flow model and did some similar calculations like injecting a gas etc and I did not observe this type of behaviour. But I did see this type of behaviour with pressure based solver. So i find it bit strange the density based solver shows this. |
|
April 15, 2021, 02:57 |
|
#12 | |
Senior Member
Join Date: May 2012
Posts: 548
Rep Power: 15 |
Quote:
I would not use an inlet velocity boundary condition. Pressure inlet or mass flow inlet should be the proper one. |
||
April 15, 2021, 04:18 |
|
#13 |
Senior Member
|
Could you give more details on the relevant settings? 2D/3D/AXI? What are all the material settings? All the other boundary conditions? Maybe a sketch of the mesh/domain? Convective scheme? Gradient? Turbulent or laminar?
Also, to which extent is this high speed region relevant? Does it seem as a physical result of the used boundary conditions or some numerical artifact? Maybe a picture would help as well. Finally, qualify converged results. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Issues on the simulation of high-speed compressible flow within turbomachinery | dowlee | OpenFOAM Running, Solving & CFD | 11 | August 6, 2021 06:40 |
Pressure based and Density based Solver | Xobile | Main CFD Forum | 39 | August 19, 2020 06:04 |
New solver with temperature-dependent density based on simpleFoam? | fanta | OpenFOAM Programming & Development | 0 | December 5, 2017 03:40 |
Pressure based and Density based Solver | taekyu8 | FLUENT | 0 | January 28, 2013 11:05 |
Density based solver with cnst density! | p_agoodboy | FLUENT | 0 | November 8, 2010 22:18 |