[Yade-users] [Question #702448]: 2PFV - update Triangulation
New question #702448 on Yade:
https://answers.launchpad.net/yade/+question/702448
Hi all,
I am investigating how the mesh evolves during the drainage using 2PFV.
I my original problem, I have a significant particle deformation/movement
during the drying process, but tracking "getPoreThroatRadiusList()" in
different instants I got the exactly same values. This also for the CellVolume,
CellPorosity. It seems that the cell is constant even the pack is changing over
time.
Any clue?
Here you have a simplified idea of my original script:
#!/usr/bin/python
# -*- encoding=utf-8 -*-
#*
# Copyright (C) 2010 by Bruno Chareyre *
# bruno.chareyre_at_grenoble-inp.fr *
#*
# This program is free software; it is licensed under the terms of the *
# GNU General Public License v2 or later. See file LICENSE for details. *
#*/
from yade import pack
from yade import bodiesHandling
from yade import export
from yade import utils
from yade import ymport
import math
import numpy
### DEFINING VARIABLES AND MATERIALS ###
# The following 5 lines will be used later for batch execution
nRead=readParamsFromTable(
num_spheres=3000,# number of spheres
compFricDegree = 1, # contact friction during the confining phase
key='_triax_base_', # put you simulation's name here
unknownOk=True
)
from yade.params import table
num_spheres=table.num_spheres# number of spheres
key=table.key
targetPorosity = 0.55 #the porosity we want for the packing
compFricDegree = table.compFricDegree # initial contact friction during the
confining phase (will be decreased during the REFD compaction process)
finalFricDegree = 30 # contact friction during the deviatoric loading
rate=0 # loading rate (strain rate)
damp=0.8 # damping coefficient
stabilityThreshold=0.01 # we test unbalancedForce against this value in
different loops (see below)
#2e4+70e4medio 1e4+70e4bom 1e4+60e4bom 3e4+90e4+w3,1,-1-the best
young=80e5 # contact stiffness200e4
young2=80e5
youngcoat=80e5
bondstr=1e3#2e7
bondstr2=1e3
bondstrcoat=1e6
## create materials for spheres and plates
mat=O.materials.append(JCFpmMat(type=1,young=young,poisson=0.3,frictionAngle=radians(compFricDegree),density=2000,tensileStrength=bondstr,cohesion=bondstr,jointNormalStiffness=0,jointShearStiffness=0,jointCohesion=bondstr,jointFrictionAngle=radians(0),jointDilationAngle=0.0,label='spheres'))
O.materials.append(JCFpmMat(type=1,young=20e7,poisson=0.3,frictionAngle=radians(0),density=2600,tensileStrength=0,cohesion=0,jointNormalStiffness=0,jointShearStiffness=0,jointCohesion=0,jointFrictionAngle=radians(0),jointDilationAngle=0.0,label='walls'))
O.materials.append(JCFpmMat(type=1,young=youngcoat,poisson=0.3,frictionAngle=radians(1),density=1500,tensileStrength=bondstrcoat,cohesion=bondstrcoat,jointNormalStiffness=0,jointShearStiffness=0,jointCohesion=bondstrcoat,jointFrictionAngle=radians(0),jointDilationAngle=0.0,label='spherescoat'))
## create walls around the packing
mn,mx=Vector3(0,0,0),Vector3(0.0015,0.0015,0.0015)
mnbox,mxbox=Vector3(-0.0003,-0.0002,0.0003),Vector3(0.0025,0.0025,0.0025)#Vector3(0.002,0.00195,0.002)
walls=aabbWalls([mnbox,mxbox],thickness=0,material='walls')
wallIds=O.bodies.append(walls)
sp=pack.SpherePack()
sp.makeCloud(minCorner=mnbox, maxCorner=mxbox,rMean=0.0001)
O.bodies.append([sphere(center,rad,material='spheres') for center,rad in sp])
#O.bodies.append(ymport.textExt("matrix_vtest6.txt", format='x_y_z_r',
shift=Vector3(0,0.0001,0), scale=1.0,material='spheres',color=(0,1,1)))#0.5
#O.bodies.append(ymport.textExt("coat_vtest6r756.txt", format='x_y_z_r',
shift=Vector3(0,0,0), scale=1.0,material='spherescoat',color=(0,1,1)))
#O.bodies.append(ymport.textExt("coat_vtest6r56.txt", format='x_y_z_r',
shift=Vector3(0,0,0), scale=1.0,material='spherescoat',color=(0,1,1)))
###
## DEFINING ENGINES ###
###
triax=TriaxialStressController(
## TriaxialStressController will be used to control stress and strain. It
controls particles size and plates positions.
## this control of boundary conditions was used for instance in
http://dx.doi.org/10.1016/j.ijengsci.2008.07.002
maxMultiplier=1.+2e4/young, # spheres growing factor (fast growth)
finalMaxMultiplier=1.+2e3/young, # spheres growing factor (slow growth)
thickness = 0,
## switch stress/strain control using a bitmask. What is a bitmask, huh?!
## Say x=1 if stess is controlled on x, else x=0. Same for for y and z,
which are 1 or 0.
## Then an integer uniquely defining the combination of all these tests is:
mask = x*1 + y*2 + z*4
## to put it
Re: [Yade-users] [Question #702425]: 2PFV - Cell Velocity
Question #702425 on Yade changed: https://answers.launchpad.net/yade/+question/702425 Status: Answered => Solved Luis Barbosa confirmed that the question is solved: Thanks Robert Caulk, 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 #702276]: Adding an attribute to a material
Question #702276 on Yade changed: https://answers.launchpad.net/yade/+question/702276 Status: Open => Solved Robert Caulk changed the question status: User solved it 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
[Yade-users] [Question #702446]: liquid migration
New question #702446 on Yade: https://answers.launchpad.net/yade/+question/702446 Hello, I am using the liquid migration feature in Yade and have a couple of questions: 1) in order to use this model only LiqControl should be added to the engine before the newton integrator, is this right? In addition, I wanted to check that the liquid bridge volume Vb will be calculated based on the given Vf and Vmin. So, we should determine Vf and Vmin for each particle. Then the quantity initialised for Vb must not be important. Is this correct ? 2) If I want to track the liquid content amount of each of the particles in each timestep, in which parameter is it stored ? Vf? Thank you in advance! Regards, Roxana Saghafian -- 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 #702418]: About flow.bndCondIsPressure and flow.bndCondValue in flowengine
Question #702418 on Yade changed: https://answers.launchpad.net/yade/+question/702418 Status: Answered => Solved Ziyu Wang confirmed that the question is solved: Thanks Robert Caulk, 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 #702331]: CHOLMOD test with GPU acceleration
Question #702331 on Yade changed: https://answers.launchpad.net/yade/+question/702331 Status: Answered => Solved Ziyu Wang confirmed that the question is solved: Thanks Robert,I will study this problem further. -- 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 #702331]: CHOLMOD test with GPU acceleration
Question #702331 on Yade changed: https://answers.launchpad.net/yade/+question/702331 Ziyu Wang confirmed that the question is solved: Thanks Robert Caulk, 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 #702421]: particle/particle Conductive heat transfer in thermal engine
Question #702421 on Yade changed: https://answers.launchpad.net/yade/+question/702421 Status: Answered => Open Rioual is still having a problem: Hello Robert, Thanks for your input. Do you have references concerning the Kern method and Hertz method ? Thanks, Fr. * François Rioual Chargé de Recherches CR1 Unité FRISE- Génie des procédés frigorifiques pour la sécurité alimentaire et l’environnement 1 rue Pierre- Gilles de Gennes - CS10030 92761 Antony Cedex Courriel: [ mailto:francois.rio...@inrae.fr | francois.rio...@inrae.fr ] Tél. : +33 (0)1 40 96 61 21 (poste : 81 90 15) *** - Mail original - De: "Robert Caulk" À: "Rioual François" Envoyé: Lundi 11 Juillet 2022 11:06:01 Objet: Re: [Question #702421]: particle/particle Conductive heat transfer in thermal engine Your question #702421 on Yade changed: https://answers.launchpad.net/yade/+question/702421 Status: Open => Answered Robert Caulk proposed the following answer: Hello, That is the thermal resistance using a cross sectional area based on the contact area between two spherical particles using the Kern method. If you prefer, you can use the hertz method with thermal.useKernMethod = False thermal.useHertzMethod = True Or set them both false to consider a wall consuming the space between the particle centers (better for continuum modeling). You can easily add in which ever relationship you want there, such as the method you refer to in CFD. Cheers, robert -- If this answers your question, please go to the following page to let us know that it is solved: https://answers.launchpad.net/yade/+question/702421/+confirm?answer_id=0 If you still need help, you can reply to this email or go to the following page to enter your feedback: https://answers.launchpad.net/yade/+question/702421 You received this question notification because you asked the 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 #702331]: CHOLMOD test with GPU acceleration
Question #702331 on Yade changed: https://answers.launchpad.net/yade/+question/702331 Status: Open => Answered Robert Caulk proposed the following answer: Hello, >>But I have a workstation with nvidia tesla T4,it seems suitable for GPU acceleration? Check the double precision on your CPU and your tesla. Which ever is higher is going to be faster. Cheers, Robert -- 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 #702418]: About flow.bndCondIsPressure and flow.bndCondValue in flowengine
Question #702418 on Yade changed: https://answers.launchpad.net/yade/+question/702418 Status: Open => Answered Robert Caulk proposed the following answer: >>If is suitable? Yes this is how you apply a pressure gradient across the specimen which will allow you to compute permeability. Cheers, Robert -- 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 #702425]: 2PFV - Cell Velocity
Question #702425 on Yade changed: https://answers.launchpad.net/yade/+question/702425 Status: Open => Answered Robert Caulk proposed the following answer: Hello, You need to provide the coordinates of the cell to this function [1]: unsat.getCellVelocity((1,1,1)) To get velocity at the location of (1,1,1). Cheers, Robert [1]https://yade- dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngine.getCellVelocity -- 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 #702421]: particle/particle Conductive heat transfer in thermal engine
Question #702421 on Yade changed: https://answers.launchpad.net/yade/+question/702421 Status: Open => Answered Robert Caulk proposed the following answer: Hello, That is the thermal resistance using a cross sectional area based on the contact area between two spherical particles using the Kern method. If you prefer, you can use the hertz method with thermal.useKernMethod = False thermal.useHertzMethod = True Or set them both false to consider a wall consuming the space between the particle centers (better for continuum modeling). You can easily add in which ever relationship you want there, such as the method you refer to in CFD. Cheers, robert -- 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
