Re: [Yade-users] [Question #695065]: ThermoHydro analyses
Question #695065 on Yade changed: https://answers.launchpad.net/yade/+question/695065 Zoheir Khademian posted a new comment: > You are describing a lot of things that are not represented in your MWE. So perhaps your implentation is wrong. No way to know without an MWE: it would just be pure speculation for me to comment on this. This is a fair point. I will make a more representative MWE and post it here. 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 #695065]: ThermoHydro analyses
Question #695065 on Yade changed: https://answers.launchpad.net/yade/+question/695065 Robert Caulk proposed the following answer: >The documentation says the minimumPorosity [1] is also included in the flowEngine [2] Right, TwoPhaseFlowEngine is not the same as FlowEngine. So you should be using the FlowEngine class (based on the MWE you provided). >If I dont remove the bonds between my grains but still compact them in a cube, the packing is normal but the TH analysis gets unstable. This may be because the boundaries are not as smooth as the first case, owing to the irregularity in the shape of my bonded grains. Solving this case is what I am looking for. Do you think alpha boundary can solve the issue in the second case? You are describing a lot of things that are not represented in your MWE. So perhaps your implentation is wrong. No way to know without an MWE: it would just be pure speculation for me to comment on this. -- 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 #695065]: ThermoHydro analyses
Question #695065 on Yade changed: https://answers.launchpad.net/yade/+question/695065 Zoheir Khademian posted a new comment: Hi Robert, > TwoPhaseFlowEngine.minimumPorosity, that you are referencing, is part of TwoPhaseFlowEngine, not ThermalEngine. So it will not change your results. The documentation says the minimumPorosity [1] is also included in the flowEngine [2] and apparently it makes the model stable at the beginning although after a while it gets unstable again. - Thanks for the offer for including alpha boundary in the ThermalEngine. I will get back to you on that. I should also clarify two points: - When I remove the bonds between my grains (clusters of spheres) and compact my sample in a cube, TH (thermo-hydro) analysis works just fine. I guess this is because the packing is normal (cubic) in this case and boundaries are pretty smooth due to the fact that clusters are dissembled and I am left with a packing of poly-disperse spheres instead of a packing of irregular-shaped, bonded grains. - If I dont remove the bonds between my grains but still compact them in a cube, the packing is normal but the TH analysis gets unstable. This may be because the boundaries are not as smooth as the first case, owing to the irregularity in the shape of my bonded grains. Solving this case is what I am looking for. Do you think alpha boundary can solve the issue in the second case? Thanks for the help, [1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngine.minimumPorosity [2] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngine -- 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 #695065]: ThermoHydro analyses
Question #695065 on Yade changed: https://answers.launchpad.net/yade/+question/695065 Status: Open => Answered Robert Caulk proposed the following answer: Hello, TwoPhaseFlowEngine.minimumPorosity, that you are referencing, is part of TwoPhaseFlowEngine, not ThermalEngine. So it will not change your results. >Can you think of any other way around except for making boundaries normal? Really depends on what it is you are simulating. I guess you could use a normal packing and then go in and control the cell/particle conditions based on coordinates. In the end you would hopefully have your irregular shaped particle embedded within a cube. It *could* work but still tough to say if it will actually be useful without knowing what your end goal is for this simulation. Otherwise, alpha boundary conditions are designed to adapt to any shape. So presumably you could just use your grain as is. But alpha boundary conditions are not currently available in ThermalEngine so you would need to go through and adapt the ThermalEngine components in the source code. The modifications are moderate, but surgical. This, it is possible, but this will require quite intimate knowledge of the FlowEngine, the triangulation behind the FlowEngine (i.e. CGAL), and ThermalEngine. So I would say it is certainly not easy. If you don't know C++ or Yade source, it might be insurmountable in less than a month or two. I guess this would take me about 8-10 hours to do properly and fully debug all the affected Thermal components. I do not have the time to do that for free, but you can hire me if you need it done. 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 #695065]: ThermoHydro analyses
Question #695065 on Yade changed: https://answers.launchpad.net/yade/+question/695065 Status: Needs information => Open Zoheir Khademian gave more information on the question: Hey Robert, Sorry for deviating from the forum rules again. Will make sure to follow them in the future. I tried a normal shape packing and all works just fine but for the problem I am studying I have irregular-shape grains and the boundaries cannot be really smooth or normal. I tired limiting the minimum porosity to 0.01 by [1]. By doing this, temperatures start updating at the beginning very nicely but after a few timestep, it gets unstable and temps keep rising indefinitely. I tried reducing timesteps but the same issue repeated. Can you think of any other way around except for making boundaries normal? Thanks, [1] https://yade- dem.org/doc/yade.wrapper.html#yade.wrapper.TwoPhaseFlowEngine.minimumPorosity -- 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 #695065]: ThermoHydro analyses
Question #695065 on Yade changed: https://answers.launchpad.net/yade/+question/695065 Status: Open => Needs information Robert Caulk requested more information: I ran your script and notice that the packing you are using does not have a normal shape (not a cuboid, cylinder, or sphere). With such a random shape, I would expect undefined behavior from the definition of the boundaries which would result in an ill posed problem. Can you make your packing into a normal shape? -- 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 #695065]: ThermoHydro analyses
Question #695065 on Yade changed: https://answers.launchpad.net/yade/+question/695065 Robert Caulk posted a new comment: Hey Zoheir, Please do not change history. It renders the thread confusing and difficult to follow. Better to post the requested material as a new comment instead of revising the question. It may seem as though we have too many rules around here, but they are to keep things organized. This is a knowledge base that we are spending our time to build. It is a waste of our time to build disorganized threads that no one can read. Thanks for considering the Yade community, 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 #695065]: ThermoHydro analyses
Question #695065 on Yade changed:
https://answers.launchpad.net/yade/+question/695065
Description changed to:
Hello,
I am trying to run a TH analysis with warm fluid going through a cold
pack of polydisperse particles using the provided example in [1].
When I use a previously compacted, poly-disperse-particle sample for
this analysis, the temperature of all particles remains constant at
first and then turns to "nan". I was wondering what needs to be adjusted
in the thermal and flow engines following the features of the granular
sample under study.
[1]
https://gitlab.com/yade-dev/trunk/-/blob/master/examples/ThermalEngine/flowScenario.py
The MWE follows:
Thanks - Zoheir
from yade import pack, ymport
from yade import timing
import numpy as np
import shutil
young=1e9
rad=0.003
mn,mx=Vector3(0.0005944,0.0005596,0.00099759),Vector3(0.00119,0.00119,0.0016) #
corners of the initial packing
thermalCond = 2. #W/(mK)
heatCap = 710. #J(kg K)
t0 = 333.15 #K
# micro properties
r = rad
k = 2.0 # 2*k*r
Cp = 710.
rho = 2600.
D = 2.*r
m = 4./3.*np.pi*r**2/rho
# macro diffusivity
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 = O.bodies.append(ymport.textExt('ShCube', 'x_y_z_r',color=(0.1,0.1,0.9),
material='spheres'))
O.bodies.append([
sphere((0.000945595,0.000784186,0.00149986),0.875846,material='spheres'),
sphere((0.000778022,0.000753066,0.00145943),0.875846,material='spheres'),
sphere((0.000888722,0.000810072,0.00133622),0.875846,material='spheres'),
sphere((0.000903823,0.000647224,0.00139896),0.875846,material='spheres'),
sphere((0.00105629,0.000841192,0.00137666),0.875846,material='spheres'),
sphere((0.00101452,0.00070423,0.00127576),0.875846,material='spheres'),
sphere((0.00084695,0.00067311,0.00123532),0.875846,material='spheres'),
sphere((0.000999422,0.000867078,0.00121301),0.875846,material='spheres'),
sphere((0.00103358,0.00108745,0.0011483),0.815788,material='spheres'),
sphere((0.000902032,0.00103675,0.00153707),0.735711,material='spheres'),
sphere((0.000986413,0.00091621,0.00153692),0.735711,material='spheres'),
sphere((0.00111325,0.000915283,0.00146233),0.735711,material='spheres'),
sphere((0.00108202,0.000666186,0.00154242),0.665643,material='spheres'),
sphere((0.00111005,0.000740653,0.00132599),0.665643,material='spheres'),
sphere((0.00080908,0.00109404,0.0014756),0.499888,material='spheres'),
sphere((0.000751015,0.00103522,0.00153186),0.499888,material='spheres'),
sphere((0.000670623,0.00108745,0.00150348),0.499888,material='spheres'),
sphere((0.000687866,0.00108202,0.00126569),0.459731,material='spheres'),
sphere((0.000703081,0.00107221,0.00135584),0.459731,material='spheres'),
sphere((0.00110443,0.00083307,0.00156066),0.420579,material='spheres'),
sphere((0.00112293,0.000765198,0.00151454),0.420579,material='spheres'),
sphere((0.00106475,0.000819695,0.0014877),0.420579,material='spheres'),
sphere((0.0011389,0.000814544,0.00144832),0.420579,material='spheres'),
sphere((0.00112143,0.000746581,0.00109922),0.439814,material='spheres'),
sphere((0.00105308,0.000800634,0.0023),0.439814,material='spheres'),
sphere((0.00110668,0.000870009,0.00110401),0.439814,material='spheres'),
sphere((0.00112389,0.000807084,0.00116302),0.439814,material='spheres'),
sphere((0.00103833,0.000924061,0.00111603),0.439814,material='spheres'),
sphere((0.000792443,0.000977976,0.00151641),0.220555,material='spheres'),
sphere((0.000760407,0.000979304,0.00148612),0.220555,material='spheres'),
sphere((0.000836178,0.000937735,0.0014684),0.220555,material='spheres'),
sphere((0.000775855,0.00092188,0.00146726),0.220555,material='spheres'),
sphere((0.000798293,0.000958519,0.00147726),0.220555,material='spheres'),
sphere((0.00073797,0.000942664,0.00147612),0.220555,material='spheres'),
sphere((0.000766257,0.000959847,0.00144696),0.220555,material='spheres'),
sphere((0.000728372,0.000980632,0.00145582),0.220555,material='spheres'),
sphere((0.00081374,0.000901096,0.0014584),0.220555,material='spheres'),
sphere((0.000842027,0.000918279,0.00142924),0.220555,material='spheres'),
sphere((0.000781704,0.000902424,0.00142811),0.220555,material='spheres'),
sphere((0.000804142,0.000939063,0.0014381),0.220555,material='spheres'),
sphere((0.000772106,0.000940391,0.00140781),0.220555,material='spheres'),
sphere((0.000734221,0.000961175,0.00141667),0.220555,material='spheres'),
sphere((0.000809991,0.000919607,0.00139895),0.220555,material='spheres'),
sphere((0.00113905,0.00100891,0.00119864),0.235576,material='spheres'),
sphere((0.000698534,0.000815074,0.00126433),0.765702,material='spheres'),
Re: [Yade-users] [Question #695065]: ThermoHydro analyses
Question #695065 on Yade changed:
https://answers.launchpad.net/yade/+question/695065
Zoheir Khademian gave more information on the question:
Sorry for the confusion.
The script [1] works perfectly but when I replace the packing with my
own packing, the body temperatures do not get calculated ("nan").
Comparing to [1], I changed the wall dimensions to fit the sample and
removed post-processing fictions ("pressureField" and "endFlux") .
Thanks for the help
Zoheir
[1] https://gitlab.com/yade-
dev/trunk/-/blob/master/examples/ThermalEngine/flowScenario.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 #695065]: ThermoHydro analyses
Question #695065 on Yade changed:
https://answers.launchpad.net/yade/+question/695065
Description changed to:
Hello,
I am trying to run a TH analysis with warm fluid going through a cold
pack of polydisperse particles using the provided example in [1].
When I use a previously compacted, poly-disperse-particle sample for
this analysis, the temperature of all particles remains constant at
first and then turns to "nan". I was wondering what needs to be adjusted
in the thermal and flow engines following the features of the granular
sample under study.
[1]
https://gitlab.com/yade-dev/trunk/-/blob/master/examples/ThermalEngine/flowScenario.py
The MWE follows:
Thanks - Zoheir
from yade import pack, ymport
from yade import timing
import numpy as np
import shutil
young=1e9
rad=0.003
mn,mx=Vector3(0.0005944,0.0005596,0.00099759),Vector3(0.00119,0.00119,0.0016) #
corners of the initial packing
thermalCond = 2. #W/(mK)
heatCap = 710. #J(kg K)
t0 = 333.15 #K
# micro properties
r = rad
k = 2.0 # 2*k*r
Cp = 710.
rho = 2600.
D = 2.*r
m = 4./3.*np.pi*r**2/rho
# macro diffusivity
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 = O.bodies.append(ymport.textExt('ShCube', 'x_y_z_r',color=(0.1,0.1,0.9),
material='spheres'))
O.bodies.append([
sphere((0.000945595,0.000784186,0.00149986),0.875846),
sphere((0.000778022,0.000753066,0.00145943),0.875846),
sphere((0.000888722,0.000810072,0.00133622),0.875846),
sphere((0.000903823,0.000647224,0.00139896),0.875846),
sphere((0.00105629,0.000841192,0.00137666),0.875846),
sphere((0.00101452,0.00070423,0.00127576),0.875846),
sphere((0.00084695,0.00067311,0.00123532),0.875846),
sphere((0.000999422,0.000867078,0.00121301),0.875846),
sphere((0.00103358,0.00108745,0.0011483),0.815788),
sphere((0.000902032,0.00103675,0.00153707),0.735711),
sphere((0.000986413,0.00091621,0.00153692),0.735711),
sphere((0.00111325,0.000915283,0.00146233),0.735711),
sphere((0.00108202,0.000666186,0.00154242),0.665643),
sphere((0.00111005,0.000740653,0.00132599),0.665643),
sphere((0.00080908,0.00109404,0.0014756),0.499888),
sphere((0.000751015,0.00103522,0.00153186),0.499888),
sphere((0.000670623,0.00108745,0.00150348),0.499888),
sphere((0.000687866,0.00108202,0.00126569),0.459731),
sphere((0.000703081,0.00107221,0.00135584),0.459731),
sphere((0.00110443,0.00083307,0.00156066),0.420579),
sphere((0.00112293,0.000765198,0.00151454),0.420579),
sphere((0.00106475,0.000819695,0.0014877),0.420579),
sphere((0.0011389,0.000814544,0.00144832),0.420579),
sphere((0.00112143,0.000746581,0.00109922),0.439814),
sphere((0.00105308,0.000800634,0.0023),0.439814),
sphere((0.00110668,0.000870009,0.00110401),0.439814),
sphere((0.00112389,0.000807084,0.00116302),0.439814),
sphere((0.00103833,0.000924061,0.00111603),0.439814),
sphere((0.000792443,0.000977976,0.00151641),0.220555),
sphere((0.000760407,0.000979304,0.00148612),0.220555),
sphere((0.000836178,0.000937735,0.0014684),0.220555),
sphere((0.000775855,0.00092188,0.00146726),0.220555),
sphere((0.000798293,0.000958519,0.00147726),0.220555),
sphere((0.00073797,0.000942664,0.00147612),0.220555),
sphere((0.000766257,0.000959847,0.00144696),0.220555),
sphere((0.000728372,0.000980632,0.00145582),0.220555),
sphere((0.00081374,0.000901096,0.0014584),0.220555),
sphere((0.000842027,0.000918279,0.00142924),0.220555),
sphere((0.000781704,0.000902424,0.00142811),0.220555),
sphere((0.000804142,0.000939063,0.0014381),0.220555),
sphere((0.000772106,0.000940391,0.00140781),0.220555),
sphere((0.000734221,0.000961175,0.00141667),0.220555),
sphere((0.000809991,0.000919607,0.00139895),0.220555),
sphere((0.00113905,0.00100891,0.00119864),0.235576),
sphere((0.000698534,0.000815074,0.00126433),0.765702),
sphere((0.000695882,0.00071517,0.00138037),0.765702),
sphere((0.000670978,0.000701228,0.00116568),0.765702),
sphere((0.000683566,0.000664341,0.00153398),0.481936),
sphere((0.00113412,0.000910585,0.00138428),0.40366),
sphere((0.00112752,0.00102435,0.0013913),0.40366),
sphere((0.0011135,0.000914201,0.0013063),0.40366),
sphere((0.00100741,0.000944425,0.00133576),0.40366),
sphere((0.00108155,0.000966389,0.00135897),0.40366),
sphere((0.0011069,0.00102797,0.00131333),0.40366),
sphere((0.00102898,0.00102219,0.00133366),0.40366),
sphere((0.00107495,0.00108015,0.00136599),0.40366),
sphere((0.00112847,0.00110573,0.00131123),0.40366),
sphere((0.00113885,0.000975778,0.00126066),0.40366),
sphere((0.00106093,0.000970005,0.001281),0.40366),
sphere((0.00105433,0.00108377,0.00128802),0.40366),
Re: [Yade-users] [Question #695065]: ThermoHydro analyses
Question #695065 on Yade changed:
https://answers.launchpad.net/yade/+question/695065
Description changed to:
Hello,
I am trying to run a TH analysis with warm fluid going through a cold
pack of polydisperse particles using the provided example in [1].
When I use a previously compacted, poly-disperse-particle sample for
this analysis, the temperature of all particles remains constant at
first and then turns to "nan". I was wondering what needs to be adjusted
in the thermal and flow engines following the features of the granular
sample under study.
[1]
https://gitlab.com/yade-dev/trunk/-/commit/f0c96c63937a666d1fa1529173b554ac74b04456
The MWE follows:
Thanks - Zoheir
from yade import pack, ymport
from yade import timing
import numpy as np
import shutil
young=1e9
rad=0.003
mn,mx=Vector3(0.0005944,0.0005596,0.00099759),Vector3(0.00119,0.00119,0.0016) #
corners of the initial packing
thermalCond = 2. #W/(mK)
heatCap = 710. #J(kg K)
t0 = 333.15 #K
# micro properties
r = rad
k = 2.0 # 2*k*r
Cp = 710.
rho = 2600.
D = 2.*r
m = 4./3.*np.pi*r**2/rho
# macro diffusivity
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 = O.bodies.append(ymport.textExt('ShCube', 'x_y_z_r',color=(0.1,0.1,0.9),
material='spheres'))
O.bodies.append([
sphere((0.000945595,0.000784186,0.00149986),0.875846),
sphere((0.000778022,0.000753066,0.00145943),0.875846),
sphere((0.000888722,0.000810072,0.00133622),0.875846),
sphere((0.000903823,0.000647224,0.00139896),0.875846),
sphere((0.00105629,0.000841192,0.00137666),0.875846),
sphere((0.00101452,0.00070423,0.00127576),0.875846),
sphere((0.00084695,0.00067311,0.00123532),0.875846),
sphere((0.000999422,0.000867078,0.00121301),0.875846),
sphere((0.00103358,0.00108745,0.0011483),0.815788),
sphere((0.000902032,0.00103675,0.00153707),0.735711),
sphere((0.000986413,0.00091621,0.00153692),0.735711),
sphere((0.00111325,0.000915283,0.00146233),0.735711),
sphere((0.00108202,0.000666186,0.00154242),0.665643),
sphere((0.00111005,0.000740653,0.00132599),0.665643),
sphere((0.00080908,0.00109404,0.0014756),0.499888),
sphere((0.000751015,0.00103522,0.00153186),0.499888),
sphere((0.000670623,0.00108745,0.00150348),0.499888),
sphere((0.000687866,0.00108202,0.00126569),0.459731),
sphere((0.000703081,0.00107221,0.00135584),0.459731),
sphere((0.00110443,0.00083307,0.00156066),0.420579),
sphere((0.00112293,0.000765198,0.00151454),0.420579),
sphere((0.00106475,0.000819695,0.0014877),0.420579),
sphere((0.0011389,0.000814544,0.00144832),0.420579),
sphere((0.00112143,0.000746581,0.00109922),0.439814),
sphere((0.00105308,0.000800634,0.0023),0.439814),
sphere((0.00110668,0.000870009,0.00110401),0.439814),
sphere((0.00112389,0.000807084,0.00116302),0.439814),
sphere((0.00103833,0.000924061,0.00111603),0.439814),
sphere((0.000792443,0.000977976,0.00151641),0.220555),
sphere((0.000760407,0.000979304,0.00148612),0.220555),
sphere((0.000836178,0.000937735,0.0014684),0.220555),
sphere((0.000775855,0.00092188,0.00146726),0.220555),
sphere((0.000798293,0.000958519,0.00147726),0.220555),
sphere((0.00073797,0.000942664,0.00147612),0.220555),
sphere((0.000766257,0.000959847,0.00144696),0.220555),
sphere((0.000728372,0.000980632,0.00145582),0.220555),
sphere((0.00081374,0.000901096,0.0014584),0.220555),
sphere((0.000842027,0.000918279,0.00142924),0.220555),
sphere((0.000781704,0.000902424,0.00142811),0.220555),
sphere((0.000804142,0.000939063,0.0014381),0.220555),
sphere((0.000772106,0.000940391,0.00140781),0.220555),
sphere((0.000734221,0.000961175,0.00141667),0.220555),
sphere((0.000809991,0.000919607,0.00139895),0.220555),
sphere((0.00113905,0.00100891,0.00119864),0.235576),
sphere((0.000698534,0.000815074,0.00126433),0.765702),
sphere((0.000695882,0.00071517,0.00138037),0.765702),
sphere((0.000670978,0.000701228,0.00116568),0.765702),
sphere((0.000683566,0.000664341,0.00153398),0.481936),
sphere((0.00113412,0.000910585,0.00138428),0.40366),
sphere((0.00112752,0.00102435,0.0013913),0.40366),
sphere((0.0011135,0.000914201,0.0013063),0.40366),
sphere((0.00100741,0.000944425,0.00133576),0.40366),
sphere((0.00108155,0.000966389,0.00135897),0.40366),
sphere((0.0011069,0.00102797,0.00131333),0.40366),
sphere((0.00102898,0.00102219,0.00133366),0.40366),
sphere((0.00107495,0.00108015,0.00136599),0.40366),
sphere((0.00112847,0.00110573,0.00131123),0.40366),
sphere((0.00113885,0.000975778,0.00126066),0.40366),
sphere((0.00106093,0.000970005,0.001281),0.40366),
sphere((0.00105433,0.00108377,0.00128802),0.40366),
Re: [Yade-users] [Question #695065]: ThermoHydro analyses
Question #695065 on Yade changed: https://answers.launchpad.net/yade/+question/695065 Status: Needs information => Open Zoheir Khademian gave more information on the question: Thanks, Robert. I revised the question following the forum 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 #695065]: ThermoHydro analyses
Question #695065 on Yade changed:
https://answers.launchpad.net/yade/+question/695065
Description changed to:
Hello,
I am trying to run a TH analysis with warm fluid going through a cold
pack of polydisperse particles using the provided example in [1].
When I use a previously compacted, poly-disperse-particle sample for
this analysis, the temperature of all particles remains constant at
first and then turns to "nan". I was wondering what needs to be adjusted
in the thermal and flow engines following the features of the granular
sample under study.
[1]
https://gitlab.com/yade-dev/trunk/-/commit/f0c96c63937a666d1fa1529173b554ac74b04456
The MWE follows:
Thanks - Zoheir
from yade import pack, ymport
from yade import timing
import numpy as np
import shutil
young=1e9
rad=0.003
mn,mx=Vector3(0.0005944,0.0005596,0.00099759),Vector3(0.00119,0.00119,0.0016) #
corners of the initial packing
thermalCond = 2. #W/(mK)
heatCap = 710. #J(kg K)
t0 = 333.15 #K
# micro properties
r = rad
k = 2.0 # 2*k*r
Cp = 710.
rho = 2600.
D = 2.*r
m = 4./3.*np.pi*r**2/rho
# macro diffusivity
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 = O.bodies.append(ymport.textExt('ShCube', 'x_y_z_r',color=(0.1,0.1,0.9),
material='spheres'))
O.bodies.append([
sphere((0.000945595,0.000784186,0.00149986),0.875846),
sphere((0.000778022,0.000753066,0.00145943),0.875846),
sphere((0.000888722,0.000810072,0.00133622),0.875846),
sphere((0.000903823,0.000647224,0.00139896),0.875846),
sphere((0.00105629,0.000841192,0.00137666),0.875846),
sphere((0.00101452,0.00070423,0.00127576),0.875846),
sphere((0.00084695,0.00067311,0.00123532),0.875846),
sphere((0.000999422,0.000867078,0.00121301),0.875846),
sphere((0.00103358,0.00108745,0.0011483),0.815788),
sphere((0.000902032,0.00103675,0.00153707),0.735711),
sphere((0.000986413,0.00091621,0.00153692),0.735711),
sphere((0.00111325,0.000915283,0.00146233),0.735711),
sphere((0.00108202,0.000666186,0.00154242),0.665643),
sphere((0.00111005,0.000740653,0.00132599),0.665643),
sphere((0.00080908,0.00109404,0.0014756),0.499888),
sphere((0.000751015,0.00103522,0.00153186),0.499888),
sphere((0.000670623,0.00108745,0.00150348),0.499888),
sphere((0.000687866,0.00108202,0.00126569),0.459731),
sphere((0.000703081,0.00107221,0.00135584),0.459731),
sphere((0.00110443,0.00083307,0.00156066),0.420579),
sphere((0.00112293,0.000765198,0.00151454),0.420579),
sphere((0.00106475,0.000819695,0.0014877),0.420579),
sphere((0.0011389,0.000814544,0.00144832),0.420579),
sphere((0.00112143,0.000746581,0.00109922),0.439814),
sphere((0.00105308,0.000800634,0.0023),0.439814),
sphere((0.00110668,0.000870009,0.00110401),0.439814),
sphere((0.00112389,0.000807084,0.00116302),0.439814),
sphere((0.00103833,0.000924061,0.00111603),0.439814),
sphere((0.000792443,0.000977976,0.00151641),0.220555),
sphere((0.000760407,0.000979304,0.00148612),0.220555),
sphere((0.000836178,0.000937735,0.0014684),0.220555),
sphere((0.000775855,0.00092188,0.00146726),0.220555),
sphere((0.000798293,0.000958519,0.00147726),0.220555),
sphere((0.00073797,0.000942664,0.00147612),0.220555),
sphere((0.000766257,0.000959847,0.00144696),0.220555),
sphere((0.000728372,0.000980632,0.00145582),0.220555),
sphere((0.00081374,0.000901096,0.0014584),0.220555),
sphere((0.000842027,0.000918279,0.00142924),0.220555),
sphere((0.000781704,0.000902424,0.00142811),0.220555),
sphere((0.000804142,0.000939063,0.0014381),0.220555),
sphere((0.000772106,0.000940391,0.00140781),0.220555),
sphere((0.000734221,0.000961175,0.00141667),0.220555),
sphere((0.000809991,0.000919607,0.00139895),0.220555),
sphere((0.00113905,0.00100891,0.00119864),0.235576),
sphere((0.000698534,0.000815074,0.00126433),0.765702),
sphere((0.000695882,0.00071517,0.00138037),0.765702),
sphere((0.000670978,0.000701228,0.00116568),0.765702),
sphere((0.000683566,0.000664341,0.00153398),0.481936),
sphere((0.00113412,0.000910585,0.00138428),0.40366),
sphere((0.00112752,0.00102435,0.0013913),0.40366),
sphere((0.0011135,0.000914201,0.0013063),0.40366),
sphere((0.00100741,0.000944425,0.00133576),0.40366),
sphere((0.00108155,0.000966389,0.00135897),0.40366),
sphere((0.0011069,0.00102797,0.00131333),0.40366),
sphere((0.00102898,0.00102219,0.00133366),0.40366),
sphere((0.00107495,0.00108015,0.00136599),0.40366),
sphere((0.00112847,0.00110573,0.00131123),0.40366),
sphere((0.00113885,0.000975778,0.00126066),0.40366),
sphere((0.00106093,0.000970005,0.001281),0.40366),
sphere((0.00105433,0.00108377,0.00128802),0.40366),
Re: [Yade-users] [Question #695065]: ThermoHydro analyses
Question #695065 on Yade changed: https://answers.launchpad.net/yade/+question/695065 Robert Caulk requested more information: >The compacted sample can be found here: External links are explicitly forbidden from the forum. An MWE is self contained in a single script. Please review our rules at [1]. Thank you! [1]https://www.yade-dem.org/wiki/Howtoask -- 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 #695065]: ThermoHydro analyses
Question #695065 on Yade changed: https://answers.launchpad.net/yade/+question/695065 Status: Open => Needs information Robert Caulk requested more information: Hello, Which example script are you referring to? You just point to a commit of many scripts... Please link to the exact script somewhere in the examples folder [1]. >When I use a previously compacted, poly-disperse-particle sample for this analysis, the temperature of all particles remains constant at first and then turns to "nan". What exactly are you changing compared to the example script? You are using a poly-disperse packing instead of a mono-disperse packing? Please try to say something like "when I change *BLANK* the script works, when I change *BLANK* the script fails." -Robert [1] https://gitlab.com/yade- dev/trunk/-/tree/master/examples/ThermalEngine -- 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
[Yade-users] [Question #695065]: ThermoHydro analyses
New question #695065 on Yade:
https://answers.launchpad.net/yade/+question/695065
Hello,
I am trying to run a TH analysis with warm fluid going through a cold pack of
polydisperse particles using the provided example in [1].
When I use a previously compacted, poly-disperse-particle sample for this
analysis, the temperature of all particles remains constant at first and then
turns to "nan". I was wondering what needs to be adjusted in the thermal and
flow engines following the features of the granular sample under study. The
compacted sample can be found here:
https://drive.google.com/file/d/1Bg-TxgnM4DNe_7t8_QOghIxuspoZiY5T/view?usp=sharing
[1]
https://gitlab.com/yade-dev/trunk/-/commit/f0c96c63937a666d1fa1529173b554ac74b04456
The MWE follows:
Thanks - Zoheir
###--TH analyses###
from yade import pack, ymport
from yade import timing
import numpy as np
import shutil
young=1e9
rad=0.003
mn,mx=Vector3(0.0003548178,0.0003386574,0.0007729061),Vector3(0.0014559856,0.0013972314,0.00188235628)
# corners of the initial packing
thermalCond = 2. #W/(mK)
heatCap = 710. #J(kg K)
t0 = 333.15 #K
# micro properties
r = rad
k = 2.0 # 2*k*r
Cp = 710.
rho = 2600.
D = 2.*r
m = 4./3.*np.pi*r**2/rho
# macro diffusivity
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 = O.bodies.append(ymport.textExt('ShCube', 'x_y_z_r',color=(0.1,0.1,0.9),
material='spheres'))
ThermalEngine = ThermalEngine(dead=1,label='thermal');
newton=NewtonIntegrator(damping=0.2)
intRadius = 1
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intRadius),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intRadius),Ig2_Box_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()],label="iloop"
),
FlowEngine(dead=1,label="flow",multithread=False),
ThermalEngine,
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
newton
]
for b in O.bodies:
if isinstance(b.shape, Sphere):
b.dynamic=False # mechanically static
flow.dead=0
flow.defTolerance=-1#0.3
flow.meshUpdateInterval=200
flow.useSolver=4
flow.permeabilityFactor= 1
flow.viscosity= 0.001
flow.bndCondIsPressure=[1,1,0,0,0,0]
flow.bndCondValue=[10,0,0,0,0,0]
flow.thermalEngine=True
flow.debug=False
flow.fluidRho = 997
flow.fluidCp = 4181.7
flow.getCHOLMODPerfTimings=True
flow.bndCondIsTemperature=[1,0,0,0,0,0]
flow.thermalEngine=True
flow.thermalBndCondValue=[343.15,0,0,0,0,0]
flow.tZero=t0
flow.pZero=0
#flow.clampKValues=False
flow.maxKdivKmean=1
flow.minKdivmean=0.0001;
thermal.dead=0
thermal.debug=False
thermal.fluidConduction=True
thermal.ignoreFictiousConduction=False#True
#thermal.minimumFluidCondDist=1.89e-5
thermal.conduction=True
thermal.thermoMech=False
thermal.solidThermoMech = False
thermal.fluidThermoMech = False
thermal.advection=True
thermal.bndCondIsTemperature=[0,0,0,0,0,0]
thermal.thermalBndCondValue=[0,0,0,0,0,0]
thermal.fluidK = 0.6069#0.650
thermal.fluidConductionAreaFactor=1.
thermal.particleT0 = t0
thermal.particleDensity=2600.
thermal.particleK = thermalCond
thermal.particleCp = heatCap
thermal.useKernMethod=True
#thermal.useHertzMethod=False
timing.reset()
O.dt=0.1e-5
O.dynDt=False
O.run(1,1)
flow.dead=0
def bodyByPos(x,y,z):
cBody = O.bodies[1]
cDist = Vector3(100,100,100)
for b in O.bodies:
if isinstance(b.shape, Sphere):
dist = b.state.pos - Vector3(x,y,z)
if np.linalg.norm(dist) < np.linalg.norm(cDist):
cDist = dist
cBody = b
print('found closest body ', cBody.id, ' at ', cBody.state.pos)
return cBody
# Finding a body at the X flow boundary
bodyOfInterest = bodyByPos(0.00035,0.001,0.001)
from yade import plot
def history():
print(bodyOfInterest.state.temp)
plot.addData(
ftemp1=flow.getPoreTemperature((0.00035,0.001,0.001)),
t=O.time,
i = O.iter,
bodyOfIntTemp = O.bodies[bodyOfInterest.id].state.temp
)
O.engines=O.engines+[PyRunner(iterPeriod=5,command='history()',label='recorder')]
plot.plots={'t':(('ftemp1','k-'),('bodyOfIntTemp','r-'))} #
plot.plot()
O.saveTmp()
def ColorScaler():
for s in O.bodies:
s.shape.color=scalarOnColorScale(s.state.temp,333,343)
O.engines=O.engines+[PyRunner(command='ColorScaler()',iterPeriod=1)]
ColorScaler()
O.run(1000)
--
You received this question notification because your team yade-users is
an answer
