Re: [Yade-users] [Question #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Robert Caulk posted a new comment: Hello, Please either indicate which answer solved your problem by marking it with a star, or describe the solution to your problem. Thank you for helping keep an organized knowledge base. -rc -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Answered => Solved Ziyu Wang confirmed that the question is solved: That solved my problem. -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Expired => Answered Robert Caulk changed the question status: Marking as answered to avoid launchpad janitor until OP returns. -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Needs information => Expired Launchpad Janitor expired the question: This question was expired because it remained in the 'Needs information' state without activity for the last 15 days. -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Expired => Needs information Robert Caulk changed the question status: Avoiding launchpad janitor -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Robert Caulk posted a new comment: Hello, Please confirm that updating your yade installation has solved your problem. If not, please indicate the exact changes you have made compared to the example script, and which change in particular is the one that results in your issue. Thank you, 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Open => Expired Launchpad Janitor expired the question: This question was expired because it remained in the 'Open' state without activity for the last 15 days. -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Solved => Open Ziyu Wang is still having a problem: Hello Robert, The situation is like this, I consulted the forum and referred to the solution of [1] to change the value of flow.useSolver. As you said,when I change the value to 0 or 3,it seems that the error not happen and everything is normal. I also refer to the explanation of flow.useSolver in [2] (ps:Maybe it’s because I don’t know what the exact result should look like) Yes,the script have thermal.advection=True for achieving advection,but when I run the script as useSolver=4,it still have the same problem although I have update the yadedaily to 20211018-5981~9b2e38e.: CHOLMOD warning: matrix not positive definite. file: ../Supernodal/t_cholmod_super_numeric.c line: 911 the plot is blank and O.bodies.state.temp=nan.. What should I do. Best regards. [1]https://answers.launchpad.net/yade/+question/697642 [2]https://yade-dem.org/doc/yade.wrapper.html?highlight=flowengine#yade.wrapper.FlowEngine.useSolver -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Robert Caulk posted a new comment: Hello, Although you have not specified what value you changed flow.useSolver to, I can tell you that anything besides 4 will not solve for advection. Hence why it is important to make sure the example scripts run as is. 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
Re: [Yade-users] [Question #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Ziyu Wang posted a new comment: Just an update: Hello everyone,I seem to have solved this problem by changing the value of flow. Usesolver(Although I'm not sure if it's the right thing to do) Best regards. -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Answered => Solved Ziyu Wang confirmed that the question is solved: Hello, It seems that the situation is very bad for me personally... Ok,I will try to solve this problem myself.. Finally, how should CHOLMOD warning be solved. Thank you for your timely and patient help again! -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Open => Answered Robert Caulk proposed the following answer: Hello, noFlowScenario.py works as expected on 18.04 using latest source code. I ran the script you provided, it also works. Since I am unable to reproduce the error that you describe, I cannot help you. 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
Re: [Yade-users] [Question #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Needs information => Open Ziyu Wang gave more information on the question: Hello, I am sorry to tell you that I do use latest source compiled on 18.04(This can be verified by the information in printallversions() ). The only change is that: shutil.copyfile(sys.argv[0],'txt'+timeStr+identifier+'/'+sys.argv[0]) (FileNotFoundError: [Errno 2] No such file or directory,So I delete this line,but this should not lead to problems?) In fact,I guess the root of the problem lies in 'CHOLMOD warning'?(Just guess) >>Yes it means you have no spheres touching one another, I rechecked thermoHydroMechanical_coupling.py and found that after makecloud, triaxial internal compaction is carried out to make the sphere in close contact.Thank you for your explanation. I have realized my mistake. (ps:I seem to remember that during the compilation and installation process (i.e. during the execution of the 'make' command), the terminal prompts warnings about Law2_ScGeom_MindlinPhys_Mindlin?,does it have anything to do with this) Best regards. -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Open => Needs information Robert Caulk requested more information: Hello, Please confirm that you have latest source compiled on 18.04 and noFlowScenario.py (no modifications) does not run. If you confirm this, then there is something else going on which may be difficult to trouble shoot. I am running 18.04 on latest source and noFlowScenario.py runs as expected. >>does the gaseous state of the sphere mean that the sphere is too scattered Yes it means you have no spheres touching one another, which would imply quite large conduction diffusivity, which would imply very small timestep. But besides this, the conduction model is not physically realistic in such a gaseous state so it should not be used. You can view your packing by clicking "Show 3D" in the GUI window. 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
Re: [Yade-users] [Question #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Needs information => Open Ziyu Wang gave more information on the question: Hello,Robert 1.I have followed your suggestion, installed ubuntu18.04 and installed Source-code.Unfortunately,the problem still exists..: CHOLMOD warning: matrix not positive definite. file: ../Supernodal/t_cholmod_super_numeric.c line: 911 /usr/lib/python3/dist-packages/matplotlib/__init__.py:831: MatplotlibDeprecationWarning: axes.color_cycle is deprecated and replaced with axes.prop_cycle; please use the latter. mplDeprecation) /usr/lib/python3/dist-packages/matplotlib/__init__.py:801: MatplotlibDeprecationWarning: axes.color_cycle is deprecated and replaced with axes.prop_cycle; please use the latter. mplDeprecation) printAllVersions(): Yade version : 2021-10-08.git-a6d7547 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: ~/.yade-2021-10-08.git-a6d7547 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 | | sphinx| 1.6.7-final-0| | | 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 (ps:other example scripts such as flowscenario and THMcoupling can run normally,only this script is unsuccessful, maybe it is my personal installation problem?) 2.>>Have you actually looked at your packing? It seems you are simulating a gaseous state of spheres (although there is no way to know without an MWE). ThermalEngine was not designed to handle this kind of packing. Sorry for my stupidity,does the gaseous state of the sphere mean that the sphere is too scattered (not dense enough?).I refer to thermoHydroMechanical_coupling.py before using makecloud. Is there anything wrong? MWE as follows: # from yade import pack, ymport from yade import timing import numpy as np import shutil timeStr = time.strftime('%m-%d-%Y') num_spheres=1000# number of spheres young=1e9 rad=0.003 mn,mx=Vector3(0,0,0),Vector3(0.05,0.05,0.05) # corners of the initial packing thermalCond = 2. #W/(mK) heatCap = 710. #J(kg K) t0 = 333.15 #K r = rad k = 2.0 # 2*k*r Cp = 710. rho = 2600. D = 2.*r m = 4./3.*np.pi*r**2/rho identifier = '-noFlowScenario' if not os.path.exists('VTK'+timeStr+identifier): os.mkdir('VTK'+timeStr+identifier) else: shutil.rmtree('VTK'+timeStr+identifier) os.mkdir('VTK'+timeStr+identifier) if not os.path.exists('txt'+timeStr+identifier): os.mkdir('txt'+timeStr+identifier) else: shutil.rmtree('txt'+timeStr+identifier) os.mkdir('txt'+timeStr+identifier) mn,mx=Vector3(0,0,0),Vector3(0.05,0.05,0.05) O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(3),density=2600,label='spheres')) O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='walls')) walls=aabbWalls([mn,mx],thickness=0,material='walls') wallIds=O.bodies.append(walls) sp=SpherePack() sp.makeCloud(mn,mx,num=100,rRelFuzz=0.333) sp.toSimulation(material='spheres') print('num bodies ', len(O.bodies)) triax=TriaxialStressController( maxMultiplier=1.+2e4/young, finalMaxMultiplier=1.+2e3/young, thickness = 0, stressMask = 7, internalCompaction=True, ) ThermalEngine = ThermalEngin
Re: [Yade-users] [Question #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Open => Needs information Robert Caulk requested more information: Hello, >>I compiled it from the source code again. The version is yade-2021-09-28.git-52f0c5a.But when I run the original script (unchanged),the same problem still exists... I cannot reproduce this issue using latest source. The example scripts run as intended. It may be an issue with distribution, you are running 20.04. Can you please install latest yade source on Ubuntu 18.04 and confirm your issue? >>I have tried the number of spheres as 100,200,300,500,1000,2000,and the O.dt as 1e-4 to 1e-6 (I also tried using the O.dt=PWaveTimeStep(), not sure if it's appropriate).The success of the simulation does not seem to have much to do with O.dt but O.dt affects the speed of the simulation(The premise is that the simulation is successful, that is, there is temperature instead of Nan) Have you actually looked at your packing? It seems you are simulating a gaseous state of spheres (although there is no way to know without an MWE). ThermalEngine was not designed to handle this kind of packing. Please provide an MWE. 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
Re: [Yade-users] [Question #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Ziyu Wang gave more information on the question: Hello Robert! I have the latest findings, which may help solve this problem: I follow your suggestion and consider changing O.dt and the number of spheres respectively.I found something interesting.. 1.I have tried the number of spheres as 100,200,300,500,1000,2000,and the O.dt as 1e-4 to 1e-6 (I also tried using the O.dt=PWaveTimeStep(), not sure if it's appropriate).The success of the simulation does not seem to have much to do with O.dt but O.dt affects the speed of the simulation(The premise is that the simulation is successful, that is, there is temperature instead of Nan) 2.As described in 1, I found that when the number of spheres is 100,200,1000,2000, the simulation is successful (i.e. the pore and body have their own temperature rather than Nan).But 300,500,the problem still exists,like the original script (unchanged).With these scripts,there is a same warning: CHOLMOD warning: matrix not positive definite. file: ../Supernodal/t_cholmod_super_numeric.c line: 911 (Maybe this is the reason of problem?) 3.As mentioned in 2, even in the successful cases, the curves are different.In the example script, I guess the ideal curve should be to keep the pore temperature at 343.15 and gradually increase the particle temperature from 333.15 to 343.15.If I'm right, this is true only when num = 100, and in other cases, the body temperature will increase to 353.15.. What is wrong with the problem? Best regards. -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Answered => Open Ziyu Wang is still having a problem: Thanks a lot Robert! >>Will you please install the latest source code version? You appear to have installed a snapshot from mid-March. Yes,I also noticed that.As you suggested[1], I used the following command: cd trunk git checkout -b working_state 61c90697 cd ../build cmake -DCMAKE_INSTALL_PREFIX=../install ../trunk make -j4 install And then the version changes to 2021-03-16. I compiled it from the source code again. The version is yade-2021-09-28.git-52f0c5a.But when I run the original script (unchanged),the same problem still exists... >>Sphere size and distribution changes the maximum allowable timestep for numerical stability. If you decrease sphere size, you will need to decrease O.dt to maintain stability. Ok,I will decrease O.dt and provide timely feedback. -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Robert Caulk proposed the following answer: Increasing sphere number in the same geometrical domain size, means smaller spheres. Smaller spheres means smaller characteristic lengths for the system. Smaller char. lengths means smaller allowable timestep. Please refer to this discussion for more details: https://answers.launchpad.net/yade/+question/698948 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
Re: [Yade-users] [Question #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Open => Answered Robert Caulk proposed the following answer: Hello, >>1.When I run the original script (unchanged), the problem is the same as before: plot is blank and body.State.Temp = Nan. Will you please install the latest source code version? You appear to have installed a snapshot from mid-March. >>So I guess the problem has something to do with the number of spheres? Sphere size and distribution changes the maximum allowable timestep for numerical stability. If you decrease sphere size, you will need to decrease O.dt to maintain stability. 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
Re: [Yade-users] [Question #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Description changed to: Dear all, As discussed in [1],I have installed the Source code and use the yade-2021-03-16.git-61c9069.I ran the script in example/thermalengine and found some new problems after my changes. The thermoHydroMechanical_coupling.py and flowScenario.py can run as expected.The problem is mainly in noFlowScenario.py. 1.When I run the original script (unchanged), the problem is the same as before: plot is blank and body.State.Temp = Nan. 2.My changes to the script are as follows: >>#sp = O.bodies.append(ymport.textExt('5cmEdge_1mm.spheres', >>'x_y_z_r',color=(0.1,0.1,0.9), material='spheres')) sp=pack.SpherePack() sp.makeCloud(mn,mx,rMean=0.0015,rRelFuzz=0.333,num=100,seed=1) sp.toSimulation(color=(0.752, 0.752, 0.752),material='spheres') (That is, change the way the sphere is generated,this is the only change) Then I ran the script again and was surprised to find that the simulation ran successfully as expected. However,I changed the num in makeCloud and found the bad thing:when the num=300,the problem mentioned above reappears:plot is blank and body.State.Temp = Nan... I also tried using randomDensePack,and the problem is the same. So I guess the problem has something to do with the number of spheres? The following is the version information in printAllVersions() Yade version : 2021-03-16.git-61c9069 Yade features : BoostLog RealHP mpmath PrecisionDouble Odeint VTK OpenMP GTS GUI-Qt5 CGAL PFVFLOW PFVFLOW LINSOLV MPI TWOPHASEFLOW FEMLIKE GL2PS LBMFLOW THERMAL PARTIALSAT PotentialParticles PotentialBlocks Yade config dir: ~/.yade-2021-03-16.git-61c9069 Yade precision : 53 bits, 15 decimal places, with mpmath Yade RealHP<…> : (15, 33, 45, 60, 120, 150, 300) decimal digits in C++, (15, 33) decimal digits accessible from python ``` Libraries used : | library | cmake| C++| | - | | -- | | boost | 107100 | 1.71.0 | | cgal | | 5.1.5 | | clp | 1.17.5 | 1.17.5 | | cmake | 3.16.3 || | coinutils | 2.11.4 | 2.11.4 | | compiler | /usr/bin/c++ 9.3.0 | gcc 9.3.0 | | eigen | 3.3.7| 3.3.7 | | freeglut | 2.8.1|| | gl| | 20190805 | | ipython | 7.13.0 || | metis | | 5.1.0 | | mpi | 3.1 | ompi:4.0.3 | | mpi4py| 3.0.3|| | mpmath| 1.1.0|| | openblas | | OpenBLAS 0.3.8| | python| 3.8.10 | 3.8.10 | | qglviewer | | 2.6.3 | | qt| | 5.12.8 | | sphinx| 1.8.5-final-0|| | sqlite| | 3.31.1 | | suitesparse | 5.7.1| 5.7.1 | | vtk | 6.3.0| 6.3.0 | ``` Linux version : Ubuntu 20.04.3 LTS Architecture : amd64 Little endian : True Thanks for your help! Best regards. [1]https://answers.launchpad.net/yade/+question/698830 -- 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