yade-dev team mailing list archive

[Bruno Chareyre <1362090@xxxxxxxxxxxxxxxxxx>]

> If I uncomment line 145: flow.boundaryUseMaxMin=[0,0,0,0,0,0]
> triangulation is not working correctly...

This is the real problem. When using max/min the pressure is weird but
it is not a bug, it is expected behavior. Basically you can use that
only when the boundaries are fixed or they have imposed pressure.

Thanks for the simplified script. I've seen the problem. You define
walls with thickness 0 while the default in FlowEngine is 0.001.
Thus the positions of flow boundaries are wrong and it breaks triangulation.
Assigning the correct value solves the problem:
flow.wallThickness=0

Clearly, having 0.001 thickness by default is stupid. It should be zero,
this is the bug. I will fix that.

Bruno

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

Title:
  problem with pressure calculation in PFV cells

Status in Yet Another Dynamic Engine:
  New

Bug description:
  Hi,

  I want to get the pore pressure of a saturated soil sample in an
  undrained triaxial test when the sample volume is decreasing.

  ...
  TriaxialStressController(thickness=0,stressMask=7,internalCompaction=False,label='triax')
  ...

  	triax.goal1 = 35000
  	triax.goal2 = 35000
  	triax.goal3 = 70000
  	triax.wall_back_activated = False
  ...
  	flow.viscosity = 1.3e-3
  	flow.fluidBulkModulus = 2e9
  	flow.bndCondValue=[0,0,0,0,0,0]
  	flow.bndCondIsPressure=[0,0,0,0,0,0]
  ...

  For one of my simulations a decrease of ca. 1,8 % in porosity (=-dV/V)
  is obtained. When I calculate corresponding pressure increase I get

  dp = K*(-dV/V) = 2e9*0.018 = 36000000 Pa

  But the output from PFV shows maximum average pressures around of 16
  Pa (not MPa!). As long as I use SI units the output of
  flow.averagePressure() should be Pa, right?

  So is there a problem in the code or did I miss something in the
  compressibility option of PFV?

  regards,

  christian

  
  ###additional info:

  I double-checked my script and did not see anything wrong. So I think
  it is on C++ side. So I was going on testing PFV method:

  I also tested a drained situation with incompressible fluid. The porosity is decreasing due to a triggered collapse of the grain structure, but the pore pressure gets negative (see attachment).
  So I think it is a bug in pressure calculation of PFV method.

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