Re: [Yade-users] [Question #698948]: Stability of advection modeling by ThermalEngine
Question #698948 on Yade changed: https://answers.launchpad.net/yade/+question/698948 Zoheir Khademian posted a new comment: BTW, my Ubuntu version is: 20.04.3 -- 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 #698948]: Stability of advection modeling by ThermalEngine
Question #698948 on Yade changed: https://answers.launchpad.net/yade/+question/698948 Zoheir Khademian posted a new comment: Hi Robert, >>This is quite odd: applyTempDeltaToSolids is deactivated by default: the user needs to manually set thermal.delT>0 for that function to activate I think the problem is delT is being used as a variable here as well [1], so during fluid-fluid conduction calculations, delT becomes positive and it activates applyTempDeltaToSolids automatically. I deleted detT variable here [1] and directly used "cell->info().temp() - neighborCell->info().temp()" for conductionEnergy calculation here[2]. Then recompiled and the results are stable. Let me know if this does not makes sense to you. [1] https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/pfv/Thermal.cpp#L530 [2] https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/pfv/Thermal.cpp#L545 In any case, here is the list of particles (all the same size) and the MWE to run the model. Thanks - Zoheir #format x_y_z_r 0.0177741 0.0293044 0.0178294 0.0015 0.0322250.0197924 0.0178472 0.0015 0.0322266 0.0291739 0.0178539 0.0015 0.0279236 0.0322490.0261029 0.0015 0.0322118 0.0302264 0.0320853 0.0015 0.0272061 0.0316643 0.0320534 0.0015 0.0293610.0178245 0.0290514 0.0015 0.0223492 0.0237293 0.0212975 0.0015 0.0322662 0.0177355 0.0290630.0015 0.0249392 0.0297701 0.0178535 0.0015 0.0278041 0.0253036 0.0320824 0.0015 0.0209031 0.0239399 0.0259082 0.0015 0.0299960.0247971 0.0227321 0.0015 0.0231086 0.0224061 0.0178694 0.0015 0.0237889 0.0267710.0280437 0.0015 0.0271796 0.0296372 0.0298686 0.0015 0.0293275 0.0207069 0.0294883 0.0015 0.0297280.0307252 0.0238577 0.0015 0.0294550.0177438 0.03204 0.0015 0.0220291 0.0251347 0.0320876 0.0015 0.0228457 0.0177657 0.0205171 0.0015 0.0247976 0.0201580.0187476 0.0015 0.0221528 0.0177930.0297048 0.0015 0.0177809 0.0178246 0.0320527 0.0015 0.0178286 0.0204492 0.0192357 0.0015 0.0191019 0.0272907 0.0227524 0.0015 0.0231820.0266745 0.0241058 0.0015 0.0228497 0.0296931 0.0320761 0.0015 0.0203349 0.0184782 0.0277075 0.0015 0.0178115 0.0232242 0.0181458 0.0015 0.0236951 0.0187262 0.0273578 0.0015 0.0178150.0178142 0.0290756 0.0015 0.0268862 0.0177484 0.0275311 0.0015 0.0293104 0.0293306 0.0320274 0.0015 0.0178224 0.0228895 0.0229322 0.0015 0.0229362 0.0182461 0.0233847 0.0015 0.0234661 0.0204362 0.0296891 0.0015 0.0308867 0.0228868 0.0282905 0.0015 0.0322099 0.0250677 0.0296246 0.0015 0.0322810.0178211 0.0320448 0.0015 0.0322146 0.0204259 0.0245334 0.0015 0.0322253 0.0322509 0.0239898 0.0015 0.0279165 0.0228391 0.0279433 0.0015 0.0234598 0.0322574 0.0219894 0.0015 0.0202221 0.0215691 0.0178190.0015 0.0292976 0.0252520.0289125 0.0015 0.0178179 0.0322929 0.0207288 0.0015 0.0219679 0.0204761 0.0198587 0.0015 0.0219003 0.0322297 0.0244164 0.0015 0.0188566 0.0202206 0.0221634 0.0015 0.0276402 0.0177882 0.0232129 0.0015 0.0178161 0.0322694 0.0269258 0.0015 0.0295711 0.0280749 0.0290676 0.0015 0.0265758 0.0178165 0.0320975 0.0015 0.0321235 0.0295753 0.0253156 0.0015 0.0203762 0.0177426 0.0220801 0.0015 0.0292814 0.0231458 0.0251311 0.0015 0.0299286 0.0268797 0.0186235 0.0015 0.0288255 0.0202145 0.0243078 0.0015 0.0178049 0.0276979 0.0202425 0.0015 0.0322722 0.0256006 0.0178161 0.0015 0.0223008 0.0202683 0.0253642 0.0015 0.0255165 0.0266773 0.0206327 0.0015 0.0199502 0.0216179 0.0289218 0.0015 0.0301260.0177624 0.0178698 0.0015 0.0254739 0.0197561 0.0231361 0.0015 0.0322638 0.0228734 0.0228553 0.0015 0.0265408 0.0222830.0178242 0.0015 0.0254871 0.0259749 0.0258332 0.0015 0.0322915 0.0227056 0.0178061 0.0015 0.0319731 0.0202977 0.0214805 0.0015 0.0235946 0.0323112 0.0178428 0.0015 0.0250886 0.0203689
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 #699071]: About the yadedaily update
Question #699071 on Yade changed: https://answers.launchpad.net/yade/+question/699071 Status: Answered => Solved Ziyu Wang confirmed that the question is solved: Hello everyone, Good news, I can also update to 9b2e38e through package manager now. Thank you for your efforts. My problem has been solved. -- 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 #699135]: How to reduce average particle velocity of a packing
Question #699135 on Yade changed: https://answers.launchpad.net/yade/+question/699135 Status: Answered => Solved Leonard confirmed that the question is solved: Thanks Jan Stránský, 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 #699157]: test a new contact model
Question #699157 on Yade changed: https://answers.launchpad.net/yade/+question/699157 Status: Answered => Solved Chenguang Liang confirmed that the question is solved: thanks,my friend. I will try to test the simulations as you said. When there is no bugs and other problems, i am happy to share it on yade to make it better. -- 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 #699157]: test a new contact model
Question #699157 on Yade changed: https://answers.launchpad.net/yade/+question/699157 Status: Open => Answered Jan Stránský proposed the following answer: Hello, TL;TR: it is very individual > Are there any other tests necessary to check my code ... ? It depends on many factors and on the definition of "necessary". How many people / how often will use the model? Just you for one article? You and your team for a longer-term project? You and your team and you want to publish it and let anybody use it? ...? The more people you want to use it, the better it should be tested (IMO). The answer strongly depends on the model, its purpose and its typical usage (please provide more information). If its primary purpose is some behavior at bouncing and rebouncing (which according to JKR and wikipedia seems yes), then probably a bouncing simulation is OK. Concerning bounding simulation, it strongly depends on what "following the tutorial file" means (please provide more information). You can test one specific case. On the other hand, even testing "just" bouncing, you can test different radii of spheres, different angles of impact, different impact velocities, impact + rotating sphere etc. etc. etc. All in all, it is very individual, depends on many factors etc. Cheers Jan -- 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 #699030]: How to set the particle size when generating particles
Question #699030 on Yade changed: https://answers.launchpad.net/yade/+question/699030 Status: Open => Answered Jan Stránský proposed the following answer: Please provide a MWE [1], W = working, i.e. we can copy-paste it and test it. Here e.g. N_particles is not defined, 'sp' is not defined etc.. M part is perfect :-) Have a look at example script [2] and let us know if it helps or if you need further help. Cheers Jan [1] https://www.yade-dem.org/wiki/Howtoask [2] https://gitlab.com/yade-dev/trunk/-/blob/master/examples/test/psd.py -- 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 #698948]: Stability of advection modeling by ThermalEngine
Question #698948 on Yade changed: https://answers.launchpad.net/yade/+question/698948 Robert Caulk requested more information: Also, please indicate your version of Ubuntu. -- 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 #692131]: acoustic emission(JCFmat)recordCrack=true
Question #692131 on Yade changed: https://answers.launchpad.net/yade/+question/692131 d sd posted a new comment: To see if the hassle is coming from the getBoundaryFlux() feature, you may take a look at the waft speed (or the flux) at different factors positioned inside the fracture. If the rate is the same everywhere, you can even use it to compute the permeability of your fracture with out the usage of the Qin and Qout. Learn more about it here https://bright- star.com.au/product-category/firecrackers/ -- 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 #698948]: Stability of advection modeling by ThermalEngine
Question #698948 on Yade changed: https://answers.launchpad.net/yade/+question/698948 Robert Caulk proposed the following answer: Hey Zoheir, This is quite odd: applyTempDeltaToSolids is deactivated by default: the user needs to manually set thermal.delT>0 for that function to activate [1]. As you see in [1], thermal.delT is set to 0 if user does nothing. So compiling without the line you indicated will change nothing in the output unless your user script manually activated thermal.delT>0. I don’t see you setting that value in any of the MWEs that you provided here. So this cannot be the issue (unless you are using a different script to the one provided here where you use delT). Can you please find the unstable packing, then export those sphere radii and location to text file, the copy paste that file to this thread so I can replicate exactly your problem. Also please copy paste the identical script that goes with running that packing. Then please type here the exact output of yade -v where yade is the executable you are using (it may be named differently depending on how you compile or which package you download). Cheers, Robert [1] https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/pfv/Thermal.hpp#L133 Le mar. 19 oct. 2021 à 23:01, Zoheir Khademian < question698...@answers.launchpad.net> a écrit : > Question #698948 on Yade changed: > https://answers.launchpad.net/yade/+question/698948 > > Zoheir Khademian posted a new comment: > Hey Robert, > Sorry for my delays. I try to give your comments enough thoughts and do > investigation before replying. > > I tried both released version [1] and development version [2] and still > observed the instabilities meaning the particle temp goes above the max > boundary temp while cell temp is still within the range. > > > Regarding the diffusivity coefficient kA/L, in the MWE I uploaded last, I > have reduced kA/L by reducing fluidk and fluidConductionAreaFactor to > arbitrarily low values. I also tried your values > (thermal.minimumFluidCondDist=rMean and > thermal.fluidConductionAreaFactor=0.1) but still the particle temp is > unstable. > > Finally I looked at the thermal.cpp and noticed that > applyTempDeltaToSolids [3] adds pore temp differences to the solid temp > for no apparent reason (particle temp is updated here [4]). I > deactivated applyTempDeltaToSolids and recompiled the source. I ran the > MWE and the results are stable now. > > Let me know if applyTempDeltaToSolids must not be deactivated or I am > missing something here > > Thanks > Zoheir > > > [1] https://launchpad.net/yade/+download > > [2] https://gitlab.com/yade-dev/trunk > > [3] https://gitlab.com/yade- > dev/trunk/-/blob/master/pkg/pfv/Thermal.cpp#L94 > > [4] https://gitlab.com/yade- > dev/trunk/-/blob/master/pkg/pfv/Thermal.cpp#L605 > > -- > 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 > -- 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
