yade-dev team mailing list archive

[Christian Jakob <1260213@xxxxxxxxxxxxxxxxxx>]

Public bug reported:

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

** Affects: yade
     Importance: Undecided
         Status: New

-- 
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