Hello Team, I installedMFiX-19.3.1 on Ubuntu 18.04.5 LTS.
Firstly, when I opened GUI, I got below error which I could troubleshoot using few answers on this forum.
Could not import VTK, please check your installation.
traceback:
Traceback (most recent call last):
.
.
.
.
When I open a 3D TFM tutorial case and click run, I end up into this strange error as attached. I cannot stop this nor kill this run: I have reboot my computer to start afresh.
Please attach the output of conda list. Or better yet, go to the main menu and click Submit bug report, this will create a ZIP file containing debugging information. I haven’t seen this error before.
Hello @cgw, I have an update. Error I posted above shows up when I click “Run” on GUI on my laptop. (I will send you conda list and debug zip in few minutes).
When I submit same case on cluster. After building custom solver, case gets into iterating and I do not see any error written to output file. But strangely, bed wont fluidize for any gas velocities. Residuals just move around E-8 as shown below.
Just to be clear, it looks like there are two separate issues now? You can’t run at all on your laptop, but you got MFiX running on a cluster… but the results on the cluster are different from what you expected - is that correct?
A few points:
I see that you are running MFiX version 19.3.1. If possible, please upgrade to version 21.4. It is difficult for us to find and fix bugs in older versions, especially if those bugs are not easily reproducible. If this bug is still present in the current version, we will investigate further.
Your report starts with “Could not import VTK”, but the traceback was omitted. That information can be helpful for debugging. Sometimes VTK issues are due to older hardware not supporting modern GL shaders, which newer versions of VTK require. If you are installing on an old laptop you may need to use something like “vglrun”
Will follow up separately regarding the cluster run.
Hi @cgw, at the moment, let us focus only on local machine installation. (At least I’ll have something to use as I’m chasing a deadline).
I started all-over again to install latest version available which is MFiX 21.4.
Towards end of “sh Anaconda3-2021.11-Linux-x86_64.sh” step, I see syntax error as shown.
Software gets into iterations without any errors. [I’d update you about results very soon].
3a. I’d setup this version on cluster and see if this works.
Case I’m testing has a point probe on centroid of YMIN region. There’s a monitor “Average” type to calculate pressure and y-velocity, besides, point probe monitor.
These both point to permissions problems, ie using sudo to install. Since this is not an MFiX-specific problem I’m afraid I can’t provide further insight on this one. It may not even be a problem (?)
conda output looks fine. But in the second screenshot, I see the message To launch MFIX: MFIX
This is not coming from our installation, this must have been added by someone at your university who previously installed MFiX. I can’t guarantee that your installation and this installation are not going to interfere, you may have to get help from a local systems person. However, if you have your own Anaconda installation activated, that should take priority over whatever is previously installed.
It runs! Hooray.
Please attach case files (“Submit bug report” from main menu) if you want us to diagnose output.
Hello Charles @cgw, please find attached zip folder to address monitor(s) issue. I tried “Average” and “Area-Weighted Average” and both are failing as shown below.
My observation is point probe is writing superficial gas velocity at its location, which is one-cell above YMIN face and other two monitors are writing bc_v_g(1) from usr1.f I think, to estimate UMF or for fluidization studies, superficial gas velocity is more appropriate. If we can get average across boundary instead of point probe, it is more acceptable
As a reminder, I’m using TFM on rectangular domain (please see above, no cut-cell, 2 cells thick).
A. VTK files, say BG_0010.vtu, are saved in binary. Can it be saved as ASCII or txt or csv?
B. How to write txt or csv table with x, y, z (of cell centroids) and pressure, voidfraction....? [I read a post where Paraview exports each vtu file into csv but I have around 400 vtu(s) to work on].
Why are you saying that “Average” and “Area-Weighted Average” are both failing? They give you the correct gas velocity (as specified in usr1.f) at y=ymin plane. This is a BC value. I am not sure what value you were expecting.
The value in BOTPOINT.csv is not superficial velocity, it is just the gas velocity at one point in the flow field. It won’t be the same as the BC value.
Can a 2 mm distance (between YMIN AVERAGE and POINT PROBE) has this different gas velocity?
Yes. Again, you should not compare the inlet velocity (where void fraction is one) with the gas velocity in the bed (where the void fraction is less than one). For example if the inlet velocity is 0.7m/s and the void fraction in the bed is 0.5, then the gas velocity at that point won’t be 0.7m/s but close to 0.7/0.5=1.4 m/s due to the volume occupied by solids.
If YMIN AVERAGE is working correct, why pressure is always 101320, when bed is vigorously bubbling?
YMIN AVERAGE gives you the velocity at y=ymin, but it correspond to the ghost cell just below y=ymin because the y-velocity is stored at the north face of the cell (staggered grid representation). The pressure is stored at the cell center. The value 101320 is not applied but used to compute gas density. The pressure you want to use for a pressure drop calculation is the pressure at the first fluid cell above the inlet and the pressure at the outlet (outflow pressure).
Zero pressure drop across bed as this is PressureOutflow value?
See answer to 2)
A. It is interesting to note that void fraction at YMIN plane is always 1.0. This confirms that imposed boundary condition value, bc_v_g(1), in GUI or UDF, is superficial gas velocity of system.
Also, this value being a constant, in my case, it is uniform across entire patch / boundary. So, average/value type(s) monitor is same.
B. How can we output an average of pressures at all first fluid cells above inlet? If I understood staggered cell arrangement correctly, red line shows the plane where this pressure is to be evaluated for pressure-drop-curve. Please see attached image below and ignore probe representation in this figure.
Hello Dr.@jeff.dietiker, @cgw, @onlyjus, I think, I figured out procedure to extract pressure-drop across bed.
I created a plane above YMIN (gas inlet) such that its position is around 75% of first-cell-height. For lower positions, say 50% or 25%, average-pressure-monitor showed e+32 garbage values.
With superficial gas velocities directly obtained from usr1.f, we see following plots.
In my opinion, minimum fluidization velocity is 0.316 m/s, occurs at 11.57s from start. Could you please confirm this data extraction from MFiX is correct?
Best,
