How to modify the porosity, permeability and wettability of the model in OpenPNM? #2357
-
|
I reproduced many examples in OpenPNM and wanted to extend on the basis of the examples, but so far only the pore diameter and throat length and diameter have been modified. How to modify its porosity, wettability, permeability and other physical parameters of porous media? If we take this as an example, how should we modify it? |
Beta Was this translation helpful? Give feedback.
Replies: 1 comment 14 replies
-
|
Well, by changing pore/throat diameter, you're essentially changing the porosity, and subsequently the permeability. You can't directly change the permeability of a pore network model (or any other pore-scale method), that's what you do in a volume-averaged continuum model. As for wettability, I'm not very familiar, @jgostick knows more. But, I guess it has to do with the phase properties (ex. you can modify contact angle in your phase, say |
Beta Was this translation helpful? Give feedback.
-
|
In principle, you can. Just copy those two functions and paste it in your own topotools. But, it might need some tweaking. Let's try one last time:
If this also doesn't work, I highly recommend that you remove your Anaconda distribution and reinstall it, as there might be some serious conflicts already in your base environment. |
Beta Was this translation helpful? Give feedback.
-
|
Sorry, I still failed. Is there any other way to calculate porosity? It can be rough, such as only calculating the ratio of pore volume to total volume |
Beta Was this translation helpful? Give feedback.
-
import numpy as np
import openpnm as op
shape = [5, 6, 7]
spacing = 1e-5
net = op.network.Cubic(shape, spacing)
geo = op.geometry.SpheresAndCylinders(network=net, pores=net.Ps, throats=net.Ts)
volume_void = net["pore.volume"].sum() + net["throat.volume"].sum()
volume_total = np.prod(shape) * spacing**3
porosity = volume_void / volume_total
print(f"Porosity: {porosity:.3f}) |
Beta Was this translation helpful? Give feedback.
-
|
Sorry, before seeing your reply. I found the definitions of these two functions on the official website of OpenPNM, so I downloaded all the files on the official website and replaced the previous files. As a result, the following problem now occurs. |
Beta Was this translation helpful? Give feedback.
-
|
It's really hard to tell. The only clean option that comes to mind is:
|
Beta Was this translation helpful? Give feedback.
Well, by changing pore/throat diameter, you're essentially changing the porosity, and subsequently the permeability. You can't directly change the permeability of a pore network model (or any other pore-scale method), that's what you do in a volume-averaged continuum model. As for wettability, I'm not very familiar, @jgostick knows more. But, I guess it has to do with the phase properties (ex. you can modify contact angle in your phase, say
op.phase.Mercury).