[Yade-users] [Question #692704]: Using flow engine got the wrong pressure field
New question #692704 on Yade:
https://answers.launchpad.net/yade/+question/692704
Hi,
I'm currently simulating a simple seepage problem. My code is as followed:
from yade import pack, plot, ymport
from numpy import genfromtxt
mn,mx=Vector3(-.25,-.25,0),Vector3(.25,.25,.5)
SoilMat=CohFrictMat(young=20e9,poisson=0.3,density=2650,frictionAngle=0.7,alphaKr=50,alphaKtw=50,momentRotationLaw=True)
O.materials.append((SoilMat))
O.materials.append(CohFrictMat(young=10e9,poisson=0.5,frictionAngle=0,density=0,label='walls',momentRotationLaw=True))
walls=aabbWalls([mn,mx],thickness=0,material='walls')
wallIds=O.bodies.append(walls)
O.bodies.append(ymport.text('/home/jp/DEMSimulation/Sample/SP-TestSample5.txt',material=SoilMat))
newton=NewtonIntegrator(damping=0.3)
#Engine of the simulation
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Box_Aabb(),Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
FlowEngine(dead=1,label="flow"),
newton,
]
#Time step
O.dt=.1*PWaveTimeStep()
O.run(1,1)
newton.damping=0
flow.dead=0
flow.defTolerance=0.3
flow.meshUpdateInterval=200
flow.useSolver=3
flow.permeabilityFactor=1
flow.viscosity=10
flow.bndCondIsPressure=[0,0,0,0,1,1]
flow.bndCondValue=[0,0,0,0,100,0]
flow.boundaryUseMaxMin=[0,0,0,0,0,0]
The problem is that the first iteration gives the right result, like [1], while
the second and the following gives a wrong pressure field. The pore pressure of
some cells may reach -1 and more, like [2].
I think the problem is with the sample. When I change the SP-TestSample5.txt to
the following code:
sp=pack.SpherePack()
sp.makeCloud((-.25,-.25,0),(.25,.25,.5),rMean=.0175,rRelFuzz=.5)
sp.toSimulation(material=SoilMat)
the result is reasonable again.
However, I don't know what's wrong with the sample of SP-TestSample5.txt. This
txt file can be found in [3]. The way I generate the SP-TestSample5.txt can be
found in [4]. Any idea why this happens?
Thank you,
Peilun
[1]https://cloud.tsinghua.edu.cn/f/d076877a726b47808591/
[2]https://cloud.tsinghua.edu.cn/f/ba447f257ddb4caf90e9/
[3]https://cloud.tsinghua.edu.cn/f/7e8ede6246554253aa64/
[4]https://cloud.tsinghua.edu.cn/f/ec0e4790e006474f9c09/
--
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 #692650]: How to use the DFNFlowEngine
Question #692650 on Yade changed: https://answers.launchpad.net/yade/+question/692650 Status: Open => Solved Huang peilun confirmed that the question is solved: Thanks Jérôme Duriez, 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 #692609]: Triaxial test of cylindrical specimens
Question #692609 on Yade changed: https://answers.launchpad.net/yade/+question/692609 chanaka Udaya posted a new comment: Paper Particle-Scale Mechanisms in Undrained Triaxial Compression of Biocemented Sands: Insights from 3D DEM Simulations with Flexible Boundary Link https://ascelibrary.org/doi/full/10.1061/%28ASCE%29GM.1943-5622.0001346?af=R Thesis Discrete Element Modeling of the Mechanical Response of Cemented Granular Materials by: Pu Yang lInk https://repository.asu.edu/attachments/211427/content/Yang_asu_0010E_18510.pdf see page 204 onwards to see the 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 #692332]: Flow around cylinder
Question #692332 on Yade changed: https://answers.launchpad.net/yade/+question/692332 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 #692613]: Order of defining object and engine is causing problems
Question #692613 on Yade changed: https://answers.launchpad.net/yade/+question/692613 Status: Needs information => Open Rohit John gave more information on the question: Dear Robert, Thanks for your reply. I am using two materials: one for the internal nodes and one for the external surface. I tried your recommendation but it did not help. The cylinder is still falling apart. I think we need all the Laws. I followed the code given here (https://gitlab.com/yade-dev/trunk/blob/master/examples/cylinders /cylinderconnection-roots.py). It says the following Law2_ScGeom_FrictPhys_CundallStrack(), # contact law for sphere-sphere Law2_ScGridCoGeom_FrictPhys_CundallStrack(),# contact law for cylinder-sphere Law2_ScGeom6D_CohFrictPhys_CohesionMoment(),# contact law for "internal" cylinder forces Law2_GridCoGridCoGeom_FrictPhys_CundallStrack() # contact law for cylinder-cylinder interaction I tried removing Law2_ScGeom_FrictPhys_CundallStrack() and I got the following error None of given Law2 functors can handle interaction #0+1, types geom:ScGeom=1 and phys:FrictPhys=3 (LawDispatcher::getFunctor2D returned empty functor) QObject::~QObject: Timers cannot be stopped from another thread Segmentation fault (core dumped) I tried removing Law2_ScGridCoGeom_FrictPhys_CundallStrack() that did not do anything. Still the only fix is to define the cylinder after the engines. I don't understand why it works though. -- 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 #692650]: How to use the DFNFlowEngine
Question #692650 on Yade changed: https://answers.launchpad.net/yade/+question/692650 Jérôme Duriez posted a new comment: Hi, In terms of articles, [Papachristos2017] from https://yade- dem.org/doc/publications.html#journal-articles should be a good starting point. I could not find examples in current source code. >From what I know (not much) from these engines, DFNFlowEngine differs from FlowEngine through the (additional ?) consideration of flow through planar fractures (fractures due to broken interactions), instead of a flow in the classical pore space. I do not know if there is a dynamic / quasi-static difference as well. -- 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
