yade-dev team mailing list archive

[Christian Jakob <1260213@xxxxxxxxxxxxxxxxxx>]

Sorry, I found the bug in my script. Set O.run(1,1) before getCell()
solves these two error messages.

But another error occurs:

ERROR /home/christian/YADE/trunk/pkg/dem/FlowEngine.cpp:189 getCellFlux:
Getting flux with cond higher than imposedP size.

-- 
You received this bug notification because you are a member of Yade
developers, which is subscribed to Yade.
https://bugs.launchpad.net/bugs/1260213

Title:
  error in FlowEngine with imposed pressure and getCellFlux()

Status in Yet Another Dynamic Engine:
  New

Bug description:
  Hi,

  I want to increase pore pressure in a PFV cell in the middle of my model containing loose granular material (sand).
  So I did:
  #################

  porePressure = 5000000

  middleCell = flow.getCell(Vector3(xDim/2,yDim/2,zDim/2)) # [1]

  print 'increased pore pressure at ',Vector3(xDim/2,yDim/2,zDim/2), '
  to ',porePressure

  flow.imposePressure(Vector3(xDim/2,yDim/2,zDim/2),porePressure) # [2]

  #I want to see which velocity corresponds to this porePressure:

  vel_fluid = flow.getCellFlux(middleCell) # [3]

  #################

  This error messages occured:

  Triangulation does not exist. Waht did you do?!
  increased pore pressure at  Vector3(0,001,0,001,0,0015)  to  5000000
  Traceback (most recent call last):
    File "<string>", line 1, in <module>
    File "scripts/definitions-final.py", line 78, in writeData
      vel_fluid = flow.getCellFlux(middleCell)
  OverflowError: can't convert negative value to unsigned

  [1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngine.getCell
  [2] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngine.imposePressure
  [3] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngine.getCellFlux

To manage notifications about this bug go to:
https://bugs.launchpad.net/yade/+bug/1260213/+subscriptions