Parameter setting problem of heating rock sample

Asked by Ziyu Wang

Hello,
I want to simulate the mechanical properties of rock under high temperature in deep underground,so I tried to use FlowEngine and ThermalEngine.But I am unclear about some of the parameters(Mainly the temperature-related part)
1.I have known that FlowEngine.thermalBndCondValue represents the pore boundary condition and ThermalEngine.thermalBndCondValue represents the particle boundary condition.
My question is :The pore here refers to the space between particles, right? If so, do the pore boundary conditions and particle boundary conditions refer to the state of pores and particles on the boundary (such as near the wall)?
2.I want to simulate the heating process of a sample that is initially at room temperature,so I set the following parameters:

##############
flow.tZero=298.15
flow.bndCondIsTemperature=[1,1,1,1,1,1]
flow.thermalBndCondValue=[573.15,573.15,573.15,573.15,573.15,573.15]
thermal.bndCondIsTemperature=[0,0,0,0,0,0]
thermal.thermalBndCondValue=[0,0,0,0,0,0]
thermal.particleT0=298.15
#############

Is it correct to set the parameters like this for my purpose?If there is something inappropriate, how to modify it.

Best regards.

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Robert Caulk
Solved:
Last query:
Last reply:
Revision history for this message
Robert Caulk (rcaulk) said :
#1

Hello,

>>The pore here refers to the space between particles, right?

We have detailed this very carefully with figures and descriptions in our paper [1]. Yes, it is the space between the particles.

>>if so, do the pore boundary conditions and particle boundary conditions refer to the state of pores and particles on the boundary (such as near the wall)?

Currently, the boundary cells can be visualized using the withBoundaries=True argument of flow.saveVTK(). In paraview they are distinguished using "isFictious" flag.

The boundary particles are the particles incident to those flow cells.

Cheers,

Robert

[1]https://www.sciencedirect.com/science/article/abs/pii/S0045782520304771

Revision history for this message
Ziyu Wang (ziyuwang1) said :
#2

Hello Robert,

Yes,I have read the paper[1] and I've learned about the definition of pore space.But I still don't know the difference between the following two boundary conditions:
flow.thermalBndCondValue
thermal.thermalBndCondValue
What's the difference between them?

Thanks again.

Revision history for this message
Best Robert Caulk (rcaulk) said :
#3

Hello,

[1] indicates the temperature to impose to the cells which comprise the boundary of the triangulation.
[2] indicates the temperature to impose to the particles which comprise the boundary of the particle packing.

Cheers,

Robert

[1]https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngine.thermalBndCondValue
[2]https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngine.thermalBndCondValue

Revision history for this message
Ziyu Wang (ziyuwang1) said :
#4

Thanks Robert Caulk, that solved my question.