Some questions about saturation states #2760
Unanswered
Tianmeng818
asked this question in
Q&A
Replies: 1 comment 3 replies
-
The software continues to evolve with time, so I am not surprized that old codes does not work with the new versions. We try to collect useful examples on the website...have you searched through all of those already? If you cannot find one that explains what you're doing, then let me know and I'll add it to my to-do list to write one. |
Beta Was this translation helpful? Give feedback.
3 replies
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
Dear Doctor:
Hello! Previously in studying this article of yours from 2016 "OpenPNM: A Pore Network Modeling Package", I have successfully reproduced it in versions 2.8.2 and 3.1.2, at least without reporting errors. However, there have been a few questions that have been bothering me and I would like to get your help.The following questions are based on version 2.8.2:
1.In this 2016 article, the final effect image is not sure how you implemented it in paraview, because the effect image I generated does not match your original image, I can only generate simple water intrusion or gas escape, not combine them, while your original image combines them.
2.In this article in 2016, about the ">8000" in this line of code, I don't quite understand it, after deleting these two lines of code, the program can run, so I don't know its exact meaning, I hope to get your answer.
phys_air [ 'throat.conductance'] = phys_air [ 'throat.diffusive_conductance '] '(OP [ 'throat.inv_Pc ] 〉 8000)
My guess is that it has something to do with the saturation state, and that this code can be applied to implement a model with a non-saturated initial state, but I have been trying to find that the function I want is not possible. My attempt is to change the value ">8000" to extract a different conductivity, and changing 8000 to 12000 will turn part of the phys_air['throat.conductance'] value to 0, thus reporting an error.
Deleting the zero element in phys_air['throat.conductance'], but cannot assign the array a to phys_air['throat.conductance'], will report an error.
So I'm not sure if my guess "the use of this line of code is related to the saturation state of the implementation model" is correct, and would appreciate your answer.Maybe you have another method or example to create a non-saturated initial state?
3.So I tried to use other methods to implement the initial state of non-saturation, based on your 2016 code, I added an attempt to set the initial saturation state by fixing the initial saturation and marking pore.occupancy as 1, indicating that the pore is occupied.
An error is reported after running the code.
So I would like to ask what is the solution for this error report? And is there any example to build a non-saturated initial state to follow?
4.Now I try to implement in openpnm2.8.2 version, first the water body invades from the upper part to the lower part, then the gas invades from the upper part to the lower part, and then after the whole process, finally get the state at a certain moment in the whole process, similar to doing different time slices for the model, I want to try the time step related method, but there is no example so I haven't implemented it yet. I don't know if there are any relevant examples to learn.
Sorry if I am disturbing you, but this has really bothered me for a long time and I hope you can help me sort it out.
Beta Was this translation helpful? Give feedback.
All reactions