Re: [Yade-users] [Question #698778]: Problem with imposing flux in a cavity in flow engine
Question #698778 on Yade changed: https://answers.launchpad.net/yade/+question/698778 Status: Answered => Solved Zoheir Khademian confirmed that the question is solved: Thanks Robert Caulk, that solved my question. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698778]: Problem with imposing flux in a cavity in flow engine
Question #698778 on Yade changed: https://answers.launchpad.net/yade/+question/698778 Zoheir Khademian posted a new comment: Thanks so much Robert. After checking the VTK output for pressure, I noticed the cell I was injecting to was not registered as a cavity. So, simply applying the flux to one of the cavity cells in the cavity list solved the problem. Now, the influx and outflux are the same. One more question though: In your flux calculation above, the unit of q the flux seams to be m/s. Do we need to have the area of the cavity boundary to obtain the flux in m3/s? if so, what is the easiest way to get that area? Thanks again for the help! -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698830]: Can't running noFlowScenario.py with yade or yadedaily
Question #698830 on Yade changed: https://answers.launchpad.net/yade/+question/698830 Status: Answered => Open Ziyu Wang is still having a problem: Hello Robert! I have run the script oedometer.py on ubuntu 20.04 as you said, and everything is normal. And then,I have installed ubuntu 18.04 and ran noFlowScenario.py, unfortunately, it still did not succeed.There are different errors on yade and yadedaily: yade:On ubuntu 18.04, the version of yade can only use 2018.02b, but it seems that ThermalEngine is not defined on this version.So I install yadedaily. yadedaily:python3.6: /build/vtk6-VHOYAG/vtk6-6.3.0+dfsg1/IO/XML/vtkXMLOffsetsManager.h:142:void OffsetsManagerGroup::Allocate(int, int):Assertion `numElements > 0' failed. Aborted (core dumped) printAllVersions() in yadedaily as follows: Yade version : 20210923-5934~44b94a0~bionic1 Yade features : LOGGER USEFUL_ERRORS VTK OPENMP GTS QT5 CGAL PFVFLOW PFVFLOW LINSOLV MPI TWOPHASEFLOW FEMLIKE GL2PS LBMFLOW THERMAL PARTIALSAT POTENTIAL_PARTICLES POTENTIAL_BLOCKS Yade config dir: ~/.yadedaily Yade precision : 53 bits, 15 decimal places, without mpmath, PrecisionDouble ``` Libraries used : | library | cmake| C++ | | - | | --- | | boost | 106501 | 1.65.1 | | cgal | | 4.11| | clp | 1.16.11 | 1.16.11 | | cmake | 3.10.2 | | | coinutils | 2.10.14 | 2.10.14 | | compiler | /usr/bin/c++ 7.5.0 | gcc 7.5.0 | | eigen | 3.3.4| 3.3.4 | | freeglut | 2.8.1| | | gl| | 20190911| | ipython | 5.5.0| | | metis | | 5.1.0 | | mpi | 3.1 | ompi:2.1.1 | | mpi4py| 2.0.0| | | mpmath| 1.0.0| | | openblas | | OpenBLAS 0.2.20| | python| 3.6.9| 3.6.9 | | qglviewer | | 2.6.3 | | qt| | 5.9.5 | | sqlite| | 3.22.0 | | suitesparse | 5.1.2| 5.1.2 | | vtk | 6.3.0| 6.3.0 | ``` Linux version : Ubuntu 18.04.6 LTS Architecture : amd64 Little endian : True Looking forward to your reply! -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698388]: MPI with Clumps Particles
Question #698388 on Yade changed: https://answers.launchpad.net/yade/+question/698388 Status: Answered => Open Nima Goudarzi is still having a problem: Hi Bruno, Thanks so much. I think both options are tedious but I am interested in trying 1 first. Are there any examples showing me how to do this for clumps? Do I need to predefined subdomains manually (I mean grouping clumps) before importing or I am able to do so in the main analysis script? For 2, I don't have access to hundreds of cores, unfortunately. I highly appreciate it if I'd receive some guidance for 1. Much Obliged Nima -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698754]: Making video of pfacet clump: clump falling apart
Question #698754 on Yade changed: https://answers.launchpad.net/yade/+question/698754 Status: Open => Answered Bruno Chareyre proposed the following answer: Hi, this ugly hack works for me. Depending on the latency of the OS you may have to adjust the sleep() times. I'm considering adding a function to redraw, since here I'm just indirectly calling redraw by pausing. This is not really using makeVideo but you can merge the files into one video indepently, or find inspiration in the makeVideo code. I hope it helps. Bruno snapN=0 def snap(): global snapN O.pause() time.sleep(1) yade.qt.views()[0].saveSnapshot("snap"+str(snapN)+".png") time.sleep(1) snapN+=1 if snapN<5: O.run() O.engines = O.engines+[PyRunner(command="snap()", iterPeriod=100)] -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698836]: temp=Nan?
Question #698836 on Yade changed: https://answers.launchpad.net/yade/+question/698836 Robert Caulk posted a new comment: Until you can properly execute the example scripts on 20.04, it is futile to try to debug unknown scripts such as this. Therefore, I propose we fix your other threads first and then come back to this one. Thanks, Robert -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698830]: Can't running noFlowScenario.py with yade or yadedaily
Question #698830 on Yade changed: https://answers.launchpad.net/yade/+question/698830 Status: Open => Answered Robert Caulk proposed the following answer: Hello, Ok, I suspect ThermalEngine has some issues on 20.04. But first can you please run oedometer.py [1] and report your results? Next, can you install ubuntu 18.04 and run the thermalengine sample script? The two quickest options for this are to install ubuntu 18.04 on a virtual machine or to boot 18.04 from a USB. At least this will verify to us that ThermalEngine has issues on 20.04. Finally, if you want to try some alternative installations on 20.04 you can: * try running docker [2] * try compiling from source [3] At the moment, I have not run ThermalEngine on Ubuntu 20.04 so I cannot bug test this directly. But I will try to spin up a virtual machine or boot from USB in the coming week. Cheers, Robert [1]https://gitlab.com/yade-dev/trunk/-/blob/master/examples/FluidCouplingPFV/oedometer.py [2]https://yade-dem.org/doc/installation.html#docker [3]https://yade-dem.org/doc/installation.html#source-code -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698778]: Problem with imposing flux in a cavity in flow engine
Question #698778 on Yade changed: https://answers.launchpad.net/yade/+question/698778 Robert Caulk posted a new comment: Also, I notice that you define the cavity cells *after* emulating action in flow and setting flow.dead=0. Please do this *before* to avoid any kind of undefined behavior in the solver. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698778]: Problem with imposing flux in a cavity in flow engine
Question #698778 on Yade changed: https://answers.launchpad.net/yade/+question/698778 Status: Open => Answered Robert Caulk proposed the following answer: Hey Zoheir, >> I injected 1e-4 m3/s into the cavity but the output of flow.getCavityFlux() [1] is about 8e-18 m3/s. Would you please let me know what I am missing? Before I dive into the script myself (I won't have time until the weekend), I will propose some possible reasons. You are using a compressible fluid, and quite a small time step with very few time steps before checking the cavity flux. So technically, the governing equations would allow for an unequal fluid flux exiting the cavity to that which is entering the cavity. However that contrast in fluxes seems quite stark and so it is very possible something else is wrong. Have you saved the VTK output to verify that the cavity is defined as you expect it to be? This will also allow you to look at the pressure gradients and verify that they are also as you expect them to be. You can also compute the flux yourself by iterating on the cavity cells. This would verify whether or not the function "getCavityFlux()" is correct. (btw it is useful to store the list of cavityCells by their ID so you can just iterate on those directly eachtime you want to make a calculation such as this). totalVolume = 0 v = np.array([0,0,0]) cavityCells = [] for i in range(0,flow.nCells()): coords = flow.getCellCenter(i) if flow.getCellPImposed(i): continue if np.linalg.norm((coords-(maxX-dx/2.,maxY-dy/2.,maxZ-dz/2.)))> By the way, does flow.imposeFlux apply volumetric flux or mass flux? In FlowEngine, it is a volumetric flux. Cheers, Robert -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp