Error with lubrication law solver
Hello,
I'm trying to simulate a viscous flow using a lubrication law, but I'm running into an issue with two of the available resolution methods. I was initially using a modified version of the latest yadedaily release (2018.02b-
I've simplified the script to try to pinpoint the exact issue:
from yade import pack,plot
import math
import random as rand
import numpy as np
from scipy.stats import gaussian_kde
O.dt=1e-6
slope = math.radians(31)
length = 9
width = 2
height = 2
gravityVector = Vector3(
O.materials.
O.materials.
DepBox = O.bodies.
partCloud = pack.SpherePack()
partCloud.
spIds = partCloud.
theta_method = 0.55
resolution = 2
O.engines=[
ForceResetter(),
InsertionSortC
[
Bo1_Sphere_
Bo1_Box_Aabb(),
Bo1_Wall_Aabb()
],
verletDist=0.1,
allowBiggerTha
InteractionLoop(
[
Ig2_Sphere_
Ig2_Box_
],
[
Ip2_FrictMat_
],
[
Law2_
activateN
activateT
activateT
activateR
theta=
resolutio
warnedOnc
maxSubSte
SolutionT
debug=True)
]),
DragEngine(
NewtonIntegrat
]
Drag.ids=spIds
Drag.Rho=1200.
Drag.Cd=0.47
This runs fine when using resolution=0 (Iterative exact resolution with substepping) in the lubrication law options.
However, when resolution=1 is used (Newton-Rafson dimentionless resolution), the following error fills the console entirely several times over:
WARN /data/trunk/
A different error appears with resolution=2 (Dichotomy dimentionless resolution), filling the console as well:
WARN /data/trunk/
WARN /data/trunk/
ERROR /data/trunk/
Is something wrong with my script?
I'm also not entirely sure I understand the difference between those methods and if I should just use the one that works. Is there a difference in precision?
Thanks,
Léo
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Léo Sénique for more information if necessary.