|
[Sponsors] |
Eulerian-Eulerian turbulent two-fluid model is not converging? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 16, 2020, 14:09 |
Eulerian-Eulerian turbulent two-fluid model is not converging?
|
#1 |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
Ansys/Fluent/Multiphase
Hello All, I am using the Eulerian-Eulerian two-fluid model (in Ansys Fluent) to simulate the turbulent solid-liquid two-phase flow in a horizontal pipe. This is the reference paper I am trying to reproduce: Eulerian-Eulerian Simulation of Particle-Liquid Slurry Flow in Horizontal Pipe Titus Ntow Ofei and Aidil Yunus Ismail Currently, I am quite sure about my physical model settings. But my simulation is not converging and please see the attachment. I think the settings that I can play with are Solution Controls, Method Controls, and mesh. In addition, except for playing with the current model, I am also thinking to first simulate a laminar flow and then add the turbulence model. Your valuable comments would be greatly appreciated! |
|
April 16, 2020, 21:53 |
Updates
|
#2 |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
For the current model:
1. I reduced the inlet velocity to have a smaller Reynolds number (but it is still turbulent flow), it didn't converge. 2. I used a very small inlet velocity to have a laminar flow, it didn't converge. Then, I thought it might be due to the mesh quality. 3. I even established a new model on another reference paper. It seems for me it won't converge (please see the attachment). In this case, I think the mesh should be fine. This is the reference: Sand Transport and Deposition Behaviour in Subsea Pipelines for Flow Assurance. Yan Yang, Haoping Peng and Chuang Wen. 4. In addition, I saw some Console outputs before crashing: Reversed flow on 2 faces of pressure-outlet 5. Stabilizing epsilon to enhance linear solver robustness. Stabilizing epsilon using GMRES to enhance linear solver robustness. turbulent viscosity limited to viscosity ratio of 1.000000e+05 in 465776 cells Currently, I am thinking I might need to degrade my case further to ONLY simulate turbulent fluid flow first, then add the solid phase. Your valuable comments would be very appreciated! Last edited by minzhang; April 17, 2020 at 02:01. |
|
April 17, 2020, 02:03 |
Updated not converging results
|
#3 |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
Please see the attachments of both the Graphics and Console output.
Your comments would be very appreciated! |
|
April 17, 2020, 06:52 |
Convergence
|
#4 |
Senior Member
|
The image you showed in the second last post shows convergence, however, the latest one does seem to have some issues. Laminar cases are more difficult to converge than the turbulent ones since the momentum does not diffuse and cause numerical instabilities that require higher order models for resolution. What you call non-convergence is the expected behavior for a transient case. The simulation is supposed to converge in each time-step, hence, the residuals appear to oscillate. In reality, these are convergence plots for each time-step joined together. Hence, they appear as not converging.
__________________
Regards, Vinerm PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority. |
|
April 17, 2020, 13:00 |
|
#5 | |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
Thank you so much for your reply!
How about the Console outputs? It seems that the epsilon equation blows up and also the turbulent viscosity calculation has some issues? Quote:
|
||
April 17, 2020, 13:21 |
Convergence
|
#6 |
Senior Member
|
As per your last post with the images, yes, the solver is facing difficulty. However, this could be due to anything, starting from a bad mesh all the way down to numerical schemes. You have to ensure that the mesh is good. All material properties, boundary conditions, phase interactions, and operating conditions are good. Then check numerical setup, such as, discretization schemes, pressure-velocity coupling, and URF, etc. You may also try a smaller time-step. If the simulation requires more than 30-40 iterations per time-step for convergence, then your time-step is large for the setup.
__________________
Regards, Vinerm PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority. |
|
April 19, 2020, 17:08 |
Try to reproduce the results in the paper
|
#7 |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
More info:
For the first reference paper: The case I am trying to reproduce is: particle diameter = 270microns, particle volume fraction = 0.1. The paper was using CFX and I am using Fluent. The settings I am using are consistent with the settings mentioned in the paper. The settings that the paper did not mention and I think they might influence my simulation results are listed below. In addition, there is a difference between their geometry and mine. In their paper, L=44D, and cell no. = 801,848. Since I am using the student free version and I have the limit of 512,000 cell number. So in my case, L=25D and my cell no. = 508,600. My mesh might be different from theirs but I think my mesh resolution should be fine compared with theirs. Please see the attached screenshot of my mesh. My Scaled Residuals plot is also attached. Your valuable comments would be very very appreciated!! Settings were not mentioned in the paper: 1. Materials - Fluid - solid: what is the viscosity value? (1e-6 kg/(m.s)) 2. Multiphase Model - Phases - Secondary Phase - Granular Conductivity: syamlal-obrien 3. Multiphase Model - Phases - Secondary Phase - Frictional Viscosity: shaeffer 4. Multiphase Model - Phases - Secondary Phase - Frictional Pressure: Johnson-et-al 5. Multiphase Model - Phases - Secondary Phase - Radial Distribution: lun-et-al 6. Multiphase Model - Phase Interaction - Forces - Drag Coefficient: gidaspow 7. Viscous Model - Turbulence Multiphase Model: Per Phase 8. Velocity Inlet - Turbulence - Turbulent Intensity: 3 Velocity Inlet - Turbulence - Hydraulic Diameter: 0.0072 Velocity Inlet - Granular Temperature: 0.007272 9. Pressure Outlet - Backflow Turbulent Intensity: 0 Pressure Outlet - Backflow Turbulent Viscosity Ratio: 0 Pressure Outlet - Backflow Granular Temperature: 0 Pressure Outlet - Backflow Volume Fraction: 0 10: Wall - Granular Conditions: Johnson-Jackson Wall - Granular Properties - Restitution Coefficient: 0.2 |
|
April 20, 2020, 04:27 |
Multiphase Model
|
#8 |
Senior Member
|
If volume fraction is only 0.1, you should not use Granular model. I suppose you have got the number wrong because 0.1 particle volume fraction does not imply a slurry. So, I'd suggest you verify this number. And if it is really less than 10%, use DPM; that's more accurate.
__________________
Regards, Vinerm PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority. |
|
April 20, 2020, 12:04 |
Volume fraction
|
#9 | |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
Hello vinerm,
The particle volume fraction is 10%-45%. 0.1 means 10%, is that right? Quote:
|
||
April 20, 2020, 12:32 |
Volume Fraction
|
#10 |
Senior Member
|
Yes, then granular makes sense. Instead of using per phase turbulence, use mixture. Do not set flow reversal parameters to 0. If there is flow reversal in the beginning, it will cause trouble. Only pressure should be set to 0. For others, use some positive values.
__________________
Regards, Vinerm PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority. |
|
April 20, 2020, 13:58 |
How to set the backflow parameters?
|
#11 |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
I am wondering how to decide the backflow parameters? What would these non-zero backflow parameters do?
Thanks! |
|
April 20, 2020, 16:30 |
Flow Reversal
|
#12 |
Senior Member
|
If there is no flow reversal at the outlet for you case, then these values are not used. However, if there is a flow reversal at the outlet, then the outlet begins to behave like inlet. And for inlet, Fluent requires inlet conditions. So, you can use conditions similar to you inlet.
__________________
Regards, Vinerm PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority. |
|
April 20, 2020, 16:38 |
|
#13 |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
If the calculated velocity at one outlet cell is -1m/s, then the outlet boundary condition will adjust this value to be backflow parameter, for example, +5m/s. Is that correct?
|
|
April 20, 2020, 16:59 |
Calculation at Outlet
|
#14 |
Senior Member
|
The velocity for the backflow does not depend on the velocity of the outgoing flow. It depends on the pressure in the cell adjacent to the outlet. If the pressure reduces below the pressure at the outlet, then the flow starts to reverse. Flow rate or velocity will depend upon the difference in the pressures at the outlet and first cell adjacent to the outlet.
__________________
Regards, Vinerm PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority. |
|
April 20, 2020, 18:39 |
|
#15 |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
If the pressure at the first cell adjacent to the outlet reduces below the pressure at the outlet, then the flow starts to reverse.
Then, what will happen if I have a backflow condition which is the same as the inlet condition? Thanks! |
|
April 21, 2020, 01:39 |
|
#16 |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
I made these changes, but it still could not converge.
|
|
April 21, 2020, 05:16 |
Pressure at the outlet
|
#17 |
Senior Member
|
When I mentioned that you can maintain similar conditions as inlet, I did not mean the pressure as well. Until unless the pressure is lower than the inlet, flow cannot pass though. Pressure at the inlet should be left at 0. For turbulence quantities, you can use similar values. Volume fraction of secondary phase can be left as 0 as well.
__________________
Regards, Vinerm PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority. |
|
April 21, 2020, 11:37 |
|
#18 | |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
I gave the "Supersonic/Initial Gauge Pressure" = 0 for my Velocity Inlet. I think this is only an initial pressure value for the velocity inlet. The solver will calculate the pressure at the inlet, which will be higher than 0.
For my Pressure Outlet, the Gauge Pressure = 0. Quote:
|
||
April 21, 2020, 11:42 |
Correct
|
#19 |
Senior Member
|
Yes, that's correct.
__________________
Regards, Vinerm PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority. |
|
April 21, 2020, 11:49 |
|
#20 |
Member
Min Zhang
Join Date: Mar 2017
Posts: 81
Rep Power: 9 |
Since I still could not have a converged result, the only option I have now is to have a refined mesh. Is that right?
In addition, I am wondering what are the possible reasons that I could not reproduce that reference paper? Thanks! |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] swakExpression not writing to log | alexfells | OpenFOAM Community Contributions | 3 | March 16, 2020 19:19 |
Velocity vector in impeller passage | ngoc_tran_bao | CFX | 24 | May 3, 2016 22:16 |
Multiphase simulation of bubble rising | Niru | CFX | 5 | November 25, 2014 14:57 |
Setting rotating frame of referece. | RPFigueiredo | CFX | 3 | October 28, 2014 05:59 |
Terrible Mistake In Fluid Dynamics History | Abhi | Main CFD Forum | 12 | July 8, 2002 10:11 |