Custom Diffusivity #27012
Answered
by
cpgr
marinsiebert
asked this question in
Q&A Modules: Porous Flow
Custom Diffusivity
#27012
-
|
Hi, what would be the most straight forward approach of creating a custom diffusivity material that takes the components k diffusion coefficients d in phase B, where the coefficients are expressed as a function of temperature and pressure? Ideally this works as an extension within the Millington Quirk material or as a modification, where the diffusion_coefficient vector accepts other materials or aux_variables as input. Right now Millington Quirk only accepts direct input of float values. Thank you |
Beta Was this translation helpful? Give feedback.
Answered by
cpgr
Mar 19, 2024
Replies: 3 comments 12 replies
-
|
You'll have to write your own material to achieve this, but it shouldn't be too difficult. There are examples of taking variables and using them inside the material, eg, PorousFlowPermeabilityConstFromVar. |
Beta Was this translation helpful? Give feedback.
0 replies
-
|
I tried writing my own Base, so that I could then use the custom diffusion coefficients in the other diffusivity objects. |
Beta Was this translation helpful? Give feedback.
12 replies
-
|
Hey @marinsiebert - if you want I can start this off and you can finish it/ test it? |
Beta Was this translation helpful? Give feedback.
-
|
As much as i like messing with the code, any help is very much appreciated @cpgr After @GiudGiud remark about pointers I was now (still) looking at std::vector<variable_value> with * inside the brackets, this was actually accepted by the compiler, but I am not sure how to populate this now. |
Beta Was this translation helpful? Give feedback.
-
|
OK. I'll make a start on it |
Beta Was this translation helpful? Give feedback.
-
|
Check out https://github.com/cpgr/moose/tree/pfdiffvar - the last two commits have the changes I made. This will enable you to use an AuxVariable in place of the diffusion coefficient, which I assumed is what you wanted. So you should be able to do something like where |
Beta Was this translation helpful? Give feedback.
Answer selected by
marinsiebert
-
|
I see, I will try my best to comprehend the involved code changes for future requirements. Thank you ! |
Beta Was this translation helpful? Give feedback.
-
|
You should be able to check out my branch, compile it and test if it works as you want. Let me know! |
Beta Was this translation helpful? Give feedback.
-
|
Works as intended. Thank you! |
Beta Was this translation helpful? Give feedback.
0 replies
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Check out https://github.com/cpgr/moose/tree/pfdiffvar - the last two commits have the changes I made. This will enable you to use an AuxVariable in place of the diffusion coefficient, which I assumed is what you wanted.
So you should be able to do something like
where
d00etc are aux variables